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1 .  INTRODUCTION 


I . 1  Background 


This  report  documents  the  results  of  our  two-year  Phase  II  effort  to 
investigate  the  application  of  data  from  the  proposed  DMSP  microwave  moisture 
sounder  (SSM/T-2)  to  the  inference  of  water  vapor  profiles,  temperature  pro¬ 
files,  cloud,  precipitation,  and  the  properties  of  the  underlying  surface. 
Although  originally  conceived  as  an  enhancement  to  the  SSM/T  temperature 
sounder  to  provide  the  capability  to  obtain  water  vapor  profiles  alone 
utilizing  the  water  vapor  resonance  feature  in  the  vicinity  of  183  GHz,  the 
SSM/T-2  sensor  channel  set  is  also  sensitive  to  surface  emission  and  clouds 
and  precipitation  within  the  f ield-of -view.  This  is  due  to  the  inclusion  of 
"window"  channels  at  91.66  and  150.0  GHz  in  addition  to  channels  near  183  GHz 
when  the  water  vapor  abundance  is  small  enough  to  permit  the  sensor  to  see 
into  middle  and  lower  troposphere. 

The  goal  of  this  research  is  to  provide  a  unified  retrieval  scheme  capa¬ 
ble  of  processing  microwave  and  millimeter  wave  sensor  (i.e.,  sounder  and 
imager)  data  which  extends  the  adopted  operational  algorithm  based  on  a  purely 
statistically  based  retrieval  of  temperature  profiles  from  the  SSM/T-1,  mois¬ 
ture  profiles  from  the  SSM/T-2,  and  surface  and  atmospheric  parameters  such  as 
precipitation  and  integrated  cloud  liquid  water  from  the  SSM/I .  The  approach 
uses  two  concepts  to  increase  retrieval  accuracy.  First,  physical  considera¬ 
tions  based  on  the  application  of  the  radiative  transfer  equation  are  used  to 
quality  control  the  statistically  derived  first  guess  retrievals.  Here  the 
key  question  is,  how  consistent  are  the  statistically  retrieved  profiles  of 
temperature  and  moisture  with  the  observed  brightness  temperatures?  Second, 
data  sources  from  colocated  microwave  mission  sensors  (i.e.,  the  SSM/T-1, 
SSM/T-2,  and  SSM/I)  are  combined  when  necessary  to  more  completely  character¬ 
ize  the  spectral  dependence  of  the  observed  brightness  temperatures  on  the 
meteorological  properties  of  the  surface  and  atmosphere.  By  combining  data 
from  all  available  sensors,  the  desired  meteorological  parameters  can  be  more 
accurately  determined.  Based  on  our  evaluation  of  cloud  effects  on  millimeter 
wave  brightness  temperatures  and  water  vapor  retrievals  (Isaacs  and  Deblonde, 
1987),  we  have  suggested  use  of  the  visible  and  infrared  Operational  Linescan 
System  (OLS)  to  provide  information  on  rloud  presence.  The  unified  retrieval 
scheme  has  been  summarized  in  the  paper  by  Isaacs  (1988). 
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In  the  Initial  Phase  I  effort  of  this  program  (Isaacs  et  al ,  1985a), 
simulation  models  were  developed  to  treat  the  instrumental  response  of  a 
sounder  operating  in  the  vicinity  of  183  GHz  to  representative  variations  of 
cloud  and  precipitation  realistically  encountered  in  the  atmosphere.  The 
Phase  I  technical  objectives  are  summarized  in  Table  11.  The  simulation 
model  developed  used  the  gaseous  absorption  algorithm  and  model  atmospheres 
from  RADTRAN  (Falcone  et  al.,  1982)  as  well  as  the  cloud  models  from  the  AFGL 
AEROSOL  subroutine  (Falcone  et  al.,  1979).  In  the  presence  of  precipitation, 
a  multiple  scattering  calculation  is  required.  The  discrete  ordinate  method 
(DOM)  (Liou,  1973)  was  chosen  to  evaluate  the  equation  of  radiative 
transfer.  The  version  of  the  DOM  chosen  for  implementation  (Liou  et  al., 

1980)  had  never  been  successfully  run  at  AFGL  and  required  ancillary 
millimeter  wave  multiple  scattering  data  for  the  precipitation  drop  size 
distribution  which  was  not  internally  available.  The  code  was  successfully 
debugged  and  run  on  the  AFGL  Cyber  and  a  parameterization  was  developed  for 
evaluating  the  required  scattering  properties  as  a  function  of  frequency, 
temperature,  drop  size  distribution,  rain  rate,  and  phase  (i.e.,  ice  or  water) 
without  the  necessity  for  on-line  Mie  theory  calculations.  The  model  is 
described  in  the  Phase  I  interim  report  (Isaacs  et  al . ,  1985a). 

As  the  result  of  limited  instrument  and  sensitivity  tests  performed  with 
this  simulation  model  on  candidate  systems,  the  potential  degradation  of  mois¬ 
ture  retrievals  due  to  the  presence  of  clouds  and  precipitation  was  demon¬ 
strated.  Retrievals  were  performed  using  a  statistical  approach  based  on  the 
"D"  matrix  formalism  (Gaut  et  al,  1975)  which  is  conceptually  quite  similar  to 
the  adopted  operational  retrieval  scheme.  Additionally,  the  dependence  of 
water  vapor  profile  retrieval  accuracy  on  the  nature  of  the  underlying  surface 
was  evaluated.  For  example,  retrievals  over  the  ocean  (assuming  calm  seas) 
were  more  accurate  than  those  over  land. 

In  order  to  improve  moisture  retrievals  in  the  presence  of  clouds  and 
additionally  exploit  the  demonstrated  sensitivity  of  the  SSM/T-2  instrument  to 
obtain  cloud  and  precipitation  properties  when  present,  it  was  suggested  that 
a  retrieval  scheme  should  be  sought  which  explicitly  treats  these  factors. 
Furthermore  it  was  suggested  that  the  optimal  retrieval  approach  should  be 
implemented  utilizing  all  available  DMSP  sensor  data.  Brightness  temperature 
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data  from  the  ensemble  of  microwave  mission  sensors  colocated  aboard  the  DMSP 
spacecraft  (i.e.,  the  SSM/T-1,  SSM/T-2,  and  SSM/I)  and  as  necessary,  visible 
and/or  infrared  radiance  data  from  the  Operational  Linescan  System  (01. S) 
should  be  integrated  to  provide  a  generalized  geophysical  parameter  retrieval. 
This  multispectral  approach  would  enable  available  information  on  surface 
properties  from  the  microwave  imager  to  be  integrated  into  the  retrievals  of 
temperature  and  moisture  profiles.  Additionally,  information  on  the  tempera¬ 
ture  dependence  of  the  water  vapor  profiles  and  water  vapor  dependence  of  the 
temperature  profiles  would  be  included  explicitly. 

This  approach  to  the  retrieval  problem,  i.e.,  as  a  multispectral/multi - 
parameter  system,  expresses  an  emergent  philosophy  within  the  state-of-the-art 
remote  sensing  community  (cf.  Susskind  et  al . ,  1984,  1985;  Smith  et  al., 

1985) . 

1 . 2  Study  Objectives 

The  technical  objectives  of  the  Phase  II  effort  are  given  in  Table  1-1. 
The  primary  objective  of  this  research  effort  was  to  develop  a  modified 
approach  to  the  retrieval  of  water  vapor  profiles  and  cloud  and  precipi tation 
properties,  employing  an  enhanced  version  of  the  millimeter  wave  moisture 
sounder  simulation  model  developed  under  the  Phase  I  effort.  As  described 
above,  this  approach  will  integrate  data  from  the  microwave  mission  sensors 
and  apply  physical  retrieval  considerations  to  enhance  the  accuracy  of  water 
vapor  and  temperature  profile  retrievals  and,  additionally,  provide  informa¬ 
tion  on  the  state  of  the  underlying  surface.  Four  subtasks  were  identified  to 
accomplish  the  retrieval  development  objective.  These  include:  (a)  imple¬ 
mentation  of  a  physical  retrieval  method  for  the  183  GHz  moisture  channels, 

(b)  development  of  a  cloud  filtering  algorithm  for  use  in  the  physical 
retrieval,  (c)  incorporation  of  microwave  imager  data  and  other  relevant  data 
(such  as  that  from  the  OLS),  and  (d)  testing,  tuning,  and  verification  of  the 
retrieval  algorithm. 

About  nine  months  into  the  performance  period  covered  by  this  report,  the 
contract  supporting  this  study  was  modified  to  increase  the  level  of  effort 
devoted  to  the  incorporation  of  other  data  within  the  retrieval  development 
(item  c  above).  Specifically,  it  was  desired  to  investigate  the  effects  of 
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the  variability  of  the  microwave  properties  of  the  earth's  surface  and 
atmosphere  within  the  relatively  large  footprint  of  the  SSM/I  .  In  analyzing 
the  SSM/I  f ield-of -view,  the  small  pixel  size  of  the  visible  and  infrared 
imagery  was  to  be  employed  to  provide  guidelines  on  when  the  SSM/1  microwave 
imager  was  viewing  nonuniform  conditions  and  what  the  contributions  of  these 
nonuniform  surfaces  is  to  the  microwave  observations.  The  study  was  to 
include  the  identification  of  rain  areas  and  include  modelling  various  types 
of  surfaces  and  estimating  their  emissi vities .  The  analysis  performed  was  to 
be  supported  by  using  Scanning  Multichannel  Microwave  Radiometer  (SMMR)  data 
along  with  available  visible  and  infrared  imagery. 

In  order  to  successfully  accomplish  the  primary  objective  described 
above,  certain  simulation  model  enhancement  objectives  were  also  identified: 

•  To  improve  the  accuracy  and  efficiency  of  the  multiple  scattering 
submodel  based  on  the  DOM  to  treat  precipitation,  it  was  desired  to 
replace  the  iterative  root  searching  eigenvalue  solution  provided  by 
Liou  et  al.  (1980)  with  a  matrix  eigensystem  approach. 

•  Computer  models  were  also  desired  to  calculate  the  exact  optical 
parameters  required  by  the  multiple  scattering  code  via  Mie  theory 
for  arbitrary  size  distributions  of  glaciated  or  liquid  hydromete¬ 
ors  . 

•  The  existing  DOM  treated  only  nadir  viewing  sensors  and  did  not 
treat  polarization.  A  complete  polarized  multiple  scattering  model 
based  on  the  vectorized  radiative  transfer  equation  was  desired. 

•  The  simulation  model  required  surface  emissivities  as  an  input 
parameter.  It  was  desired  to  develop  physical  models  of  surface 
emission  for  a  variety  of  terrain  surfaces  including  snow  and  sea 
ice . 

This  has  been  a  fruitful  research  program  with  much  interaction  with  AFGL 
and  other  interested  parties.  Much  of  the  work  described  here  has  been 
presented  at  conferences  (Appendix  A)  or  published  (Appendix  B) . 

1 . 3  Overview  of  Report 

This  report  is  divided  into  four  major  sections.  In  Section  2  to  follow, 
the  development  of  a  physically  based  retrieval  approach  using  a  statistical 
first  guess  is  described.  This  algorithm  is  capable  of  treating  brightness 
temperature  data  from  a  variety  of  instruments  and  retrieves  simultaneously  to 


1 


obtain  temperature  profile,  water  vapor  profile,  and  surface  properties  As 
yet,  the  algorithm  treats  clear  skies  only.  Our  work  on  the  effect  of  clouds 
on  millimeter  wave  retrievals  is  also  described  here  including  results  of 
calculations  showing  the  sensitivity  of  SSM/T-2  channel  brightness  tempera¬ 
tures  to  cloud  integrated  water  amounts.  Finally,  we  outline  procedures  to 
use  OLS  imager  data  in  a  unified  retrieval  with  the  microwave  mission  sen¬ 
sors.  A  flowchart  of  the  procedure  is  provided  which  is  a  simple  extension  of 
the  clear  sky  algorithm  described  which  employs  image  processing  of  the  OLS 
data . 

In  Section  3,  the  simulation  model  enhancements  are  described.  These 
topics  include  the  implementation  of  the  eigensystem  solution  for  the  DOM 
multiple  scattering  code,  modification  of  Mie  scattering  codes  for  use  at  AFGL 
which  will  provide  scattering  properties  for  arbitrary  size  distributions  of 
hydrometeors,  development  of  a  vectorized  radiative  transfer  model  to  treat 
the  anisotropic  multiple  scattering  of  microwaves  by  precipitation  in  noniso- 
thermal  atmospheres,  and  the  development  of  a  modularized  set  of  microwave 
surface  models  for  use  in  the  simulation  studies. 

Section  4  discusses  the  work  to  date  on  the  characterization  of  the 
uniformity  of  SSM/I  fields-of -view  based  on  the  analysis  of  high  resolution 
visible  and  infrared  imagery  data. 

2.  TASK  1  -  CLEAR  SKY,  MULTIPLE  INSTRUMENT  RETRIEVAL  ALGORITHM  DEVELOPMENT 

2 . 1  Overview 

The  significance  of  improved  capabilities  to  specify  the  global  moisture 
field  (including  water  vapor,  cloud,  and  precipitation)  potentially  afforded 
by  an  optimally  designed  microwave  moisture  sounder  cannot  be  underestimated. 
Implications  for  numerical  weather  prediction  are  particularly  promising 
(Kaplan  et  al,  1983;  Isaacs  et  al .  ,  1986).  Radiatively,  water  in  its  many 
forms  is  the  most  active  constituent  of  the  troposphere.  Clouds  and  water 
vapor  profoundly  affect  the  infrared  temperature  retrieval  problem  (Chahine, 
1974,  1982)  and  the  parameterization  of  radiative  heating  (WMO,  1978).  A 
variety  of  studies  have  shown  that  updating  the  moisture  field  can  improve 
short  range  precipitation  forecasts  (Perkey,  1980;  Maddox  et  al . ,  1981)  and 
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that  latent  heat  release  associated  with  condensation  can  have  large  short 
term  influence  on  storm  development  (Anthes  et  al.,  1983).  Knowledge  of  the 
moisture  field  is  also  of  great  importance  for  meeting  the  operational 
requirements  for  better  cloud  forecasts  and  predicting  optical  and  millimeter- 
wave  propagation. 

An  observing  system  to  retrieve  moisture  consists  of  both  the  sensor 
hardware  to  obtain  the  data  and  the  software  algorithms  to  retrieve  the 
desired  profiles.  Operational  moisture  sounders  aboard  polar  orbiters  have 
been  infrared  instruments  such  as  the  HIRS-2  component  of  the  TIROS-N 
Operational  Vertical  Sounder  (TOVS)  (Smith  et  al .  ,  1979)  and  the  Defense 
Military  Satellite  Program  (DMSP)  SSH/2  (Barnes  Engineering,  1978).  These 
instruments  have  exploited  the  infrared  rotati on-vibration  bands  near  6.7  ^m 
(three  channels)  and  the  pure  rotational  bands  near  20  pm  (seven  channels) , 
respectively.  In  general,  the  success  in  obtaining  operationally  useful  water 
vapor  profiles  from  these  instruments  has  been  limited.  Goals  for  the  TOVS 
retrievals  using  a  statistical  inversion  method  are  limited  to  precipitable 
water  values  in  three  broad  vertical  layers  at  accuracies  of  30%  (NOAA,  1981) . 
Investigators  using  retrievals  based  on  physical  approaches  have  fared  some¬ 
what  better.  For  example,  Reuter  and  Susskind  (1986)  show  that  an  iterative 
approach  can  improve  over  the  first  guess  based  on  a  forecast  field, 
particularly  in  the  upper  levels  (i.e,,  300-500  mb). 

Microwave  water  vapor  remote  sensing  in  the  past  has  focused  on  the  use 
of  the  weak  rotational  line  at  22.235  GHz  for  obtaining  precipitable  water. 
This  line  does  not  have  sufficient  opacity  to  obtain  a  water  vapor  profile.  A 
much  stronger  feature  available  for  profiling  is  the  ( 3 i 3  -  220)  millimeter 
wave  rotational  feature  at  183.31  GHz  (Schaerer  and  Wilheit,  1979).  A  set  of 
channels  has  been  proposed  for  the  1990' s  N0AA-NEXT  series  called  the  Advanced 
Microwave  Sounding  Unit  B  package  (AMSU-B)  (Schultz,  1982).  A  similar  sensor 
called  the  SSM/T-2  will  enhance  the  SSM/T  microwave  temperature  sounder  aboard 
DMSP.  A  number  of  studies  have  investigated  the  simulated  retrieval  of  mois¬ 
ture  profiles  using  a  hypothetical  183  GHz  sensor  (Schaerer  and  Wilheit,  1979; 
Rosenkranz  et  al ,  1982;  Wang  et  al.,  1983;  Kakar,  1983;  Kakar  and  Lambrigsten, 
1984) .  For  the  most  part  these  studies  have  focused  on  retrievals  in  clear 
sky  conditions.  Recently,  the  effects  of  clouds  on  statistical  retrievals  of 
water  vapor  at  183  GHz  have  been  evaluated  (Isaacs  et  al  1985a, b).  These 


rasa 


f, 


Va'.V 


J 


m 


a 


J 


studies  noted  that  the  presence  of  cloud  can  seriously  degrade  the  accuracy  of 
the  retrieved  profiles.  Furthermore,  it  was  noted  that  surface  emissivity 
plays  a  critical  role,  i.e.,  retrievals  over  emissive  land  areas  are  not  as 


accurate  as  those  obtained  over  the  less  emissive  oceans. 


The  proposed  method  to  retrieve  water  vapor  from  the  DMSP  SSM/T-2  will 
consist  of  a  statistical  approach  based  on  the  "D"  matrix  methodology  to 
obtain  water  vapor  abundances  and  relative  humidity  at  each  of  the  mandatory 
levels  (A.  Stogryn,  personal  communication).  The  proposed  "D"  matrix  will  be 
globally  applicable  and  will  treat  both  clear  and  cloudy  f ields-of -view.  A 
two  stage  retrieval  is  suggested,  one  applying  to  low  and  moderate  water  vapor 
amounts,  and  another  to  large  water  vapor  amounts.  The  primary  justification 
for  the  approach  is  computational  efficiency  (i.e.,  the  "D"  matrix  retrieval 
uses  little  computer  time  and  the  current  operational  unavailability  of  data 
to  support  a  cloud/no  cloud  discrimination  step  within  the  retrieval).  The 
approach  will  use  both  SSM/T-1  and  SSM/T-2  data,  but  no  other  supporting  data. 


For  a  variety  of  reasons,  the  proposed  operational  approach  described 
above  could  be  improved.  The  two  main  criticisms  are  its  reliance  on 
statistical  data  and  the  neglect  of  other  potentially  available  data  sources 
Statistical  methods  provide  excellent  retrievals  of  conditions  near  that  of 
the  mean,  but  biased  retrievals  of  extreme  conditions.  Mean  conditions  are 
generally  of  less  meteorological  interest  and  certainly  of  less  operational 
interest  than  those  which  are  away  from  the  mean.  Furthermore,  physical 
retrievals  have  been  shown  to  outperform  statistical  retrievals  of  temperature 
(Halem  and  Susskind,  1985;  Susskind  and  Reuter,  1985).  The  "D"  matrix  ap¬ 
proach  is  an  inherently  linear  treatment,  while  nonlinearities  can  be  directly 
treated  by  physical  approaches.  To  zeroth  order,  the  temperature  retrieval 
problem  is  linear.  The  water  vapor  problem  is  inherently  nonlinear,  and 
therefore,  physical  methods  should  be  even  more  useful. 


It  would  also  seem  that  the  use  of  a  global,  clear/cloudy  "D"  matrix 
would  degrade  the  accuracy  of  the  water  vapor  retrievals  when  clouds  are  not 
present.  When  clouds  are  present,  this  approach  will  not  improve  the  retriev¬ 
als.  Therefore,  results  are  suboptimal  at  all  times.  Physical  approaches 
should  be  able  to  treat  clouds  explicitly  as  part  of  the  inversion  process. 

The  same  applies  to  the  treatment  of  surface  emissivity.  Eventually,  cloud 
fields  from  the  DMSP  OLS  (using  an  appropriately  spatially  averaged  subset  of 
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visible  and  infrared  imagery)  and/or  over  the  oceans  cloud  liquid  water  data 
from  the  SSM/I  microwave  imager  should  be  available  to  aid  in  cloud/no  cloud 
discrimination. 


Finally,  it  is  not  clear  how  to  quality  control  statistical  retrievals  in 
real  time.  With  a  global  "D"  matrix  all  water  vapor  retrievals  must  he  con¬ 
sidered  to  have  the  same  accuracy.  With  the  physical  retrievals,  the  residu¬ 
als  which  are  evaluated  as  part  of  the  retrieval  process  (i.e.,  the  differ¬ 
ences  between  the  calculated  brightness  temperatures  and  the  observed  bright¬ 
ness  temperatures)  can  be  used  to  flag  potentially  bad  retrievals  (and  good 
ones)  for  the  user.  For  example,  when  residuals  are  small  and  confidence  in 
the  retrievals  is  high  the  retrieved  profiles  can  be  given  a  higher  priority 
in  data  assimilation  schemes  such  as  those  interfacing  between  input  fields 
and  numerical  weather  prediction  models. 

Based  on  the  considerations  outlined  above,  this  study  has  investigated 
an  alternative  retrieval  approach  for  the  SSM/T-2  and  SSM/T-1  which  both  em¬ 
ploys  physical  considerations  and  allows  for  the  incorporation  of  other  data 
sources  such  as  that  from  the  SSM/I .  The  approach  builds  on  the  proposed 
statistical  retrieval  scheme  by  using  "D"  matrix  retrievals  as  first  guess 
profiles  (Section  2.2).  However,  brightness  temperatures  are  then  evaluated 
to  examine  the  consistency  of  the  first  guesses  with  the  observations.  When 
the  residuals  are  small,  the  process  terminates.  However,  when  residuals  are 
larger  than  a  preset  tolerance,  the  procedure  goes  on  to  adjust  the  first 
guess  profiles  using  the  procedure  described  in  Section  2.3.  Results  for 
clear  sky  situations  are  presented  in  Section  2.4.  Our  simulations  suggest 
that  cloud  (and  precipitation)  will  have  an  effect  on  millimeter  wave 
brightness  temperatures  (Section  2.5).  A  simple  way  of  using  the  code  in 
partly  cloudy  cases  would  be  to  calculate  equivalent  clear  column  brightness 
temperatures  from  the  observed  data  set  using  visible  or  infrared  data  to 
identify  the  cloud  coverage  fraction.  This  could  be  used  as  a  first  guess  to 
work  the  cloud  problem  as  an  integral  part  of  the  retrieval  process.  This 
approach  has  been  formulated  and  an  overall  unified  retrieval  incorporating 
all  data  sources  is  proposed  in  Section  2.6. 


9 


2 . 2  Statistical  Retrieval  Results 

As  a  first  step  we  have  investigated  application  ol  statistical  re¬ 
trievals.  In  our  earlier  work  we  showed  application  ol  this  approach  (see 
Section  2.5.3  for  a  description)  to  the  retrieval  of  water  vapor  profiles 
using  simulated  SSM/T-2  data  alone  (Isaacs  and  Deblonde,  1985b).  Here  we 
apply  the  technique  to  a  unified  retrieval  using  all  the  microwave  mission 
sensors . 

This  section  illustrates  the  application  of  the  statistical  retrieval 
approach  to  assessment  ot  the  accuracy  of  temperature  profile  (T) ,  integrated 
’  ater  vapor  profile  (U) ,  surface  temperature  (SURF  T) ,  and  emissivity  (EM1SS) 
retrieval  using  data  simulated  for  the  DMSP  millimeter/microwave  sensor 
suite.  In  this  example  individual  sensor  characteristics  such  as  noise  and 
number  and  location  of  channels  are  kept  fixed  and  the  specific  choice  of 
instruments  included  in  the  retrieval  is  varied.  The  sensor  payload  of  the 
Defense  Meteorological  Satellite  Program  (DMSP)  spacecraft  of  the  1990' s  will 
consist  of  a  visible/infrared  imager  (the  operational  linescan  system  or  OLS) , 
a  microwave  temperature  sounder  (SSM/T-1),  a  millimeter  wave  water  vapor 
sounder  (SSM/T-2),  and  a  microwave  imager  (SSM/I).  The  current  OLS  imager 
provides  high  spatial  resolution,  global  cloud  imagery.  Notably,  all  the 
other  sensors  are  millimeter/microwave  instruments.  Of  these  microwave  mis¬ 
sion  sensors,  the  SSM/T-1  and  SSM/I  are  currently  operational  (Falcone  and 
Isaacs,  1987).  The  SSM/T-2  is  scheduled  for  launch  in  the  early  1990s. 
Instrumental  characteristics  for  these  sensors  including  the  noise  equivalent 
brightness  temperatures  used  in  the  data  simulations  are  given  in  Table  2-1. 
The  use  of  one  of  these  instruments,  combinations  of  two  of  them,  or  the  use 
of  data  from  all  three  were  examined  as  possible  approaches  to  obtain  the  most 
accurate  results.  Data  was  simulated  for  these  sensors  based  on  the  RADTRAN 
(Falcone  et  al.,  1982)  computer  algorithm  and  using  ensembles  of  midlatitude 
and  tropical  atmospheres  from  a  radiosonde  data  set  consisting  of  a  total  of 
400  soundings.  Emissivity  values  were  calculated  from  an  ocean  surface 
reflectance  model  assuming  calm  seas.  The  resultant  values  were  varied 
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Table  2-1.  Defense  Meteorological  Satellite  Program  (DMSP)  Sensors 


assuming  a  Gaussian  distribution  with  the  calculated  value  as  mean  and  an  as¬ 
sumed  standard  deviation.  While  surface  emissivity  per  se  is  not  a  required 
geophysical  parameter,  it  is  emissivity  which  provides  information  on  surface 
winds  and  sea  ice  over  the  ocean,  and  vegetation,  soil  moisture,  and  snow 
cover  over  the  land.  The  surface  temperature  was  taken  as  the  same  as  the 
1000  mb  atmospheric  temperature.  The  simulations  assumed  cloud  free  condi¬ 
tions  over  the  ocean. 

The  retrieval  results  for  various  instrument  combinations  are  shown  in 
Tables  2-2  and  2-3,  for  midlatitude  and  tropical  cases,  respectively.  Tabu¬ 
lated  are  fraction  of  unexplained  variance  (FUV)  values  for  each  instrument 
combination.  Perfect  retrievals  have  an  FJV  of  0.0  while  those  which  perform 
no  better  than  climatology  have  FUVs  of  1.0.  (FUVs  greater  than  1.0  mean  that 
climatology  is  better  than  the  retrieval.)  The  FUV  for  each  parameter  is 


SR K! 


related  to  the  absolute  accuracy  of  that  parameter  in  engineering  units. 

These  are  directly  comparable  to  the  DMSP  data  acquisition  requirements 

These  results  suggest  methods  to  combine  existing  sensor  data  in  order  to 
improve  the  baseline  DMSP  capabilities  to  infer  desired  parameters.  For 
example,  Table  2-2  shows  that  the  combination  of  T  and  T-2  data  (cf  column, 
1&2),  provide  a  more  accurate  integrated  water  vapor  profile  than  the  use  of 
T-2  data  alone  (column  SSM/T2).  Especially  note  the  potential  improvement  in 
temperature  profile  retrieval  near  the  surface  (T  1000)  obtained  when  com¬ 
bining  SSM/T  data  with  that  from  the  SSM/T-2  (column  1S2)  or  SSM/1  (l&I).  The 
importance  of  the  SSM/I  data  combined  with  the  SSM/T  data  for  surface  para¬ 
meter  retrievals  is  also  reflected  in  the  surface  temperature  and  emissivity 
results.  These  midlatitude  results  can  be  compared  to  those  for  the  tropics 
(Table  2-3).  Here,  for  example,  it  can  be  seen  that  the  advantage  to  incor¬ 
porating  SSM/I  data  along  with  SSM/T  data  for  temperature  profile  retrieval 
near  the  surface  (T  1000)  demonstrated  for  midlatitudes  is  lost  in  the  tro¬ 
pics.  This  is  largely  due  to  the  smaller  variance  of  tropical  temperature 
profiles  near  the  surface  compared  to  midlatitude  counterpar ts .  (Compare 
columns  PS.  D. ,  i.e.  parameter  standard  deviation,  for  T  1000  in  Tables  2-2 
and  2-3.)  Note  that  in  the  tropics  the  combination  of  SSM/T  and  SSM/I  data 
still  provide  the  most  accurate  retrievals  of  surface  emissivity.  For  clarity 
we  have  plotted  the  retrieval  results  for  temperature  and  integrated  water 
vapor  profiles  as  RMS  errors  in  degrees  Kelvin  and  g/cm  ,  respectively,  as  a 
function  of  pressure.  These  are  illustrated  in  Figures  2-la,  b  and  2-2a,  b 
for  midlatitude  winter  and  tropical  atmospheres,  respectively. 

2 . 3  Theoretical  Basis  of  Physical  Retrieval  Approach 

Current  operational  analysis  procedures  for  the  OLS,  SSM/T- 1,  and  SSM/I 
treat  each  sensor  data  stream  independent  of  the  others.  OLS  imagery  is 
processed  into  global  cloud  property  fields  using  an  automated  nephanalysis 
algorithm  (Fye,  1978).  The  retrieval  scheme  for  temperature  sounding  by  the 
SSM/T-1  is  based  on  regression  of  SSM/T-1  brightness  temperature  data  against 
desired  mandatory  level  temperatures  (Rigone  and  Stogryn,  1977)  and  an  analo¬ 
gous  statistical  approach  will  be  used  in  the  determination  of  meteorological 
parameters  such  as  cloud  liquid  water  content  and  precipitation  (among  others) 
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Table  2-2.  Fraction  of  Unexplained  Variance  for  Statistical  Retrievals.  ■] 

Midlatitude  Summer  Atmospheres  !' 


1 


PARAM. 

SSM/T1 

SSM/T2 

SSM/I 

1  A  2 

1  &  1 

2  A  1 

1,2  A 1 

P  S  .  0  . 

P  MEAN 

-  i 

T  10 

0.0571 

1.0451 

1.0487 

0.0576 

0.0628 

1.0825 

0.0992 

6 . 7678 

218.9875 

4 

1 

T  50 

0.0588 

1.1534 

1.0212 

0.0608 

0.0632 

1  .  14  72 

0.0718 

6 . 1646 

212.3423 

* 

T  100 

0.0810 

0.9592 

0.9662 

0.0931 

0  1017 

1 .0057 

0. 1149 

5 . 8427 

212.8881 

4 

T  300 

0.2151 

0.7252 

0.7649 

0.2125 

0.2491 

1  .  326') 

0.6011 

4 .2547 

224 .5277 

T  500 

0.1763 

0.5801 

0.7562 

0.2034 

0.1503 

0.8866 

0.4175 

6.1279 

249.7131 

T  700 

0.1475 

0.4791 

0.7577 

0. 1280 

0.1932 

0.4808 

0.2393 

6.2591 

265.9886 

T  850 

0.2597 

0.5396 

0.6479 

0.2818 

0.3243 

0.4993 

0.2925 

5.8369 

273.1358 

\ 

m. 

T  1000 

0.4154 

0.6958 

0.7190 

0.3583 

0.2734 

0.9720 

0.5340 

6.8623 

277 . 7990 

H 

U  300 

0.7515 

0.5441 

0  7247 

0.1068 

0.6382 

0.6703 

0.0993 

0.0026 

0.0068 

U  500 

0.7656 

0.5694 

0. 7163 

0.1478 

0.6339 

0.6601 

0.1329 

0.0400 

0.0662 

U  700 

0.7743 

0.5545 

0.6579 

o.  m9 

0.6053 

0.7806 

0.1489 

0.1859 

0.2948 

J 

U  1000 

0.6667 

0.3476 

0.3022 

0.1998 

0.3653 

0.5747 

0.3294 

0.5531 

1.2548 

SURF  T 

0.4154 

0.6958 

0.7190 

0.3583 

0.2734 

0.9720 

0.5340 

6.8623 

277 . 7990 

J 

EMISS 

0.0269 

0.0319 

0.0237 

0.0192 

0.0109 

0.0368 

0.0225 

0.1232 

0.8041 

■ 

from  the 

SSM/I 

data  (Lo, 

1983)  . 

These 

approache 

s  share 

<3  COJMnor 

heri tag 

o  in 

» 

the  "D" 

matrix 

technique 

described  by  Gaut  et  al 

.  (1975) 

Data 

from  the 

SSM/T-1 

and  SSM/T-2  will 

be  integrated 

together 

in  a  statistical 

retrieval  of 

1 

4 

water  vapor  profiles  (cf.  Isaacs,  1987). 


For  a  variety  of  reasons,  the  operational  approach  described  above  could 
be  improved.  The  main  criticisms  are:  (1)  the  lack  of  a  mul t i spectral  per¬ 
spective,  (2)  reliance  on  statistical  retrieval  approaches,  which  produce 
retrieval  fields  with  reduced  variance  properties,  fail  to  treat  inherent 
problem  nonlinearities,  and  provide  little  opportunity  to  monitor  retrieval 
quality,  and  (3)  neglect  of  some  physical  aspects  of  the  retrieval  problem, 
such  as  the  effect  of  cloud  on  millimeter  wave  brightness  temperatures.  To 
address  these  issues,  an  alternative  retrieval  scheme  has  been  developed  under 
this  effort.  The  approach  outlined  below  is  by  no  means  statistically  opti¬ 
mal.  However,  it  does  attempt  to  address  the  difficulties  cited  above  and,  in 
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particular,  integrates  available  data  sources  in  a  unified,  multispectral 
retrieval  constrained  by  radiative  transfer  principles. 

Table  2-3.  Fraction  of  Unexplained  Variance  for  Statistical  Retrievals: 
Tropical  Summer  Atmospheres 


PARAM. 

SSM/T1 

SSM/T2 

SSM/I 

1  &  2 

1  &  I 

2  &  I 

1,2  &I 

P  S.D. 

P  MEAN 

T 

10 

0.0155 

1.1274 

0.9887 

0.0193 

0.0157 

1.1363 

0.0184 

6.9724 

226 . 3941 

T 

50 

0.1072 

1.0430 

0.9434 

0.1184 

0.1034 

1.0122 

0.1134 

7.1249 

204 . 7910 

T 

100 

0.1492 

0.7357 

0.8642 

0.1528 

0.1624 

0.7120 

0.1511 

5.8208 

195 . 5323 

T 

300 

0.2291 

0.6402 

0.6248 

0.2414 

0.2255 

0.7083 

0.2475 

3.1300 

240.5698 

T 

500 

0.1866 

0.6352 

0.6538 

0.2803 

0.1966 

0.5569 

0.2913 

3.4249 

265.9098 

T 

700 

0.2236 

0.3367 

0.6573 

0.1651 

0.2635 

0.2508 

0.1850 

3.8224 

280.7237 

T 

850 

0.2222 

0.3520 

0.5909 

0.1712 

0.2191 

0.2705 

0.1719 

4.6325 

288.0277 

T 

1000 

0.1784 

0.3640 

0.5450 

0.2107 

0.1787 

0.4027 

0.1991 

5.4435 

296.2430 

U 

300 

0.8514 

0.3737 

0.9687 

0.0929 

0.8678 

0.2939 

0.1102 

0.0042 

0.0116 

U 

500 

0.8684 

0.4197 

0.9950 

0.1286 

0.8185 

0.3360 

0.1546 

0.0848 

0.1515 

u 

700 

0.8632 

0.3982 

0.9431 

0.1168 

0.7694 

0.2876 

0.1389 

0.4503 

0.7835 

u 

1000 

0.5887 

0.4425 

0.4700 

0.1896 

0.4502 

0.2734 

0.2108 

1.2724 

3.7535 

SURF  T 

0.1784 

0.3640 

0.5450 

0.2107 

0.1787 

0.4027 

0.1991 

5.4435 

296.2430 

EMISS 

0.0200 

0.0419 

0.0089 

0.0125 

0.0056 

0.0094 

0.0078 

0.1232 

0.8041 

The  approach  is  based  on  the  simultaneous  retrieval  method  of  Smith  et 
al.  (1985)  and  is  illustrated  schematically  in  Figure  2-3.  The  first  guess 
geophysical  parameters  are  obtained  using  a  multiple  linear  regression  method 
(D-matrix  approach  using  the  technique  described  by  Gaut  et  al.,  (1975)).  Two 
data  sets  of  geophysical  parameters  as  well  as  brightness  temperatures  are 
needed  to  compute  the  first  guess.  One  set  is  used  to  compute  the  D-matrix 
and  the  other  to  compute  the  first  guess.  The  brightness  temperatures  are 
evaluated  using  a  modified  version  of  the  RADTRAN  code.  This  program  computes 
brightness  temperatures  for  microwave  frequencies.  The  versions  used  at  AFGL 
are  called  RAD9  and  RADLIB.  To  do  the  physical  retrieval  (i.e.,  to  correct 
the  first  guess) ,  use  is  made  of  the  radiative  transfer  equation  (RTE) .  The 


RTE  is  varied  with  respect  to  integrated  water  vapor  profile  U,  temperature 
profile  T,  surface  temperature  Tg ,  and  surface  emissivity  <s.  AU  is  expanded 
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Figure  2-3.  Unified  retrieval  algorithm  for  microwave  mission  sensor  data. 


in  terms  of  the  eigenvectors  of  the  standardized  (by  mean)  covariance  matrix 
of  U  of  the  first  data  set.  Similarly,  AT  is  expanded  in  terms  of  the  eigen¬ 
vectors  of  the  standardized  (by  the  mean)  matrix  of  covariance  of  T.  The  un 
knowns  then  are  the  coefficients  of  the  expansions  as  well  as  a  scalar  for  Ts 
and  «s<  After  solving  for  these  coefficients,  it  is  then  possible  to  correct 
the  integrated  water  vapor  and  temperature  profiles,  as  well  as  the  surface 
temperature  and  the  surface  emissivity. 


The  radiative  transfer  equation  for  microwave  frequencies  is: 


Tfci/  —  (^gTs  +  (1  -  £s)  J*q  T(p)dr^)  fj/Cps)  +  /  T(p)drl/ 


where 


r  (p) 
1/ 


-  exP  (-  in  k  (‘'•P'  )dP'/M] 


r  J  ( p)  -  exp  [-  JpS  k(i/,p' )dp'//^l  •  (2- 

Here,  /j  is  the  cosine  of  the  path  zenith  angle,  r ^  and  are  the  upward  and 
downward  transmission  functions,  respectively,  and  pg  is  the  surface  pres¬ 
sure.  Specular  surface  reflection  is  assumed. 

Differentiating  equation  (1)  with  respect  to  the  desired  variables  U,  T, 
Tg  and  «s  and  dropping  the  frequency  indices,  one  obtains: 


<3Tb  3Tb  <?Tb  <9Tb 

577  +  ar  «TS  +  w  6T  +  —  6U 


where : 


3TS  "  ^rPs 


.Pr 

gf  -  rps(1  -  «s>  j  f0  «Tdr  ' 


+  f  6  Tdr 
Ps 


narKmiging 


5U  w  “  (Vl)  !Ts-  Jo  T(P)dr'  +  rsT(0)l  w  6U 


+  V  llu +  rs(£s-1)  fr]  6UdT 


Parts  of  SU  (3T{j/3U)  have  been  obtained  using  integration  by  parts.  Note  that 
Or(ps)/3T)  and  the  variation  with  respect  to  surface  pressure  have  been 
ignored. 

The  quantities  6U  and  6T  are  expanded  in  series  of  the  eigenvectors  of 
the  covariance  matrices  (EOFs)  of  U  and  T,  with  Nu  and  Nt  terms,  respectively. 


N  N 
u  +  t 


SU(p)  -  2  A  4>  (p);  6T(p)  -  2  A.*.(p) 

j-1  J  J  j-N  +1  J  J 


(2.9-2 . 12) 


6Ts  “  %  +N  +1 ’  6es  “  +N  +2 
u  t  u  t 


Upon  substitution  of  (2.9  -  2.12)  into  (2.9)  a  linear  equation  in  the 
coefficients  Aj  is  obtained: 


N  +N  +  2 
u  t 

An  -  2  A.*. 

j-1  J  J 


(2.13) 


where  the  <Dj's  are  functions  of  the  terms  in  equations  (2. 5-2. 8).  The  desired 
difference  terms,  6,  in  the  relaxation  equation: 


n  n-1  ,  . 
p  -  p  +  6 


(2.19) 


are  available  by  solving  (2.13)  for  the  Aj ' s 

The  ridge  stabilized,  least  squares  solution  is  given  by: 


*  ,.T  ,  U.-1.T  n 

A  —  (4>  <t>  +  oH)  4>  A 


(2.15) 


where  a  is  the  ridge  parameter  and  the  diagonal  elements  of  the  H  matrix  are 
the  inverse  of  the  fractional  variance  due  to  each  EOF. 


The  quantities  6U  and  6T  are  expanded  in  terms  of  the  eigenvectors  of  the 
covariance  matrix  of  U  and  T,  respectively.  The  covariance  matrix  is  given  by 


<UUT>  -  [(U.  .  -  U.)/U.  .  -  U.)/U.]T 
1 ,  j  1  l.j  l  l 


(2.16) 


where  T  is  equal  to  the  transpose  and  -  —  E  tL  . ,  i  is  an  index  running 

j  1,J 

over  the  number  of  pressure  levels,  and  the  j  index  is  summed  over  the  number 
of  radiosondes,  N.  The  covariance  matrix  for  temperature  is  the  same,  replac¬ 
ing  U  by  T.  Eigenvectors  and  eigenvalues  are  obtained  by  diagonalizing  the 
above  covariance  matrix.  These  are  stored  in  the  computer  code  as: 


RUUi  , ,  i  -  l.NPEI  (2.17) 

H  l.NU 

RTTj  i .  i  -  l.NPEI 
j  -  l.NT 

The  eigenvalues  are  stored  in  LAHBDU(j),j  -  1,NU  and  LAMBDT(j),j  -  1,NT.  The 
eigenvalues  need  to  be  known  in  the  problem  so  as  to  be  able  to  compute  how 
many  eigenvectors  have  to  be  kept  in  the  computations.  There  are  NPEI  eigen¬ 
vectors  for  both  U  and  T,  but  only  the  first  few  (i.e.,  with  largest  eigen¬ 
values)  are  kept. 

The  perturbations  for  each  quantity  6U,  6T,  6es>  6Tg  (see  eqs .  (2.5)  - 
(2.8))  are  then  expanded  in  the  appropriate  eigenvectors: 

Let  6U£  -  £  RUUt  .*A, ;  j  -  1 ,  NU 

j  l  -  1,  NPEI-64  levels  (2.18) 


6T£  -  £  RTTt  ,*  A.;  j  -  1  +  NU.NU  +  NT 

j  i  -  l.NPEI  (2.19) 

*Tsi  "  ANU+NT+1  (2.20) 


5(s~  ANU+NT+2 


(2.21) 


where  NU  —  number  of  eigenvectors  used  to  compute  6U,  NT  -  number  of  eigen¬ 
vectors  used  to  compute  6T.  Note  that  in  this  problem,  it  has  been  assumed 


that : 


5TS  -  Z  RTT  *  Aj  ^NU+NT+l 
J  J  J 


(2.22) 


Therefore,  there  is  really  one  less  unknown.  Note  that  one  can  also  run 
the  program  if  6TS  is  in  general  different  from  6T  (surface).  Putting 
equation  (2.1)  in  finite  difference  form  and  therefore  using  equations  (2.18- 
2.21),  one  obtains: 


A„  -  6TB*  +  STBS  +  6TB «  +  6TB, 


(2.23) 


l  ~  l.NFREQ. 


where 


hjL  -  6TB if  AT.  ATS,  At  s  -  0 

a"  -  6TB|  if  AU,  AT,  At s  -  0 
T 

Aj  -  6TBis  if  AU,  AT,  A«s  -  0 
A"  -  6TB^*S  if  AU,  AT,  ATS  -  0 


Equation  (2.23)  can  be  solved  for  different  combinations  of  geophysical  para- 

<  s 

meters.  For  example,  if  ts  is  kept  fixed,  then  ATB^  -  0  and  is  not 
computed.  i.  -  l.NFREQ  is  the  number  of  frequencies  used  to  do  the  retrievals 
one  can  choose  between  the  3  instruments  SSM/T-1,  SSM/T-2  and  SSM/I . 

The  desired  quantities  U,  T,  Tg ,  and  fg  are  obtained  by  solving  for  the 
Aj  coefficients  (see  Eqs.  (2.18-2.21)  where  j-1,  NU+NT+2.  From 
equations  (2.5)  through  (2.8),  we  have: 


(a)  5Tb*  -  <V/Ai 


for  i  -  NU+NT+2 
t  -  1 ,  NFREQ 
NLEV  ^  surface  level 


(2.24 


fiitt  "  *■  (NLEV,£)(Ts-  Tf/) 


where : 


T  -  v  T  (t  -  T  J ;  k=l ,  NLEV-1 

1  i  .  k  1  k+l,i  k,i'’ 

k  i-1,  NFREQ 


(b)  «T^S  - 


for  i=NU+NT+ 1 
i- 1 ,  NFREQ 


(2.25) 


B.  -  r (NLEV ,  £)  f 
Ki,£  ! 


<c)  (Tb^  - 


B  -  2  (RTT )T  •  r(NLEV,i)(l-«s) . 

1  ’ 1  k  ’ 1  l 

[ r ' (k+1 , i  )  -  r ' (k, i) ]  +  (r(k,i)  -  r(k+l,i; 


for  i-NU+1,  NU-t NT 
i-1,  NFREQ 

k-1,  NLEV-1 


(2.26) 


(d)  «IB"  -  Wt  y\ 


for  i-1,  NU 

i-1,  NFREQ 
k-1,  NLEV-1 


(2.27) 


Vi  “  RUUNLEV,i  (au}  NLEV-1,  t(TS'1) 

■  t (NLEV, i)  Tx  +  Tg  -  T1£ 

- 


t£(RUU  )  <Tltl-  Tji 
k 


•  (— ) 

^airk, . 


+  t (NLEV, i) 


(V1)  (au  ^  k, i 


where  -  ^R-l 


f,  -  f  S  Tdr* 
1 ,  £  Jo 


-  (r' 

(au  ;h,i  ^  k, i 


-  Tk-i,i)/(uk  -  LW 


y.:  < 


i  -  1 , NU+NT+2 
St  -  1  ,  NFREQ 


“  (rk ,1  '  rk+l,i)/(uk  '  °k+l) 
define  PHIf _ t  -  [Py^ 

Least  squares  solution  of  6TB^  -  Z  PHI^  A^; 


Ai  -  ipBI,.i  <phij,iIt  *  ‘’Hi.i''1ipHI«.iiT  STBf  :  ‘ 

M 

where  H.  .  -  1  /  [A.  /  (  E  A.)l  ;  i  -  l.NU+NT 
l ,  i  '  ‘  1  '  x  1 


1 , NU+NT 
1 , NFREQ 


and  Hnu+nx+x >NU+NX+1  -  1.0;  Hjnj+NT+2 _ NU4NT+2 


(2.28) 

(2.29) 

(2.30) 


2 . 4  Results  for  Retrieval  of  Temperature  and  Integrated  Water  Vapor  Using 
SSM/T-1  and  SSM/T-2 

In  this  section,  we  describe  results  for  the  simultaneous  retrieval  of 
water  vapor  and  temperature  profiles  using  the  SSM/T-1  microwave  temperature 
sounder  and  SSM/T-2  millimeter  wave  moisture  sounder  instruments.  The  channel 
frequencies  and  noise  equivalent  brightness  temperatures  used  in  our  simula¬ 
tions  are  provided  in  Table  2-1.  The  noise  values  were  assumed  Gaussian  with 
a  zero  mean,  i.e.,  no  instrumental  bias  was  considered.  The  adopted  noise 
values  in  Table  2-1  represented  design  goals  and  probably  underestimate  actual 
sensor  noise  including  scene  noise.  Therefore,  the  retrieval  results 
presented  should  be  viewed  as  best  case  retrievals. 

The  atmospheric  profiles  of  temperature  and  water  vapor  used  in  the 
brightness  temperature  simulations  are  from  a  radiosonde  data  set  collected  by 
Phillips  et  al.,  1988.  Subsets  of  the  original  data  set  were  created 
consisting  of  200  radiosondes  each  corresponding  to  tropical  winter  profiles 
over  the  ocean  (TWO)  and  midlatitude  winter  profiles  over  the  ocean  (MWO) , 
respectively.  One  half  of  each  set  (denoted  set  1)  was  used  to  develop  first 
guess  and  EOF  statistics,  while  the  other  half  (i.e.,  set  2)  was  used  to 
perform  retrievals. 

Below  we  describe  the  results  for  the  simultaneous  retrieval  approach. 

The  first  guess  was  derived  either  from  statistics  or  using  the  ensemble 
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mean.  The  corrections  to  the  first  guess  were  obtained  using  the  eigenvector 
method  with  1  iteration  only  (in  general). 

2.4.1  Statistical  Estimates  as  First  Guess 

Multiple  linear  regression  (i.e.,  the  D-matrix  approach)  was  used  to 
obtain  first  guess  estimates  of  the  desired  geophysical  parameters.  To  obtain 
these  estimates,  all  the  channels  of  the  SSM/T-1  and  SSM/T-2  are  used.  To 
correct  the  first  guesses,  two  options  were  considered.  The  first  used  a 
selective  number  of  channels  for  SSM/T-1  and  a  scalar  to  compute  the  correc¬ 
tions  rather  than  eigenvectors.  The  second  method  used  all  the  channels  of 
the  SSM/T-1  and  the  eigenvectors  to  compute  the  corrections. 

2 . 4 . 1 . 1  Selective  Channels  for  the  SSM/T-1 

Retrievals  for  temperature  and  integrated  water  vapor  have  been  done 
using  all  channels  for  SSM/T-i.  me  corrections  obtained  to  the  statistical 
first  guesses  give  better  results  than  the  first  guess  for  surface  variables 
only  (i.e.,  the  layer  1-1000  mb  for  water  vapor  and  the  air  temperature  at 
1000  mb  as  well  as  the  ground  temperature). 

Therefore,  the  eigenvector  method  for  the  temperature  profile  retrieval 
has  been  changed  to  a  scalar  by  only  correcting  the  surface  temperature.  For 
the  Integrated  water  vapor  corrections,  the  eigenvectors  of  the  covariance 
matrix  for  water  vapor  are  still  used.  For  these  computations  and  the  ones 
following,  two  new  rejection  criteria  for  integrated  water  vapor  have  been 
added.  These  are  number  2  and  3  in  the  list  following. 

Criteria  for  rejection  of  updating  the  first  guess  of  an  integrated  water 
vapor  profile: 

Reject  updated  profile  if: 

1.  The  corrected  water  vapor  profile  has  negative  values  of  integrated  water 

vapor  (at  one  of  the  64  levels). 

2.  The  integrated  water  vapor  for  the  layers  1  - 300mb , 1  - bOO , 1 - 700  and 

l-1000mb  is  nonmonotonic. 


3. 
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The  water  vapor  for  the  layer  l-700mb  is  larger  than  0.4  times  the  water 
vapor  for  the  layer  l-1000mb.  Table  2-4  contains  the  ratios  of  the  mean 
of  integrated  water  vapor  for  all  the  layers  --  U(  1  -  700mb)/U(  1  -  lOOOmb)  is 
around  20%  for  both  data  sets. 

If  a  particular  integrated  water  vapor  profile  satisfies  one  of  these 
criteria,  then  the  newly  computed  values  of  integrated  water  vapor  are  set 
back  to  the  first  guess. 


Table  2-4.  Ratios  of  Mean  Integrated  Water  Vapor 


TW01 

TU02 

MW01 

MWO? 

U(l)/U(4) 

0.3% 

0.3% 

0.54% 

0.54% 

U(2)/U(4) 

4.0% 

4.0% 

5.0% 

5.3% 

U(3)/U(4) 

20.1% 

20.9% 

22.3% 

23.5% 

U(l)-U(l- 300MB) ,  U(2)-U(l- 500MB ) ,  U( 3 )-U( 1 - 700MB) ,  U(4 )=U( 1  -  1000MB) 


V 

i 
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The  results  for  the  selective  channel  method  are  shown  in  Table  2-5.  The 
layer  l-500mb  shows  a  small  improvement  and  the  total  integrated  water  vapor 
shows  almost  no  improvement.  The  layer  500-700mb  shows  no  improvement  at  all 
The  surface  temperature  shows  a  slight  improvement  for  the  midlatitude  winter 
set.  (Note  that  all  temperature  levels  except  for  the  surface  have  been  kept 
fixed. ) 


2. 4. 1.2  All  Channels  for  SSM/T-1 
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The  first  method  used  to  do  retrievals  of  integrated  water  vapor 
and  temperature  used  eigenvectors  for  both  temperature  and  water  vapor.  How¬ 
ever,  the  number  of  eigenvectors  used  for  temperature  was  high  (i.e.,  6  to  7). 
The  results  obtained  were  disappointing  and  therefore  the  method  of  using  a 
scalar  only  for  the  temperature  correction  was  tried.  It  turns  out  that  using 
eigenvectors  for  both  temperature  and  integrated  water  vapor  gives  better 
results,  provided  a  lower  number  of  eigenvectors  is  used  for  the  temperature. 
The  results  improve  only  for  the  surface  temperature  (ground  temperature)  and 
the  total  integrated  water  vapor.  Tables  2-6  and  2-7  show  results  for  dif- 
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ferent  values  of  sigma,  H  and  number  of  eigenvectors  for  tropical  and  mid- 
latitude  cases,  respectively.  Both  noise  and  no  noise  cases  are  shown. 


Table  2-5.  Simultaneous  Retrieval  Results  Using  Statistical  First  Guess  and 
Selective  Channels. 


Data  set: 

TV02 

MW02 

H2C:  #  of  eig's 

4 

3 

T:  #  of  eig's 

0 

0 

SSM/T-1  channels 

50.5,53. 

2  GHz 

50.5,53 

.2,54. 

35  GHz 

sigma 

le-05 

le-05 

U  (gr/sq  cm) 

RMS0E(*) 

RMSTE 

RMSCLI 

RMSOE 

RMSFE 

RMSCLI 

l-300mb 

.00129 

.00129 

.00423 

.00085 

.00085 

.00261 

1  - 500mb 

.0304 

.0267 

.0848 

.0154 

.0132 

.0400 

1-7  OOmb 

.  154 

.240 

.450 

.0675 

.0824 

.  186 

l-1000mb 

.554 

.543 

1.27 

.248 

.230 

.  555 

T  (K) 

10  mb 

.97 

6.97 

1.62 

6.77 

50  mb 

2.45 

7.13 

1.52 

6.16 

100  mb 

2.28 

5.82 

1.78 

5.84 

300  mb 

1.63 

3.31 

1.56 

4.25 

500  mb 

1.81 

3.43 

2.78 

6.13 

700  mb 

1.55 

3.82 

2.24 

6.26 

850  mb 

1.92 

4.63 

3.10 

00 

lT> 

lOOOmb 

2.50 

2.82 

5.44 

4.11 

3.31 

6.86 

(*)  Tolerance  -  tolerance  criterion  to  stop  computations. 

RMSOE  -  rms  error  between  first  guess  and  exact  solution. 
RMSFE  -  rms  error  between  final  and  exact  solution. 

RMSCLI  -  rms  error  of  climatology  alone. 
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Table  2-6.  Simultaneous  Retrieval  Results  for  Tropical  Data  Set:  (a)  Noise, 
(b)  No  Noise. 


(a)  Instrument 

Noise  Added 

H  #  of  eig's: 

4 

4 

4 

T  #  of  eig’s: 

4 

3 

4 

sigma 

le  -  5 

.  1 

.  1 

H 

A 

I 

I 

U(gr/sq  cm) 

RMS0E 

RMSFE 

RMSFE 

RMSFE 

1 - 300mb 

.00129 

.00129 

.00129 

.00129 

1 - 500mb 

.0304 

.0302 

.0303 

.0310 

1- 700mb 

.154 

.216 

.308 

.259 

1 - lOOOmb 

.554 

.  531 

.439 

.493 

T  (K) 

10  mb 

.97 

1.34 

1.52 

1.46 

50  mb 

2.45 

2.45 

2.55 

2.51 

100  mb 

2.28 

3.21 

2.30 

3.76 

300  mb 

1.63 

2.29 

1.92 

1.98 

500  mb 

1.81 

1.95 

2.13 

1.96 

700  mb 

1.55 

1.89 

1.99 

2.11 

850  mb 

1.92 

2.15 

2.33 

2.23 

lOOOmb 

2.50 

2.32 

2.36 

2.21 

'b)  No  Instrument  Noise 


H2O  #  of 

eig's: 

4 

4 

4 

4 

4 

T  #  of 

eig's: 

4 

4 

3 

3 

3 

sigma 

le-5 

le  -  3 

le  -05 

le  -  03 

le  - 10 

H 

A 

A 

A 

A 

A 

U(gr/sq  cm)  RMS0E 

RMSFE 

RMSFE 

RMSFE 

RMSFE 

RMSFE 

1 - 300mb 

.00100 

.00100 

.00100 

.00100 

.00100 

.00100 

1 - 500mb 

.0256 

.0264 

.0257 

.0253 

.0251 

.0253 

1- 700mb 

.116 

.159 

,214 

.223 

.195 

.227 

1- lOOOmb 

.425 

.405 

.419 

.374 

.409 

.374 

T  (K) 

10  mb 

.526 

1.04 

1.02 

1.11 

1.08 

1.11 

50  mb 

1.81 

1.73 

1.97 

1.92 

2.02 

1.91 

100  mb 

1.57 

2.40 

3.51 

2.28 

1.81 

2  .  31 

300  mb 

1.06 

2.10 

1.56 

1.52 

1  .  29 

1.53 

500  mb 

1.31 

1.56 

1.57 

1.82 

1.71 

1  .  83 

700  mb 

1.33 

1.68 

2  .  15 

2.00 

1.85 

2.01 

850  mb 

1.64 

1.83 

2.12 

2.08 

2.10 

2.08 

lOOOmb 

2.31 

2.00 

1.97 

2.18 

2.25 

2.  19 

% 
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Table  2-7.  Simultaneous  Retrieval  Results  for  Midlatitude  Data  Set: 


(a)  Noise 

(b)  No 

'Joise 

(a)  Instrument  Noise 

Added 

sigma : 

le-10 

le-05 

le-05 

le-02 

1  e  -  0 1 

H: 

A 

A 

A 

A 

A 

H20  #  of 

eig 

's: 

3 

3 

3 

3 

3 

T  #  of 

eig 

's: 

4 

3 

4 

4 

4 

U  (gr/sq 

cm) 

RMSOE 

RMSFE 

RMSFE 

RMSFE 

RMSFE 

RMSFE 

1  - 300mb : 

.00085 

.00085 

.00085 

.00085 

.00085 

.00085 

1  - 300mb : 

.01538 

.01601 

.01503 

.01539 

.01416 

.01417 

1 - 700mb : 

.06752 

.11226 

.13287 

.11208 

.12035 

.  12056 

1- lOOOmB 

.24814 

.27720 

.  28392 

.27360 

.25913 

.25983 

T  (K) 

10  mb: 

1.62 

1.67 

. 

1.71 

1 .84 

1.85 

50  mb: 

1.52 

1.54 

- 

1.56 

1.64 

1.64 

100  mb: 

1.78 

2.  18 

- 

2.40 

2.06 

2.06 

300  mb: 

1.96 

3.40 

- 

3.68 

3.02 

3.01 

500  mb : 

2.78 

3.70 

- 

3.63 

2.99 

3.00 

700  mb: 

2.24 

3.18 

- 

3.39 

3.16 

3.16 

850  mb: 

3.10 

3.69 

- 

3.80 

3.83 

3.83 

lOOOmb : 

4.11 

3.63 

3.67 

3.44 

3.35 

3.36 

(b)  No  Instrument  Noise 


sigma : 

le-10 

le-05 

le-05 

le-02 

le-01 

H20  #  of 

eig's: 

3 

3 

3 

3 

3 

T  *  of 

eig's: 

4 

3 

4 

4 

4 

H: 

A 

A 

A 

A 

A 

RMSOE 

RMSFE 

RMSFE 

RMSFE 

RMSFE 

RMSFE 

U(gr/sq 

cm) 

1  - 300mb : 

.00071 

.00071 

.00071 

.00071 

.00071 

.00071 

1  - 500mb : 

.01413 

.01425 

.01384 

.01438 

.01306 

.01301 

1  - 700mb : 

.06464 

.08961 

.10642 

.09317 

.09355 

.09571 

1- lOOOmb 

:  .23345 

.  19612 

.20735 

.  19482 

.19483 

.  19571 

T  (K) 

10  mb 

1.42 

1.41 

1  .  53 

1.44 

1.53 

1 .  53 

50  mb 

1.13 

1 . 12 

1.23 

1 .16 

1.28 

1.28 

100  mb 

1.46 

1.64 

1  .  70 

1  .  56 

1.66 

1.67 

300  mb 

1.55 

2.41 

1.59 

2.41 

2.02 

2.02 

500  mb 

2.35 

2.98 

2.86 

2.99 

2.61 

2.61 

700  mb 

1.57 

2.56 

2.53 

2.75 

2.33 

2.33 

850  mb 

2.85 

3.42 

3.41 

3.50 

3.49 

3.49 

lOOOmb 

2.80 

2.62 

2.52 

2.41 

2.29 

2.30 

The  results  depend  a  lot  on  the  number  of  eigenvectors  used  as  well  as  on 
the  value  of  sigma. 

A  method  using  a  different  criterion  for  convoigence  has  also  been  used. 
The  criterion  used  so  far  is  based  on  the  total  rms  of  the  brightness  tempera¬ 
ture  residuals.  By  the  total,  we  imply  the  residuals  for  all  the  frequencies, 
SSM/T-2  as  well  as  SSM/T-1. 

The  new  criterion  consist  of  3  instead  of  1  criterion.  The  first  one 
consists  of  the  rms  for  the  residuals  of  brightness  temperatures  for  the 
SSM/T-2.  The  second  one  consists  of  the  rms  for  the  residuals  of  brightness 
temperatures  for  the  SSM/T-1.  The  third  one  consists  of  the  brightness 
temperature  for  the  50.5  GHz  channel. 

Conditions  for  possible  update  of  first  guess: 

TBRMS(O)  -  RMS  of  residuals  of  first  guess  -  synthetic 
TBRMS(O)  SSM/T-2  >  IK  then  possibly  update  water  vapor 
TBRMS(O)  SSM/T-1  <  0.5  K  and  TBRMS  50.5  GHz  >  0.6  K  then  possibly 
update  surface  temperature 

TBRMS(O)  SSM/T-1  >  0.5  K  then  possibly  update  surf.  temp,  and  temp. 

Above  possibility  rejected  if: 

TBRMS(l)  =  rms  of  residuals  of  corrected  -  synthetic 
If  TBRMS  SSM/T-2  diverges  then  do  not  update  water  vapor 
If  TBRMS  50.5  GHz  diverges  then  do  not  update  sur.  temp. 

If  TBRMS  SSM/T-1  diverges  then  do  not  update  temperature  profile 
(does  not  include  surface  temperature) 

Results  are  shown  in  Table  2-8a,b  for  noise  and  no  noise  cases  and  2-9a,b  for 
tropical  midlatitude  noise  and  no  noise  cases. 
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Table  2-8.  Results  for  New  Convergence 

Criteria , 

Tropica 

Atmosph 

(a)  Noise,  (b) 

No  Noise 

(a)  With  Instrument  Noise 

H2O:  4  eigenvectors 

T:  4  eigenvectors 

First  guess  -  statistical 
T0L(1)  -  TOL(SSM/T- 1)  -  0. 

5 

T0L(2)  -  TOL( SSM/T- 2 )  -  1. 

0 

Tolerance  on  surface  channel  *  50.5  GHz 

channel  — 

0.6  K 

X 

I 

sigma  - 

sigma  - 

H 

le  -05 

0.1 

U  (gr/sq  cm) 

RMS0E 

RMSFE 

RMSFE 

1  - 300mb : 

.00129 

.00129 

.00129 

1 - 500mb : 

.03040 

.03119 

.03062 

1 - 700mb : 

.15387 

.24773 

.31013 

1  - lOOOmb : 

.55407 

.52818 

.51102 

T  (K) 

10  mb : 

.96716 

1.22472 

1.33096 

50  mb : 

2.45061 

2.40394 

2.55526 

100  mb: 

2.27512 

2.67906 

2.85289 

300  mb: 

1.62588 

2.07338 

1.85120 

500  mb- 

1.81382 

1.89326 

1.92000 

700  mb: 

1.55331 

1.72099 

1.88397 

850  mb: 

1.91596 

2.03671 

2.09029 

lOOOmb : 

2.49872 

1.97584 

2.17187 

V’  \  \  % 


Table  2-8.  Results  for  New  Convergence  Criteria,  Tropical  Atmospheres: 
(a)  Noise,  (b)  No  Noise  (continued) 


(b)  Without  Instrument  Noise 

First  guess  -  statistical 
T0L(1)  -  TOL( SSM/T- 1 )=0 . 5 
T0L(2)  -  T0L(SSM/T-2)“1 .0 
NITER  -  1.0 

Tolerance  on  surface  channel  convergence  -  0.6,  surface  channel  -  50.5  GHz 


H: 

X 

X 

X 

X 

H20  #  of  eig’s: 

4 

4 

4 

4 

T  #  of  eig' s : 

3 

3 

3 

4 

Sigma: 

.1 

le-5 

le-10 

le  -  05 

RMS0E 

RMSFE 

RMSFE 

RMSFE 

RMSFE 

U(gr/sq  cm) 

1 - 300mb : 

.00100 

.00100 

.00100 

.00100 

.00100 

1  -  500mb : 

.0256 

.0250 

.0259 

.0259 

.0262 

l-700mb: 

.116 

.208 

.222 

.226 

.290 

l-1000mb: 

.425 

.399 

.387 

389 

.380 

T  (K) 

10  mb : 

.526 

1.04 

1.13 

1.14 

.89 

50  mb: 

1.81 

2.00 

1.92 

1.92 

1.75 

100  mb: 

1.57 

1.66 

2.33 

2.35 

1.91 

300  mb: 

1.06 

1.13 

1.49 

1.50 

1.73 

500  mb: 

1.31 

1.52 

1.72 

1.73 

1.48 

700  mb: 

1.33 

1.54 

1.81 

1.82 

1.63 

850  mb: 

1.64 

1.88 

1.80 

1.80 

1.78 

lOOOmb : 

2.31 

2.31 

1.82 

1.82 

1.78 

Table  2-9.  Results  for  New  Convergence  Criteria,  Midlatitude  Atmospheres: 
(a)  Noise,  (b)  No  Noise 


1 

Es 


With  Instrument  Noise 


Midlatitude  winter  ocean  set  2 
H?0  #  of  eigenvectors:  3 
T  #  of  eigenvectors:  4 
First  guess  -  statistical 
T0L(1)  -  TOL(SSM/T- 1) 

T0L(2)  -  TOL(SSM/T- 2 ) 

NITER  -  1 

Tolerance  on  surface  channel  - 


50.5  GHz  channel  « 


T0L(1) : 
T0L(2) : 

H: 

sigma: 

U(gr/sq  cm) 
1 - 300mb 
1 - 500mb 
1-7  OOmb 
l-1000mb 


1.0 

0.5 

0.0 

1.0 

A 

A 

le-05 

le-05 

RMSOE 

RMSFE 

RMSFE 

00085 

.00085 

.00085 

01538 

.01565 

.01656 

06752 

.11077 

.12348 

24814 

.26563 

.26780 

■Vi'Li't'ilt'VtVl'. 


Table  2-9.  Results  for  New  Convergence  Criteria,  Midlatitude  Atmospheres 
(a)  Noise,  (b)  No  Noise  (continued) 


Without  Instrument  Noise 


T  (K) 


:  1.42 

1.41 

1  .46 

3 

:  1.13 

1.13 

1  .  12 

:  1.46 

1.50 

1.58 

5 

:  1.55 

1.73 

1.97 

y 

:  2.35 

2.48 

2.74 

< 

V 

V 

:  1.57 

1.75 

2.04 

:  2.85 

2.95 

3.16 

:  2.80 

1.84 

2.14 

n 

.4.2  Mean  Values  of  First  Set 

as  First  Guess 

If  full  statistical  data  to  develop  regression  coefficients  for  a 
statistical  retrieval  is  not  available,  a  climatological  mean  profile  might  be 
an  appropriate  first  guess.  In  this  section  we  investigate  the  use  of  a  mean 
first  guess.  The  mean  of  data  set  1  for  both  temperature  and  integrated  water 
vapor  profiles  is  used.  The  means  are  given  at  4  levels  for  water  vapor  and  8 
levels  foi  temperature.  These  values  are  then  interpolated  to  the  64  levels 
used  by  "RADTRAN." 

Since  the  first  guess  is  set  equal  to  the  mean  of  set  1 ,  the  RMS  error  of 
set  2  is  approximately  equal  to  the  standard  deviation  of  set  2.  The  eigen¬ 
vector  method  was  used  for  both  temperature  and  integrated  water  vapor.  This 
permits  corrections  at  all  levels  for  both  temperature  and  integrated  water 
vapor.  Results  are  shown  in  Table  2-10  for  one  iteration  and  adding  instru¬ 
ment  noise.  In  this  calculation  only  the  total  RMS  convergence  criteria  was 


used.  Significant  improvement  is  obtained  for  the  total  integrated  water 
vapor  for  both  data  sets.  For  the  tropical  winter  over  ocean  data  set,  the 
percentage  of  improvement  over  the  first  guess  is  22%  and  for  midlatitude 
winter  ocean  it  is  14% .  The  temperatures  for  both  data  sets  also  show 
improvement.  Levels  with  the  largest  improvements  are  10  and  50  mb  and  the 
surface  for  the  midlatitude  winter  over  ocean  data  set.  These  results  show 
that  the  physical  retrieval  is  working  since  the  first  guess  is  quite  bad. 


Table  2-10.  Retrieval  Results  Using  Mean  as  a  First  Guess  (Single  Convergence 
Criteria  Only) 


Data  set: 

TWO  2 

MW02 

H20  #  of  eig's: 

4 

3 

T  #  of  eig's: 

3 

3 

sigma: 

le-05 

le-05 

H: 

A 

X 

#  of  iter's: 

1 

1 

Tolerance : 

1.0 

1.0 

U(gr/sq  cm) 

RMS0E 

RMSFE 

RMSCLI 

RMS0E 

RMSFE 

1 - 300mb : 

.00423 

.00423 

.00423 

.00263 

.00263 

1 - 500mb : 

.08483 

.07983 

.08479 

.04064 

.03846 

l-700mb: 

.45093 

.52203 

.45030 

.18775 

.17852 

l-1000mb: 

1.2729 

.99557 

1.2724 

.55725 

.47698 

T  (K) 

10  mb : 

6.98 

5.57 

6.93 

5.57 

50  mb: 

7.14 

5.0 

6.31 

4.81 

100  mb: 

5.87 

6.2 

5.97 

4.55 

300  mb: 

3.33 

3.15 

4.26 

4.00 

500  mb: 

3.52 

3.39 

6.17 

5.42 

700  mb: 

3.88 

3.78 

6.34 

5.60 

850  mb: 

4.64 

4.8 

5.86 

5.37 

1000  mb: 

5.46 

5.43 

6.87 

5.73 

A  selective  criterion  of  convergence  has  also  been  used  with  the  first  guess 
equal  to  the  mean.  The  results  are  shown  in  Table  2-11.  The  criterion  is  the 
same  as  the  one  described  above. 


Table  2-11.  Retrieval  Results  Using  Selective  Criteria  for  Convergence 
(a)  Noise,  (b)  No  Noise 


Midlatitude  winter  ocean  set  2 
First  guess  -  mean  of  data  set  1. 

Instruments:  SSM/T-1  and  SSM/T-2 

(Use  three  criteria  of  rejection  for  water  vapor  profiles) 
T0L(1)  -  TOL(SSM/T - 1)  -  0.5 
T0L( 2 )  -  T0L(SSM/T-2)  -  1.0 


(a)  Noise  Added 

H: 

1 

A 

sigma: 

le-05 

le-05 

H2O  #  of  eig's: 

3 

3 

T  #  of  eig's: 

4 

3 

U(gr/  sq  cm) 

RMSOE 

RMSFE 

RMSFE 

1- 300mb 

.00263 

.00263 

.00263 

1  - 500mb 

.04042 

.03954 

.03876 

1  - 700rab 

.18775 

.23193 

.17993 

1  - lOOOmb 

.55725 

. 50205 

.47380 

Table  2-11.  Retrieval  Results  Using  Selective  Criteria  for  Convergence: 
(a)  Noise,  (b)  No  Noise  (continued) 

(b)  No  Noise 


H:  X  X 


sigma : 

le-5 

le-5 

H20  #  of  eig's: 

3 

3 

T  #  of  eig's: 

4 

3 

RMS0E 

RMSFE 

RMSFE 

U(  gr/  sq  cm) 

1 - 300mb 

.00263 

.00263 

.00263 

1 - 500mb 

.04042 

.03996 

.03786 

1 - 700mb 

.18775 

.23977 

.17884 

l-1000mb 

.55725 

.48121 

.46649 

T  (K) 

10  mb 

6.93 

4.55 

5.23 

50  mb 

6.31 

4.05 

4.55 

100  mb 

5.97 

4.63 

4.21 

300  mb 

4.26 

6.04 

3.91 

500  mb 

6.17 

7.85 

5.38 

700  mb 

6.34 

7.39 

5.54 

850  mb 

5.86 

6.16 

5.28 

lOOOmb 

6.87 

5.28 

5.67 

2.4.3  Results  Summary 


As  indicated  in  the  previous  sections,  a  physical  retrieval  step  such  as 
that  described  in  Section  2.3  can  improve  the  accuracy  of  the  statistically 
based  first  guess  retrieval  results.  The  tabulated  results  show  that  the 
degree  of  improvement  is  a  function  of  the  parameters  of  the  physical 
retrieval  step  such  as:  (a)  the  number  of  eigenvectors  used  to  express  the 
temperature  and  water  vapor  profiles,  (b)  the  ridge  parameter  (i.e.  a)  used  in 
the  retrieval,  (c)  the  form  of  the  smoothing  matrix  (i.e.  either  the  identity 
matrix  or  that  based  on  eigenvalue  scaling),  (d)  and  the  constraints  used  to 
terminate  the  iterations  based  on  the  residuals  (A1?  ) .  In  general  it  is  noted 
that  a  few  eigenvectors  (3  or  4  for  temperature,  and  2  or  3  for  water  vapor), 
a  ridge  parameter  of  about  10'^,  the  eigenvalue  smoothing,  and  separate 
(rather  than  total  RMS)  constraints  seem  to  perform  best. 


After  concluding  the  series  of  calculations  described  above  an  error  in 
the  computer  code  was  corrected  which  improved  the  results  reported  above. 

Results  ot  the  physical  retrieval  ior  the  tropical  temperature  and 
integrated  water  vapor  profiles  are  shown  in  Figures  2-4a,b,  respectively. 

The  statistical  first  guesses  described  in  the  previous  section  were  employed 
and  the  necessary  EOFs  of  temperature  and  water  vapor  are  evaluated  from  the 
radiosonde  data  set  used  in  the  generation  of  the  first  guess  retrieval 
statistics.  Physical  adjustment  of  the  temperature  profile  results  in  a 
slight  increase  of  RMS  error  in  the  vicinity  of  the  tropopause  and  a  small 
improvement  in  accuracy  near  the  surface.  The  tropopause  problem  is  attri¬ 
butable  both  to  EOF  truncation  representation  errors  and  the  generally  weak 
contribution  functions  in  this  region.  There  is  a  much  more  noticeable 
improvement  on  the  accuracy  of  water  vapor  results.  Retrieval  of  low  level 
moisture  is  considerably  improved  by  the  physical  adjustment  applied  to  the 
statistical  first  guess  results.  This  is  understandable  since  the  physical 
retrieval  process  directly  treats  the  inherently  nonlinear  dependence  of 
channel  brightness  temperatures  on  water  vapor  which  is  not  well  represented 
by  the  first  guess  linear  regression. 
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Figure  2-4.  Physical  retrieval  results  for  the  tropical  cases. 
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To  investigate  the  cloud  signal  on  the  various  DMSP  SSM/T-2  channels, 
both  sensitivity  and  retrieval  calculations  were  done.  Channel  weighting 
function  and  brightness  temperature  sensitivities  to  a  cloud  model  were 
evaluated  by  inserting  a  cloud  within  a  tropical  model  atmosphere  with  mono- 
tonically  increasing  values  of  the  integrated  liquid  water  content.  Five 
different  cloud  models  were  investigated.  These  included  altocumulus, 
stratocumulus ,  stratus,  nimbostratus ,  and  cumulus.  Additionally,  two  values 
of  surface  emissivity  were  used,  07  and  1.0,  to  simulate  the  ocean  surface 


and  land,  respectively.  Here  the  potential  effects  of  beam  filling  clouds  on 
millimeter  wave  moisture  retrieval  accuracy  are  investigated  and  results 
obtained  in  ~loudy  conditions  are  compared  with  thrir  clear  "mintern.irt'’ . 

2.5.1  Cloudy  Moisture  Sounder  Channel  Weighting  Functions 

The  effects  of  cloud  on  the  optical  properties  of  the  atmosphere  for  the 
millimeter  wave  sounder  channel  frequencies  were  assessed  by  evaluating  sensor 
weighting  functions.  These  calculations  were  performed  both  with  and  without 
cloud  so  that  the  cloud  effect  could  be  quantified.  Profiles  of  atmospheric 
absorption  coef f icients ,  k(i/,p),  for  each  frequency  of  interest  were  cal¬ 
culated  using  the  RADTRAN  simulation  algorithm  (Falcone  et  al. ,  1982).  These 
provided  the  transmission  profiles  (Eqns.  2.2  and  2.3)  necessary  to  evaluate 
the  weighting  function  which  was  defined  as  the  change  of  transmittance  with 
height  in  the  atmosphere  with  unit  transmittance  taken  at  the  top  of  the 
atmosphere  (in  the  case  of  this  calculation,  at  40  km).  The  channel  set 
selected  consisted  of  the  five  frequencies  corresponding  to  those  of  the  DMSP 
SSM/T-2  moisture  sounder.  These  frequencies  and  the  design  noise  equivalent 
brightness  temperatures  (NEAT's)  for  the  SSM/T-2  are  given  in  Table  2-1.  To 
simulate  sensor  and  scene  noise  for  retrieval  studies,  a  random  signal  of  IK 
(Gaussian  distribution  with  zero  mean)  was  added  to  each  simulated  brightness 
temperature.  To  simulate  the  effect  of  cloud  within  the  f iel d-of -view  of  the 
radiometer  a  suitable  set  of  cloud  models  was  incorporated  within  the  evalua¬ 
tion  of  atmospheric  attenuation  profiles.  The  cloud  models  chosen  were  taken 
from  the  AFGL  FASCODE  model  (Falcone  et  al.,  1979).  Cloud  attenuation, 
a(i/, T,z),  was  calculated  using  a  simple  relationship  (Staelin,  1^66)  between 
cloud  liquid  water  content  LWC  and  a  frequency  and  temperature  dependent  pro¬ 
portionality  constant  f(i/,T)  dependent  on  cloud  index  of  refraction. 


Index  of  refraction  data  for  water  are  from  (Ray,  1972)  while  those  for  ice 
are  from  (Warren,  1983).  Cloud  liquid  water  content  and  vertical  extent  for 
the  five  cloud  models  chosen  from  (Faicone  et  a  1.,  1979)  are  given  in 
Table  2-12. 


Table  2-12.  Cloud  Types,  Liquid  Water  Contents,  and  Vertical  Extents  Used 
in  Cloudy  Brightness  Temperature  Simulations 


The  effect  of  an  altostratus  cloud  (model  3)  on  the  profiles  of  atmo¬ 
spheric  optical  properties  at  these  millimeter  wave  frequencies  is  illustrated 
in  Figures  2-6  through  2-11.  Plotted  are  the  weighting  functions  (-  for 
six  model  atmospheres  for  the  channel  frequencies  given  in  Table  2-1  both 
without  (a)  and  with  (b)  an  altostratus  cloud  in  the  field-of -view.  The  six 
model  atmospheres  are:  (1)  U.  S.  standard,  (2)  tropical,  (3)  mid-latitude 
summer,  (4)  mid-latitude  winter,  (5)  sub-arctic  summer,  and  (6)  sub-arctic 
winter.  The  weighting  function  indicates  the  relative  contribution  of  each 
atmospheric  level  to  the  brightness  temperature  at  a  given  frequency.  It  can 
be  seen  from  these  figures  that  the  cloud  seriously  restricts  the  sensor's 
ability  to  "sound"  the  atmosphere  by  constraining  brightness  temperature 
contributions  to  levels  in  the  vicinity  of  the  cloud  (2.5  -  3.0  km  for  the 
altostratus  cloud) . 

Figures  2-6a  through  2- 11a  provide  a  useful  comparison  of  the  effect  of 
model  atmosphere  on  the  SSM/T-2  weighting  functions.  The  change  in  height 
distribution  of  each  channel  with  temperature  and  water  vapor  profile  can  be 
quite  drastic.  As  an  example,  compare  the  location  of  the  peak  of  the 
183 . 3 1± 3  GHz  channel  weighting  function  (labeled  as  180)  for  tropical 
(Fig.  2-7a)  and  midlatitude  winter  (Fig.  2-9a)  atmospheres .  The  peak  shifts 
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from  7.5  km  for  the  tropical  atmosphere  to  4 . 5  for  the  midlatitude  primarily 
due  to  the  change  in  water  vapor  amount.  Similar  behavior  can  be  noted  for 
the  other  channels. 

Significantly,  the  effect  of  the  cloud  layer  on  the  weighting  functions 
is  apparent  from  Figures  2-6b  through  2-llb.  Since  the  product  of  the  weight¬ 
ing  function  and  the  temperature  profile  determine  the  contribution  of  each 
atmospheric  level  to  the  brightness  temperature  at  a  given  frequency,  the 
highly  cloud  peaked  weighting  fun onion  suggests  an  important  impact  on  the 
SSM/T-2  channels  due  to  clouds.  These  figures  provide  an  indication  of  which 
channels  will  be  affected  by  clouds  at  various  altitudes.  Due  to  the  shift  of 
weighting  functions  with  the  model  atmosphere  noted  above,  the  number  of  af¬ 
fected  channels  will  change.  For  example,  since  channels  peak  higher  in  the 
tropics  they  will  be  least  impacted  by  low  cloud  with  the  converse  being  time 
at  higl  a*-  latitudes  where  integrated  water  vapor  amounts  are  less. 

2.5.2  Brightness  Temperature  Simulations 

Here  we  investigate  the  effect  of  cloud  liquid  water  on  brightness  tem¬ 
peratures.  The  brightness  temperatures  were  evaluated  from  Eqn.  (2.1)  using  a 
radiosonde  ensemble  corresponding  to  either  a  midlatitude  or  tropical  sounding 
over  both  land  and  ocean  surfaces.  Surface  temperatures,  T  ,  for  each  simula¬ 
tion  were  assumed  equal  to  the  radiosonde  temperature  at  1000  mb  which  was 
taken  as  the  surface  pressure,  ps .  Two  values  of  surface  emissivity,  «s,  were 
used:  one  for  land  backgrounds  (1.0)  and  one  for  oceanic  backgrounds  (0.7). 

Clear  air  attenuation  due  to  absorbing  gases  was  evaluated  using  RADTRAN  as 
described  in  the  previous  section.  Clouds  were  incorporated  into  the  bright¬ 
ness  temperature  simulation  process  by  randomly  selecting  from  among  one  of 
the  cloud  types  in  Table  2-12.  Cloud  attenuation  calculated  from  Eqn.  2.16 
was  then  added  to  the  clear  sky  absorption  for  each  atmospheric  layer  within 
the  vertical  domain  specified  for  the  selected  cloud  type  model.  It  was  as¬ 
sumed  that  all  cloud  t}  js  were  equally  probable,  completely  filled  the  field- 
of-view,  and  that  only  a  single  cloud  layer  was  present.  The  second  assump¬ 
tion  tends  to  overestimate  cloud  impact,  while  the  last  results  in  an  under¬ 
estimate.  Since  the  profiles  of  temperature  and  water  vapor  were  given  at 
constant  pressure  levels  rather  than  at  constant  heights  as  are  the  cloud 
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Figure  2  -  6 .  Z 3M/T - 2  weighting  functions  for  U.S.  standard  profile: 
(a)  clear,  (b)  cloudy. 


Figure  2-7.  SSM/T-2  weighting  functions  for  tropical  profile: 
(b)  cloudy. 
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Figure  2-8.  SSM/T-2  weighting  functions  for  midlatitude  summer  profile: 

(a)  clear,  (b)  cloudy. 
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Figure  2-9.  SSM/T-2  weighting  functions  for  midlatitude  winter  profile: 
(a)  clear,  (b)  cloudy. 
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models,  cloud  vertical  extent  had  to  be  interpolated  to  the  appropriate  pres¬ 
sure  level  within  individual  profiles.  Within  cloudy  layers  relative  humidity 
was  not  adjusted  to  saturation.  For  the  tropical  atmospheric  retrievals 
reported,  requiring  saturation  within  individual  model  layers  had  little 
effect  on  retrieval  results. 

The  impact  of  this  cloud  effect  on  calculated  brightness  temperatures  is 
illustrated  in  Figures  2-12  through  2-16,  Each  figure  is  annotated  by  cloud 
model  and  surface  emissivity  (ocean,  a,  and  land,  b) .  Plotted  is  channel 
brightness  temperature  vs.  integrated  liquid  water  content.  The  dashed  lines 
correspond  to  the  dependence  of  each  channel's  response  to  changes  in  liquid 
water  content  for  the  particular  cloud  chosen.  A  solid  horizontal  line 
appears  in  each  figure.  This  corresponds  to  a  brightness  temperature  equiva¬ 
lent  to  ♦‘he  physical  emission  temperature  of  the  cloud  top.  This  is  the 
temperature  which  would  be  sensed  by  an  infrared  radiometer,  i.e.,  that  cor¬ 
responding  to  an  optically  thick  cloud.  The  cloud  top  temperature  is  in  the 
upper  left  hand  corner  of  each  figure. 

Examining  these  results,  a  number  of  themes  emerge: 

(a)  For  low  clouds,  the  183  GHz  channels  are  virtually  unaffected  by 
increases  in  cloud  liquid  water.  This  is  apparent  from  the  flat  line 
behavior  of  the  184,32,  186.31,  and  190.31  GHz  curves.  This  does  not 
mean  that  there  is  no  cloud  impact  on  the  water  vapor  retrieval,  however. 

(b)  Although  the  183  GHz  channels  are  unaffected,  the  window  channels  at  90 
and  150  GHz  do  respond  to  the  cloud  water  presence.  Water  vapor 
retrieval  accuracy  is  highly  dependent  on  these  channels.  The  90  GHz 
channel  is  more  sensitive  to  cloud  than  the  150  GHz  channel. 

(c)  It  is  interesting  to  compare  the  figures  for  the  same  cloud  type  with 
different  backgrounds,  i.e.,  land  vs.  ocean.  There  is  virtually  no  dif¬ 
ference  in  the  signatures  of  the  183  GHz  channels  since  they  do  not  see 
the  surface.  The  90  and  150  GHz  channels  differ  considerably  depending 
on  background.  Over  the  ocean  these  channels  are  cold  (the  surface 
temperature  is  about  300  K) .  Over  land  it  is  a  few  degrees  cooler  than 
the  surface  temperature  due  to  water  vapor  effects.  The  presence  of  the 


cloud  introduces  a  relatively  warm  emission  source  over  the  ocean  and  a 
somewhat  cooler  one  (due  to  its  altitude  over  the  land).  Thus  increasing 
integrated  cloud  liquid  water  warms  the  signal  at  90  GHz  over  the  ocean 
and  cools  it  over  land.  Over  land,  behavior  is  the  same  as  that  in  the 
infrared  region  except  that  the  cloud  is  not  opaque. 

(d)  Sensitivity  to  cloud  presence  is  a  function  of  cloud  top  altitude,  opti¬ 
cal  thickness  (here  cloud  liquid  water  content) ,  and  the  underlying 
surface.  Over  the  ocean,  lower  clouds  are  more  detectable,  since  they 
are  warmer.  Over  land,  higher  clouds  are  more  detectable.  Lower  clouds 
look  radiometr ically  like  the  surface  (as  they  do  in  the  infrared) . 

(e)  As  the  cloud  liquid  water  content  increases,  the  brightness  temperatures 
at  90  and  150  GHz  asymtote  to  the  cloud  top  temperature.  Thus  with 
sufficient  water  content,  the  window  channel  signals  should  provide  an 
indication  of  cloud  top  height. 

(f)  For  higher  clouds  (see  the  cumulus  cloud  model  result),  a  response  from 
the  183  GHz  channels  is  evident.  This  is  manifested  as  a  cooling,  as  the 
cloud  blocks  emission  from  warmer  levels  below  the  cloud  top.  The  cumu¬ 
lus  cloud  quickly  becomes  optically  thick.  Over  the  ocean,  the  bright¬ 
ness  temperature  at  90  GHz  quickly  saturates  and  then  begins  to  cool. 
These  curves  resemble  brightness  temperature  vs.  rainfall  rate 
relationships  at  lower  frequencies. 

In  summary,  these  calculations  indicate  that  cloud  affects  millimeter 
wave  sensor  channels  and  that  there  is  potentially  information  on  cloud 
presence  and  height  in  the  brightness  temperature  data  obtainable  from  the 
millimeter  wave  moisture  channels. 

2.5.3  Water  Vapor  Profile  Retrievals  in  the  Presence  of  Clouds 


Vertical  moisture  profile  retrievals  were  obtained  from  independent  sam¬ 
ples  of  simulated  sensor  brightness  temperatures  using  a  statistical  inversion 
technique  (Rodgers,  1976;  Gaut  et  al.,  1975).  The  essential  element  of  the 
scheme  is  to  choose,  in  a  statistical  sense,  the  most  probable  combination  of 
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atmospheric  and  surface  properties  which  produces  the  set  of  measured  radio- 
metric  data  values.  It  is  a  general  regression  technique  which  minimizes  the 
mean  square  error  between  the  estimated  and  observed  values  of  the  parameter 
of  interest.  Here  the  desired  parameters  are  the  absolute  abundance  of  water 

.  9 

vapor  (molecules  cm  )  in  six  layers  corresponding  to  0-200,  200-300,  300  500 , 
500-700  mb,  and  850-1000  mb,  respectively .  Parameter  and  resultant  data  vec¬ 
tor  ensembles  are  divided  into  two  groups,  designated  A  and  B,  corresponding 
to  a  priori  and  retrieval  subsets,  respectively.  Water  vapor  profiles  are 
retrieved  from  the  B  data  vectors  based  on  statistics  developed  from  the  a 
priori  A  set.  An  estimate  of  the  parameter  vector  p  is  obtained  from  the 
simulated  data  vector  [i.e.  the  channel,  brightness  temperatures 
from: 

p  -  D  •  4,  (2.17) 

where  the  D-matrix  is  computed  from  the  a  priori  parameter/data  vector 
ensemble  A  using: 

D  -  C(P,^)C'1(^,^)  (2.18) 

where  p  is  the  parameter  vector  from  the  a  priori  set,  ^  is  a  generalized  data 
basis  function  defined  as 

,  d, -  <d,>  d  -  <d  > , 

4>  «  column  ■••••  <d'>  (2.19) 

1  1  m  1 

for  m  channel  radiometric  data  d  in  the  dependent  set  A  (d^ ,  i  -  1 . m) ,  the 

brackets  denote  the  mean  value  operator,  and  C(<f>,<f>')  and  C(p  ,<f>)  are  the 
diagonalized  autocorrelation  and  crosscorrelation  matrices  of  (<f>,4>)  and  (p,^) 
respectively.  For  the  SSM/T-2  channel  set  specified  in  Table  2-1,  m«5 .  The 
effect  of  cloud  on  water  vapor  profile  retrieval  accuracy  is  assessed  by  in¬ 
cluding  cloud  in  the  data  vector  simulation.  Independent  retrieval  statistics 
were  then  evaluated  for  clear  and  cloudy  situations.  Thus,  clear  cases  are 
retrieved  using  a  "clear"  D  matrix  and  cloudy  cases  using  a  "cloudy"  D 
matrix.  This  approach  avoids  the  application  of  nonrepresentat i ve  statistics 
for  either  clear  or  cloudy  cases.  It  should  be  noted,  however,  that  in  prac- 
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Figure  2-14.  Brightness  temperature  (K)  versus  integrated  liquid  water 

content  (g/cm  )  over:  (a)  ocean  and  (b)  land  for  altostratus 
cloud. 
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Figure  2-15.  Brightness  temperature  (K)  versus  integrated  liquid  water 

content  (g/cm  over:  (a)  ocean  and  (b)  land  for  strat ocumulus 
cloud . 
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Figure  2-16.  Brightness  temperature  (K)  versus  integrated  liquid  water 

content  (g/cm  )  over:  (a)  ocean  and  (b)  land  for  nimbostratus 
cloud. 


tice  it  may  not  always  be  possible  to  identify  clear  and  cloudy  cases  from  the 
data  itself.  Retrieval  accuracy  was  assessed  by  comparing  inferred  layer 
water  vapor  abundances  to  those  in  the  actual  profiles  (i.e.  those  used  in  the 
B  set  simulations)  and  then  evaluating  the  fractional  root  mean  square  (RMS) 
error  over  the  ensemble  of  retrievals.  The  RMS  fractional  error  for  each 
layer  k  evaluated  for  the  set  of  N  -  100  independent  soundings  was  defined  as: 


RMS (k)  - 


1 

<u(k)> 


(u(k,j) 


J-l 


u(k, j) 


1/2 


(2.20) 


where  u(k,j)  and  u(k,j)  are  the  retrieved  and  actual  water  vapor  amounts  for 
the  kth  layer  and  j th  sounding,  respectively,  and  <u(k)>  is  the  layer  mean 
value.  For  comparison,  the  same  statistic  was  evaluated  assuming  the  mean  of 
the  ensemble  retrieval  for  each  sounding,  i.e.  by  replacing  u  by  <u>  itself. 
This  provides  a  measure  of  the  RMS  fractional  error  obtained  when  only  cli¬ 
matological  mean  estimates  are  available. 
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A  comparison  of  statistical  water  vapor  profile  retrieval  results  for 
clear  and  cloudy  cases  is  illustrated  in  Figures  2-17  and  2  18.  Plotted  are 
layer  fractional  RMS  errors  as  a  function  of  mean  layer  pressure  for  clear  and 
cloudy  tropical  soundings  over  oceanic  (Figure  2-17)  and  land  (Figure  2-18) 
backgrounds.  Each  figure  also  shows  the  RMS  fractional  error  obtained  when 
climatology  (i.e.  the  layer  mean  value)  is  used  in  Eqn.  (2.20)  as  the  re¬ 
trieved  layer  value.  This  result  is  related  to  the  climatological  variance  of 
the  sample  and  provides  a  measure  of  retrieval  skill.  The  ratio  of  the  re¬ 
trieval  RMS  fractional  error  to  that  of  the  climatology  curve  in  each  figure 
is  a  measure  of  the  remaining  unexplained  variance  after  the  retrieval  pro¬ 
cess.  Thus  small  values  of  this  ratio  are  deemed  skillful  retrievals  and  no 
skill  is  demonstrated  as  this  ratio  approaches  unity. 

It  can  be  seen  that  the  selected  millimeter  wave  channel  set  should 
provide  water  vapor  retrievals  which  are  quite  skillful  in  clear  situations 
over  ocean  backgrounds,  particularly  in  the  middle  troposphere.  Over  land, 
the  retrievals  are  still  skillful  although  improvement  over  climatology  is  not 
quite  as  distinctive.  The  enhanced  retrieval  skill  over  the  low  emissivity 
oceanic  background  can  be  extrapolated  to  other  low  emissivity  surfaces  such 
as  dry  snow  which  provide  high  contrast  for  atmospheric  emission. 

The  clear  atmosphere  retrieval  results  can  be  contrasted  with  those 
obtained  in  cloudy  situations.  The  fractional  RMS  error  of  the  simulated 
retrievals  over  both  oceanic  and  land  backgrounds  is  increased  considerably 
due  to  the  presence  of  beam  filling  cloud.  The  effect  is  particularly 
noticeable  over  the  oceans  where  the  clear  retrievals  appear  to  be  quite 
good.  Notably,  at  upper  tropospheric  levels  (pressures  below  about  600mb) , 
cloudy  retrieval  accuracies  are  about  the  same,  regardless  of  background. 

Thus  it  appears  that  the  presence  of  an  obscuring  cloud  over  the  ocean  reduces 
the  apparent  advantage  of  millimeter  wave  sounding  against  the  comparatively 
low  emissivity  ocean  surface  in  clear  situations. 

While  the  retrieval  approach  adopted  here  to  demonstrate  the  effect  of 
clouds  is  quite  simple,  it  is  analogous  to  that  proposed  for  operational 
implementation  with  SSM/T-2  millimeter  wave  sounder  data  (Stogryn,  personal 
communication).  Within  the  context  of  such  statistical  retrievals,  the 
treatment  of  cloud  is  problematical.  Given  the  potential  degradation  of  water 
vapor  profile  retrievals  due  to  cloud  presence,  a  method  which  explicitly 
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treats  cloud  within  the  retrieval  algorithm  may  be  usefully  applied.  Physi¬ 
cally  based  approaches  used  for  infrared  sounder  data  may  provide  appropriate 
treatments . 

2 . 6  A  Unified  Retrieval  Methodology  for  the  DMSP  Meteorological  Sensors 

The  potential  for  improving  statistical  retrievals  of  temperature  and 
moisture  profiles  using  a  physical  retrieval  was  demonstrated  in  Sections  2.3 
and  2.4  In  Section  2.5,  the  degradation  of  moisture  retrieval  results  due  to 
cloud  presence  was  illustrated.  Here  we  propose  a  method  to  incorporate  all 
DMSP  meteorological  sensor  data  including  that  from  the  SSM/T,  SSM/T-2,  SSM/I , 
and  OLS  (see  Table  2-1)  into  a  unified  retrieval. 

2.6.1  Proposal  Scheme 


The  unified  retrieval  approach  for  the  DMSP  illustrated  in  Figure  2-19 
employs  physical  considerations  and  allows  for  the  incorporation  of  all  data 
sources.  Recognizing  potential  operational  constraints,  an  attempt  has  been 
made  to  build  on  the  attributes  of  the  existing  DMSP  retrieval  capability  and 
experience.  The  microwave  sensor  data,  T^0,  Is  employed  with  the  "D"  matrix 


Figure  2-19.  Schematic  of  unified  DMSP  retrieval  scheme. 


statistical  retrieval  to  provide  first  guesses  for  the  desired  parameters, 

P°:  temperature,  T(p),  and  water  vapor  profiles,  q(p),  surface  emissivity, 

<s,  and  temperature,  T  .  Simulated  brightness  temperatures,  T^n,  are  then 
evaluated  to  examine  the  consistency  of  the  first  guesses  with  the  observa¬ 
tions.  The  forward  problem  calculation  (denoted  by  F  in  Fig.  2-19)  is  ac¬ 
complished  using  the  RADTRAN  simulation  code  (Falcone  et  al,  1982)  as  modified 
by  Isaacs  et  al  (1985a).  When  the  residuals,  An,  (i.e.  differences  between 
simulated  and  observed  brightness  temperatures)  are  small,  the  process 
terminates.  However,  when  residuals  are  larger  than  a  preset  tolerance  (usu¬ 
ally  determined  by  the  sensor  noise  equivalent  brightness  temperatures  (NEAT) 
and  scene  noise),  the  procedure  goes  on  to  adjust  the  first  guess  profiles. 
This  adjustment  is  accomplished  by  using  the  residuals  in  a  simultaneous 
physical  retrieval. 

Cloud  or  precipitation  in  the  field-of -view  (FOV)  of  the  microwave  sen¬ 
sors  can  be  problematic.  Precipitation  will  generally  preclude  soundings  of 
temperature  and  moisture  and  the  determination  of  surface  properties.  Quality 
control  flags  for  precipitation  (as  well  as  precipitation  amounts)  can  be  ob¬ 
tained  from  the  SSM/I  statistical  retrieval  (Lo,  1983;  Jin  and  Isaacs,  1987). 
Isaacs  and  Deblonde  (1987)  have  discussed  the  potential  impact  of  cloud  on 
statistical  millimeter  wave  water  vapor  retrievals  and  evaluated  the  sensi¬ 
tivity  of  these  channels  to  cloud  presence.  Cloud  fields  from  the  DMSP  OLS 
(using  an  appropriately  spatially  averaged  subset  of  visible  and  infrared 
imagery)  aid  in  cloud/no  cloud  discrimination.  To  determine  first  guess  cloud 
properties  within  the  FOV  necessary  to  accomplish  the  physical  retrieval  step, 
the  high  spatial  resolution  OLS  imager  data  is  utilized.  Image  processing  of 
this  data  within  the  relatively  larger  microwave  footprint  provides  first 
guess  cloud  coverage,  N,  and  equivalent  brightness  temperatures  (EBTs)  for 
cloud  top  and  surface,  T^c,s.  Cloud  top  brightness  temperature  along  with  the 
first  guess  temperature  profile  yields  a  first  guess  cloud  top  pressure,  pc . 
Over  the  oceans,  cloud  properties  derived  from  the  OLS  imager  data  are  supple¬ 
mented  by  information  on  cloud  integrated  liquid  water  content  (ILWC)  avail¬ 
able  from  the  SSM/I.  Cloud  ILWC  provides  a  parameterization  of  cloud  optical 
thickness  and  emissivity.  These  cloud  properties  are  required  to  treat  the 
effect  of  cloud  on  the  SSM/T-2  sensor  data  and  therefore  are  input  to  the 
forward  problem. 
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In  cloudy  areas,  cloud  coverage  and  cloud  top  height  first  guesses  are 
necessary  for  the  physical  retrieval  step  (Equation  2.1  is  modified.)  Co- 
located  with  the  microwave  sensors  aboard  the  DMSP  spacecraft,  the  Operational 
Linescan  System  (OLS)  provides  both  visible  and  infrared  imagery  at  high 
spatial  resolution.  With  much  higher  spatial  resolution,  the  visible  and 
infrared  data  from  the  OLS  imagery  can  be  used  to  characterize  the  uniformi  ty 
of  the  much  larger  microwave  footprints.  In  those  areas  where  the  contribu¬ 
tions  from  the  atmosphere  to  microwave  brightness  temperature  are  small  (i.e., 
nonprecipitating  situations) ,  visible  or  infrared  data  is  able  to  provide 
guidelines  on  the  uniformity  of  the  surface  observed  within  a  field-of -view. 
When  clouds  obscure  portions  of  the  microwave  field-of-view,  the  imager  data 
provides  the  complementary  capability  of  cloud  property  determination. 

Classically,  techniques  to  infer  FOV  non-uniformity  have  been  referred  to 
as  texture  analysis  methods,  A  number  of  approaches  can  be  used  for  texture 
analysis  including:  (a)  examination  of  the  spatial  power  spectrum  of  an  im; ge 
through  Fourier  decomposition,  (b)  edge  enhancement,  and  (c)  spatial  co¬ 
herence.  We  have  chosen  the  spatial  coherence  approach  (Coakley  and  Brether- 
ton,  1982)  for  the  determination  of  both  cloud  and  surface  properties  from  OLS 
data.  The  statistics  evaluated  are  the  local  mean  (I)  and  local  standard 
deviation  (LSD)  of  radiance  (or  gray  shade)  values.  The  LSD  is  calculated  for 
n  x  n  sets  of  pixels. 

A  plot  of  LSD^  vs.  1^  gives  the  cloud  coverage  fraction  within  the  micro - 
wave  footprint  and  the  EBT  of  the  surface  and  effective  cloud  top.  Figures  2- 
20a, b  illustrate  LSD^  vs.  1^  frequency  plots  for  a  partially  cloudy  microwave 
FOV  over  the  ocean  evaluated  using  GOES  visible  and  infrared  imagery.  The 
spike  in  the  visible  result  (Fig.  2-20a)  denotes  the  surface  reflectance  value 
due  to  clear  visible  pixels.  The  higher  reflectivity,  signatures  with  nonzero 
standard  deviation  result  from  partially  cloudy  pixels.  In  the  infrared  data 
(Fig.  2-20b),  the  highest  and  lowest  EBTs  correspond  to  the  surface  and  cloud 
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top  emission,  respectively.  The  cloud  top  EBT  is  used  with  the  first  guess 
temperature  profile  to  obtain  the  first  guess  effective  cloud  top  pressure, 
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Figure  2-20.  Spatial  coherence  results  for  partially  cloudy  microwave 
field-of -view  using:  (a)  visible  and  (b)  infrared  GOES 
imager  data. 


3.  TASK  2  -  SIMULATION  MODEL  ENHANCEMENTS 

Four  subtasks  were  completed  during  the  first  year's  effort  to  enhance 
the  simulation  model  performance.  These  are  described  in  this  section. 

3 . 1  Multiple  Scattering 

The  multiple  scattering  submodel  to  treat  large  cloud  droplets  and  liquid 
and  glaciated  precipitation  used  in  the  Phase  I  effort  was  based  on  the  dis¬ 
crete  ordinate  method  as  described  by  l.iou  (1973)  and  l.iou  et  al .  (1980).  In 
this  approach  the  scalar  radiative  transfer  equation  for  brightness  tempera¬ 
ture  within  the  scattering  cloud  is  reduced  by  Gaussian  quadrature  of  degree  n 
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(in  fact  the  method  is  referred  to  in  some  circles  as  the  Gaussian  quadrature 
method)  to  a  set  of  2n  coupled  differential  equations: 


Ui^Ti(r)  -  T . ( r )  -  X  (  C.  .T  .+  C.  .T.)  -  (1-v)  T 
dr  i  ,  j  i.-J  *J  x.J  J 


(3.1) 


where  is  the  brightness  temperature  at  the  i th  Gauss  zenith  angle  cosine  u 


i  2n+ 1 

C  -  1/2  E  a  (21+l)u  P  (M  )P  (M.) 
X  •  J  i  J  i  I  1  i  j 


(3.2) 


In  the  expression  above,  the  a^-e  "‘-.efficients  of  the  Legendre  polynomial 

expansion  of  the  phase  function  and  the  aj  are  the  Gaussian  quadrature  weights 

(a.  -  a  i  ;  n .  —  -p  )  .  The  solution  for  each  is  given  by: 

J  J  J  J 


T.—  Z  L.g.(u.)  exp  (-k.r)  +  T 
1  ■  J  J  J  J  c 

j«-n  J  J  J  J 


(3.3) 


where  the  gj  and  the  kj  are  the  eigenvectors  and  eigenvalues  of  the  system  of 

equations,  respectively,  and  T  is  the  cloud  temperature.  The  L-’s  are  con- 

c  J 

stants  determined  by  the  boundary  conditions  of  the  problem. 

Liou  (1973)  originally  discussed  a  variety  of  numerical  stability  prob¬ 
lems  associated  with  his  application  of  the  polynomial  root  searching  method 
and  subsequent  evaluation  of  the  corresponding  eigenvectors.  This  method  was 
also  implemented  within  the  microwave  program  supplied  to  AFGL  (Liou  et  al., 
1980) .  We  implemented  a  new  method  to  evaluate  eigenvalues  and  eigenvectors 
based  on  available,  standardized  algebraic  eigensolution  methods.  The  ad¬ 
vantage  of  these  approaches  is  that  they  increase  accuracy  and  decrease  the 
time  required  to  find  the  solutions  to  these  systems  of  equations.  This  ap¬ 
proach  has  already  been  tried  to  solve  analogous  equations  for  the  scattering 
of  solar  radiation  (Stamnes  and  Swanson,  1981;  Stamnes  and  Conklin,  1984) 
using  algorithms  available  from  the  IMSL  (1975)  library. 

Tills  approach  was  implemented  within  the  discrete  ordinate  multiple  scat¬ 
tering  code  for  thermal  radiation  using  matrix  eigensystem  routines  from  the 
EISPACK  (Smith  et  al.,  1975)  program  set  available  on  AER's  Harris  H800  compu¬ 
ter.  Table  3-1  compares  the  eigenvalues  obtained  using  both  Liou's  root 
searching  method  and  the  EISPACK  code.  These  were  generated  for  a  case  treat 
ing  the  evaluation  of  brightness  temperatures  at  50.5  GHz  looking  at  a  rain- 
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rate  of  5  mm/h  using  16  streams  (eight  up  and  eight  down) .  Both  positive  and 
negative  eigenvalues  are  used  for  a  total  of  16  roots.  The  root  searching 
calculations  had  to  be  performed  in  quad  precision  on  the  Harris  (approximate¬ 
ly  double  precision  on  other  machines)  in  order  to  insure  convergence.  By  us¬ 
ing  the  EISPACK  eigensystem  solution  in  place  of  the  root  searching  methods, 
accuracies  were  matched  and  exceeded  and  a  time  saving  of  about  a  factor  of 
ten  was  realized  in  the  section  of  the  discrete  ordinate  code  which  evaluated 
the  eigenvalues  and  eigenvectors. 

Aside  from  the  method  of  eigensystem  solution  this  code  operates  identi¬ 
cally  to  and  has  the  same  restrictions  as  that  used  in  the  Phase  I  study.  For 
example,  it  does  not  treat  polarization  (see  section  3.3)  and  looks  only  at 
nadir.  Additionally,  the  code  must  be  run  along  with  the  attachment  of  appro¬ 
priate  subroutines  from  the  EISPACK  library. 


Table  3-1.  Comparison  of  Eigenvalues  From  Liou's  Root  Searching  Method  for 
the  Discrete  Ordinate  Method  and  those  from  the  EISPACK  Routines 


Eigenvector  # 

Liou 

EISPACK 

1 

0.948 

0.947994 

2 

1.030 

1.029818 

3 

1.116 

1.116219 

4 

1.274 

1.273931 

5 

1.552 

1.552476 

6 

2.089 

2.088644 

7 

3.390 

3.389966 

8 

10.04 

10.036796 

3.2  Mie  Scattering  Parameters 


The  discrete  ordinate  code  discussed  in  the  previous  section  requires, 
among  its  various  inputs,  the  Mie  scattering  parameters  for  the  particular 
frequency  and  rainfall  droplet  size  distribution  desired.  Unfortunately,  this 
data  is  not  calculated  on-line  with  the  program  and  must  he  supplied  by  the 
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user.  In  order  to  circumvent  this  inconvenience,  a  subroutine  was  written 
which  provided  the  necessary  input  data  on  line  with  specification  of  a  few 
simple  meteorological  parameter  inputs  such  as  rainfall  rate,  desired  size 
distribution  (Marshall - Palmer  or  Best),  and  phase  of  the  precipitation  (i.e., 
water  or  ice).  This  parameterization  is  based  on  the  published  values  of 
Savage  (1978)  and  is  briefly  discussed  below  and  in  Isaacs  et  al.  (1988a). 

The  algorithm  is  valid  between  the  frequencies  of  19.35  and  230  GHz  and  for 
the  temperature  range  between  263  and  283  K.  An  algorithm  suitable  for 
completing  exact  Mie  theory  calculations  is  also  described. 

3.2.1  Precipitation  Scattering  Property  Parameterization 

The  RADTRAN  computer  code  (Falcone  et  al.,  1982)  was  developed  by  the  Air 
Force  Geophysics  Laboratory  (AFGL)  to  provide  atmospheric  attenuation  and 
brightness  temperature  calculations  for  typical  atmospheric  paths  over  the 
frequency  range  from  1  to  300  GHz.  RADTRAN  provides  a  design  tool  which  can 
readily  be  used  to  assess  potential  environmental  impacts  on  microwave  and 
millimeter  wave  sensors.  The  atmospheric  attenuation  submodels  of  the  clear 
atmosphere,  fog,  cloud,  and  rain  used  in  RADTRAN  have  been  thoroughly  docu¬ 
mented  (Falcone  et  al.,  1979).  The  models  for  fog,  cloud,  and  rain  attenua¬ 
tion  have  also  been  incorporated  into  the  commonly  used  FASCODE  accelerated 
line-by-line  code  (Clough  et  ai.,  1986).  These  codes  can  be  used  directly  to 
evaluate  atmospheric  brightness  temperatures  in  the  absence  of  significant 
contributions  to  path  radiance  due  to  multiple  scattering.  This  is  generally 
the  case  for  non-precipitating  environments . 

When  precipitation  is  present,  however,  evaluation  of  the  multiply  scat¬ 
tered  brightness  temperature  field  requires  both  a  generalization  of  the  radi¬ 
ance  solution  and  appropriate  supplementary  data  on  precipitation  scattering 
properties.  To  generalize  the  radiance  solution,  the  multiple  scattering 
version  of  RADTRAN  uses  a  numerical,  Gaussian  quadrature  approach  to  the 
radiative  transfer  equation  (Jin  and  Isaacs,  1985).  The  FASCODE  approach  is 
based  on  an  approximation  of  the  multiple  scattering  source  function  (Isaacs 
et  al.,  1987).  Both  treatments  require  the  specification  of  precipitation 
scattering  optical  properties  including  the  extinction  coefficient  (km ’^), 
single  scattering  albedo,  and  the  angular  scattering  function.  These  data  are 
not  currently  implemented  within  RADTRAN  and  must  be  supplied  by  the  user 
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Precipitation  scattering  properties  are  generally  available  via  standard 
Mie  theory  calculations.  The  Mie  theory  formalism  requires  a  knowledge  of 
particle  size  distribution  and  index  of  refraction  (Falcone  et  al.,  1979). 

The  index  of  refraction,  in  turn,  is  dependent  on  frequency,  phase  (i.e.  ice 
or  water),  and  temperature  To  avoid  the  cumbersome  necessity  of  performing 
on-line  Mie  theory  calculations  to  support  each  possible  combination  of  these 
model  variables  within  multiple  scenario  brightness  temperature  simulations,  a 
parameterization  has  been  developed  based  on  the  existing  Mie  theory  calcula¬ 
tions  of  Savage  (Savage,  1978).  This  parameterization  is  available  for  imple¬ 
mentation  within  both  RADTRAN  and  FASCODE. 

The  resultant  subroutine  provides  an  efficient  method  to  obtain  the  ex¬ 
tinction  coefficient,  single  scattering  albedo,  and  angular  scattering  func¬ 
tion  over  the  frequency  domain  from  19  to  290  GHz.  The  angular  scattering 
function  is  given  in  terms  of  its  first  eight  Legendre  polynomial  expansion 
coefficients.  Precipitation  angular  scattering  functions  are  not  highly 
anisotropic  at  microwave  frequencies  and  this  number  of  terms  usually  suffices 
to  describe  them.  Furthermore,  this  number  of  terms  is  consistent  with  the 
Gaussian  quadrature  required  to  specify  the  brightness  temperature  field  (Jin 
and  Isaacs,  1985;  Savage,  1978).  The  scattering  function  asymmetry  factor 
used  in  the  FASCODE  multiple  scattering  approximation  is  easily  obtained  from 
a  knowledge  of  the  second  Legendre  coefficient. 

Input  parameters  to  the  subroutine  consist  of  frequency  (GHz),  rainfall 
rate  (mmhr  *) ,  temperature  (deg  C) ,  phase  (i.e.  water  or  ice),  and  size 
distribution.  The  subroutine  allows  the  user  to  choose  either  the  Marshall- 
Palmer  or  Best  size  distribution  (Savage,  1978).  Details  of  the  selected  size 
distribution  are  varied  with  rainfall  rate;  higher  rainfall  rates  skewing  the 
distributions  to  larger  drop  sizes.  Various  checks  are  made  on  the  input  data 
to  insure  internal  consistency  with  the  stored  data  domain.  A  binary  search 
routine  is  executed  and  a  simple  interpolation  is  performed  to  provide  desired 
scattering  properties  at  intermediate  temperatures  and  frequencies  not  covered 
by  the  internal  data  set. 

We  have  compared  scattering  parameters  obtained  using  this  subroutine 
with  the  results  of  exact  Mie  theory  calculations.  Comparisons  were  done  at 
millimeter  wavelengr  -ince  precipitation  angular  scattering  functions  are 
more  pronounced  at  t . ,  e  frequencies  ami  thus  provides  a  worse  case  assess- 
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ment .  Fig.  3-1  compares  extinction  coefficient  (BEXT)  and  single  scattering 
albedo  (OMEGAO)  calculated  by  the  subroutine  (lines)  and  exact  Mie  theory 
calculations  (symbols)  throughout  the  millimeter  wave  region  from  89.5  to  183 
GHz  at  two  values  of  the  rainfall  rate  (2  and  5  mmhr "^).  Fig.  2  illustrates 
much  the  same  comparison  at  selected  frequencies  for  the  first  eight  Legendre 
polynomial  coefficients  of  a  Marshal 1  -  Palmer  size  distribution  of  rain  drops 
with  a  rainfall  rate  of  5  mmhr  ^ .  In  general  the  subroutine  produces  scat¬ 
tering  properties  quite  comparable  to  the  more  computationally  demanding  exact 
Mie  theory  calculations . 

To  assess  the  accuracy  of  the  precipitation  scattering  properties  for 
RADTRAN  brightness  temperature  simulations,  calculations  were  performed  using 
an  exact  Gaussian  quadrature,  numerical  solution  to  the  multiple  scattering 
radiative  transfer  equation  (Jin  and  Isaacs,  1985)  using  both  exact  Mie  theory 
data  and  those  obtained  from  the  subroutine  described  above.  Brightness  tem¬ 
perature  differences  were  at  most  one  or  two  tenths  of  a  degree  over  the  fre¬ 
quency  range  from  89.5  to  183  GHz.  This  is  less  than  the  noise  equivalent 
brightness  temperature  achievable  by  radiometers  operating  in  this  frequency 
range.  The  subroutine  thus  provides  prec i pi tation  scattering  properties  which 
are  of  considerable  practical  value  for  brightness  temperature  simulation 
purposes . 

When  less  accurate  brightness  temperature  field  simulations  are  adequate, 
considerable  computer  time  can  be  saved  '  y  using  an  approximate  multiple  scat¬ 
tering  approach.  We  have  compared  exact  numerical  results  (Jin  and  Isaacs, 
1985)  and  those  obtained  using  the  FASCODE  multiple  scattering  approximation 
(Isaacs  et  al.,  1987).  The  model  consisted  of  a  3  km  rain  layer  with  rainfall 
rate  of  15  mmhr  .  Optical  depths  and  single  scattering  albedos  at  19.35  and 
37  GHz  were  (0.91,  0.21)  and  (3.20,  0.39),  respectively.  The  surface  tem¬ 
perature  was  assumed  to  be  287.0  K  and  the  rain  layer  was  nonisothermal  with  a 
mean  layer  temperature  of  271.0  K  and  temperature  gradient  of  19.5  K  over  the 
3  km  layer.  Results  of  the  comparison  are  illustrated  in  Fig.  3-3  as  a 
function  of  viewing  angle.  Brightness  temperature  differences  between  exact 
(lines)  and  approximate  (symbols)  results  are  generally  a  few  degrees  except 
near  the  horizon  (i.e.  90  degrees).  This  behavior  is  due  to  the  multiple 
scattering  approximation  which  is  more  accurate  for  angles  near  the  zenith. 
Even  at  the  horizon,  the  percent  error  for  the  37  GHz  case  (Fig.  3  -  3b )  is 
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about  5%.  This  is  well  within  the  20%  RMS  error  cited  for  the  approximation 
(Isaacs  et  al.,  1987). 

This  approach  should  efficiently  provide  the  required  microwave 
scattering  properties  of  precipitation  with  sufficient  accuracy  to  support 
brightness  temperature  simulation  calculations  such  as  those  in  the  multiple 
scattering  version  of  RADTRAN . 

3.2.2  Exact  Mle  Scattering  Algorithm 

For  situations  not  represented  by  the  parameterization  discussed  above, 
it  is  necessary  to  perform  the  Mie  theory  calculations  for  the  desired  condi¬ 
tions  to  obtain  the  necessary  input  data  for  the  discrete  ordinate  method 
multiple  scattering  code.  Specifically,  the  data  necessary  are  the  extinction 
coefficient,  the  albedo  of  single  scattering,  and  the  angular  scattering 
function.  The  angular  scattering  function  is  input  as  coefficients  w^  of  the 
normalized  Legendre  polynomial  series,  i.e.: 

N 

P(0)  -  S  v  P  (6)  (3.4) 

1-1 

where  the  are  the  Legendre  polynomials.  The  number  of  terms  required  is 
equal  to  the  number  of  streams  desired  in  the  brightness  temperature  calcula¬ 
tion.  For  example,  a  16  stream  calculation  will  require  sixteen  coefficients 
and  yield  a  phase  function  with  an  accuracy  of  N-16  terms  in  the  summation  of 
equation  3.4. 

In  order  to  provide  the  capability  to  generate  the  required  Mie  paramet¬ 
ers  for  arbitrary  size  distributions,  temperatures,  and  updated  complex  index 
of  refraction  data,  a  set  of  computer  codes  has  been  modified  to  run  on  the 
AFGL  Cyber.  These  codes  called  RADI  and  RAD2  are  taken  from  the  algorithms 
developed  by  Dave  (1972).  The  first  program,  RADI,  calculates  the  optical 
parameters  and  coefficients  of  the  Legendre  series  representing  the  scattering 
properties  of  a  single  spherical  particle.  The  second  program,  RAD2 ,  calcu¬ 
lates  the  analogous  properties  for  a  chosen  size  distribution  of  particles. 

The  size  distribution  can  be  chosen  from  among  the  discontinuous  power  law 
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Comparison  of  extinction  coefficient  (BEX!)  and  single  scattering 
albedo  (OMEGAO)  calculated  by  the  subroutine  (lines)  and  exact  Mie 
theory  calculations  (symbols)  throughout  the  millimeter  wave  region 
from  89.5  to  183  GHz  at  two  values  of  the  rainfall  rate  (2  and  5 
mmhr~  ) . 


N 

Comparison  of  the  first  eight  Legendre  polynomial  coefficients  (N) 
of  a  Marshall-Palmer  size  distribution  of  rain  drops  with  a  rainfall 
rate  of  5  mmhr calculated  by  the  subroutine  (lines)  and  exact  Mie 
theory  calculations  from  89.5  to  183  GHz. 


Fig.  3-3.  Comparison  between  exact  (lines)  and  approximate  (symbols)  bright¬ 
ness  temperature  calculations  as  a  function  of  viewing  angle  (THETA) 
for  a  3  km  rain  layer  with  rainfall  rate  of  15  mmhr'  at  (a)  19.35 
GHz  and  (b)  37.0  GHz. 


distribution,  the  modified  gamma  distribution,  and  the  log-normal  distribu¬ 
tion.  The  modified  gamma  distribution,  for  example,  given  by: 


n(r)  -  a  rQ  exp  (-  b  r^) 


(3.5) 


can  be  used  to  express  the  commonly  used  Marshall - Palmer  distribution  for 


precipitation : 


n(r)  -  NQ  exp  (  -b  r  ) 


(3. A) 


where  NQ  -  8.0  x  106  m'4  and  b  -  8200  R’®'^,  and  R  is  the  rainfall  rate  in 
mm/h.  The  log-normal  distribution  is  convenient  to  use  to  fit  experimental 
drops  size  data. 

Table  3-2  summarizes  the  file  structure  for  the  RADI  and  RAD2  programs 
and  their  relationship  to  one  another  For  RADI,  a  single  input  file  is 
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required  and  a  single  output  file  results.  The  output  file  from  RADI  is  used 
as  one  input  file  to  RAD2.  A  second  RAD2  input  file  specifies  the  parameters 
defining  the  desired  size  distribution. 


Table  3-2.  File  Structure  Relating  RADI  and  RAD2  Programs 


|  RAD1IN  | 


Tape  15 


RADI 


RAD10UT  | 


Tape  21 


|  RAD10UT  |-- 


Tape  15 


RAD2 


|  RAD2IN  | 


|  RAD20UT  | 


Tape  22 


Tape  21 


The  input  file  RAD1IN  is  a  single  record  which  specifies  the  real  and 
imaginary  parts  of  the  complex  index  of  refraction,  n]  and  ,  respectively, 
the  initial  Mie  parameter,  x,  the  Mie  parameter  increment,  Sx,  the  total 
number  of  x's  desired,  Nx,  and  an  integer  specifying  the  desired  program 
option,  NOPT  (—1  or  2).  To  obtain  the  desired  output  file  for  RAD2  the  option 
NOPT  —  1  is  selected.  The  input  file  and  format  is  illustrated  in  Table  3-3. 


Table  3-3.  RAD1IN  Input  File  and  Format 


Card  1 


N1  N2  XI  DE1.X  NX  NOPT 
D15.5  D15.5  D15.5  D15.5  15  15 


The  Mie  parameter  x  is  a  dimensionless  function  of  the  droplet  radius,  r,  and 
wavelength,  A,  given  by: 


x  -  2*r  /  A 


(3.7) 


The  specification  of  the  input  file  for  RADI  is  not  independent  of  that  for 
RAD2  which  requires  the  upper  and  lower  limits  in  radius  to  the  integration  of 
optical  properties  over  the  desired  size  distribution.  The  range  of  Mie  para¬ 
meters  calculated  bv  RADI  must  at  least  span  that  asked  for  in  the  integration 
of  the  RAD2  program.  For  example,  if  the  desired  wavelength  of  the  calcula¬ 
tion  is  A,  and  the  upper  and  lower  radius  limits  are  r?  and  ri ,  respectively, 
the  initial  Mie  parameter  value,  x,  specified  in  the  RAD1IN  file  must  be  less 
than  or  equal  to: 


x  -  2irr  /  A 
min  2  ' 


(3.8) 


Furthermore  the  largest  Mie  parameter  calculated  by  RADI  xmax  must  be  greater 
than  or  equal  to: 


x_„„  -  x  +5x  *  Nx  -  2xr  /A 
max  i 


(3.9) 


The  main  input  file  to  PAD2  in  the  output  file  of  RADI  which  we  have 
called  RAD10UT.  This  file  contains  Nx  records  corresponding  to  the  Mie  effi¬ 
ciency  factors  for  extinction,  absorption,  and  scattering  and  the  Legendre 
polynomial  expansion  coefficients  for  each  size  specified  in  the  input  to  RADI 
in  the  interval  from  x  to  x  +  6x*Nx.  The  other  input  file  RAD2IN  consists  of 
three  records  illustrated  in  Table  3-4  for  the  modified  gamma  distribution 
(see  equation  (3.5)). 


Table  3-4.  RAD2IN  Input  File  and  Format 


Card 

1 

a 

a 

b 

7 

blank 

2 

Card 

2 

X 

rmin 

rmax 

Card 

3 

blank 

D15.5 

D15.5 

D15 . 5 

D15.5 

15 

15 

The  parameter  JDIS  determines  the  type  of  size  distribution  selected  with 
values  1,  2,  and  3  corresponding  to  the  discontinuous  power  law  distribution, 
modified  gamma  function,  and  log-normal  distribution,  respectively.  Specifi¬ 
cations  for  the  other  size  distribution  options  can  be  found  in  Dave  (1972). 

Output  from  RAD2  is  illustrated  in  Table  0  9.  Two  calculations  are 
illustrated,  both  for  Marshall- Palmer  size  distributions  with  5  mm/h  of  rain 
(see  equation  (3.6)).  The  results  at  the  top  of  the  table  are  for  a  frequency 
of  89.5  GHz  (i.e.,  0.335  cm)  while  those  at  the  bottom  are  for  183  GHz  (i.e., 
0.164  cm).  The  values  numbered  from  1  to  N  concluding  each  output  are  the 
coefficients  of  the  Legendre  polynomial  expansion  of  the  phase  function, 

(see  equation  (3.4)).  The  first  16  of  these  coefficients  would  be  required  as 
input  for  a  sixteen  stream  run  of  the  discrete  ordinate  multiple  scattering 
code . 

3 . 3  Thermal  Radiative  Transfer  for  Polarized  Inhomogeneous  Nonisothermal 
Atmospheric  Precipitation  and  Statistical  Retrieval  of  Related 
Parameters 


3.3.1  Introduction 


A  numerical  approach  has  been  developed  to  treat  the  polarized 
thermal  radiative  transfer  theory  applicable  to  inhomogeneous  (stratified), 
nonisothermal  layers  composed  of  different  spherical  Mie  scatterers.  Calcula¬ 
tions  are  performed  to  obtain  the  microwave  brightness  temperatures  at  the 
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Table  3-5.  Sample  Output  Flies  from  the  RAD2  Code  for  a  Marshall  -  Palmer  Rain 
Distribution  at  5  mm/h  for  89.5  (top)  and  183  (bottom)  GHz 
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SSM/'i  frequencies  19.35,  22.235,  37,  and  85.5  GHz  in  the  presence  of 
atmospheric  precipitation  with  a  top  ice  crystal  layer  or  inhomogeneous 
structure.  The  surface  emission  from  various  underlying  media  including  calm 
and  rough  ocean,  sea  ice,  snowpack,  and  vegetation  are  taken  into  account. 
Functional  relationships  between  emitted  thermal  radiation  from  atmospheric 
precipitation  and  rainfall  parameters  are  discussed.  By  using  a  statistical 
regression  technique  (the  D-matrix  method),  the  geophysical  parameters,  rain 
rate,  and  emissivity  from  the  underlying  surface,  etc.,  are  retrieved  from 
simulated  data.  The  error  analysis  is  discussed. 

Radiative  transfer  theory  applicable  to  the  multiple  scattering  of  micro- 
waves  has  been  extensively  applied  to  the  study  of  atmospheric  precipitation 
and  clouds  (Wilheit  et  al,  1977;  Ishimura  and  Cheung,  1980;  Tsang  and  Kong, 
1977;  Jin  and  Kong,  1983).  The  previous  modeling,  however,  has  focused  on  the 
single  homogeneous,  isothermal  layer  of  Mie  scatterers.  In  wel 1  - devel oped 
precipitating  clouds,  it  has  been  observed  (Wilheit  et  al . ,  1982)  and  calcu¬ 
lated  by  the  one-layer  model  (Wilheit  et  al.,  1982;  Huang  and  Liou,  1983)  that 
an  ice  crystal  layer  over  precipitation  significantly  cools  the  emergent 
microwave  thermal  emission,  especially  at  higher  frequencies.  Here  a 
numerical  approach  is  described  to  solve  the  coupled  radiative  transfer 
equations  for  two-stratified,  nonisothermal  layers  composed  of  different 
(component  and  size  distributed)  spherical  Mie  scatterers.  It  is  applied  to 
calculations  of  emitted  thermal  radiation  from  atmospheric  precipitation  with 
a  top  ice  crystal  layer  or  with  nonhomogeneous  structure  (e.g.,  two  layers  of 
precipitation  with  different  rainfall  rates).  Two  radiative  transfer 
equations  are  solved  by  the  Gaussian  quadrature  method.  The  coefficients  of 
the  eigenvectors  in  the  two  respective  regions  are  determined  by  solving  four 
coupled  boundary  conditions.  The  emission  effects  from  various  underlying 
media,  including  calm  and  rough  sea  surfaces,  dry  and  wet  snowpacks ,  first- 
and  multi-year  sea  ice,  and  land  with  varying  soil  moisture  content  and 
vegetation  are  modeled  by  a  layer  of  scattering  medium  (Jin,  1984). 

Geophysical  parameter  retrieval  is  an  inverse  problem  for  remote  sensing 
of  atmospheric  and  earth  terrain  surface  properties.  The  retrieval  of 
atmospheric  temperature  and  water  vapor  profiles  have  been  studied  for  non¬ 
precipitating  atmospheres  based  upon  the  scalar  nonscattering  radiative 
transfer  equations  (Rodgers,  1976;  Gaut  et  al.,  1  973;  Rosenkranz  et  al  .  , 
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1982).  In  this  paper,  we  apply  the  statistical  D-matrix  retrieval  technique 
to  retrieve  the  geophysical  parameters  of  interest  in  remote  sensing  of  pre¬ 
cipitation  over  ocean  and  land.  By  a  Monte  Carlo  method,  the  most  probable 
combination  of  the  atmospheric  and  surface  properties  for  available  measure¬ 
ments  is  selected  in  a  statistical  sense.  Then,  by  our  simulation  modeling 
calculation  described  above,  a  set  of  simulated  radiometric  data  is  pro¬ 
duced.  By  using  a  statistical  regression  technique,  the  appropriate  D-matrix 
is  obtained  to  retrieve  the  desired  parameters,  such  as  rain  rate,  surface 
emissivities  for  both  vertical  and  horizontal  polarizations,  atmospheric 
opacity,  temperature  lapse  rate,  and  surface  temperature.  Results  indicate 
that  the  statistical  D-matrix  method  is  working  properly  for  the  retrieval  of 
these  geophysical  parameters  in  the  presence  of  atmospheric  precipitation. 
Some  restrictions  and  error  analyses  are  also  discussed. 

3.3.2  Thermal  Radiative  Transfer  Theory 

The  theoretical  model  is  shown  in  Figure  3-4.  Two  layers,  region 
f  and  region  1,  are  composed  of  different  spherical  Mie  scatterers  with  com¬ 
plex  dielectric  constants  and  respectively,  and  appropriate  size 

distributions.  Temperature  profiles  are  assumed  to  be  linear  with  the  lapse 
rate  and  T ^  (°c/km)  in  respective  regions.  From  the  vector  radiative 
transfer  theory  and  Gaussian  quadrature  method  (Tsang  and  Kong,  1977;  Jin  and 
Kong,  1983)  we  have  the  radiative  transfer  equations  as  follows: 


cos*i  d£  V*i-Z>  "  <kal  +  kgl>  CT1<Z>  -  <kel  +  kgl>V*i-Z> 

+J_N  ‘V^i-ViVV^  +  (^i*aj)iIcti(ej  -z)i- 


COS*if  dz  I^f(^if’z)  "  (kaf  1  kgf)  CTf(z)  *  (kef  +  kgf }  1/3f  (^if ,Z) 
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where  the  subscripts  1  and  f  denote  the  functions  or  parameters  taken  in  the 
respective  region  1  and  region  f;  0,  a  denote  vertical  or  horizontal  polari¬ 
zation  v  or  h;  ka,  ke  are  the  absorption  and  extinction  coefficients  of  the 
scatterers;  kg  is  the  atmospheric  gas  (oxygen  and  water  vapor)  absorption 
coefficient;  T(z)  is  the  physical  temperature  profile  of  atmosphere;  (0. ,  0- ) , 
(/?£,  a j ) ,  etc.,  are  the  scattering  phase  functions  calculated  by  Mie 
scattering  theory;  aj  are  Christofell  numbers,  N  is  the  number  taken  by 
Gaussian  quadrature.  The  temperature  profiles  are  assumed  as 

T1(z)  -  T^0)  -  (z+d)  •  r  (3.11a) 

Tf(z)  -  t£0)  -  z  •  Tf,  (3.11b) 

where  t(0) ,  t£0)  are  the  temperature  at  z  -  -d,  0  (Fig.  3-4).  In  the  case  of 
atmospheric  precipitation,  where  the  fractional  volume  of  scatterers  is  very 
small,  the  mean  permittivities  of  region  1  and  f  are  not  quite  different. 

Thus,  we  can  assume  that 


efi  -  Si  and  Nf  -  N 


(3.12) 


Therefore,  aj£  -  aj  .  Otherwise,  cost)  £^/cos0  ^  can  be  decided 
law,  and  a^£  can  be  calculated  by  the  trapezoidal  rule.  The 
conditions  are  obtained  as  follows: 
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where  the  underlying  surface  reflectivity  7^(0)  is  determined  by  the 
calculations  given  in  (Smith  et  al . ,  1979)  for  underlying  scattering  medium. 
From  the  Gaussian  quadrature,  we  obtain  the  solutions  of  equations  (3.10)  and 
(3.10)  in  terms  of  the  respective  eigenvectors,  corresponding  eigenvalues  and 
unknown  coefficients  which  will  be  determined  by  the  coupling  boundary  condi¬ 
tions  ((3.13a)  through  (3.13d).  Substituting  these  solutions  into  (3.13a)  - 
(3.13b)  we  obtain  the  desired  coefficients  by  matrix  inversion.  Then, 
finally,  we  obtain  the  observed  brightness  temperature  in  region  0  as 
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F.p  «  2N  x  2N  eigenvector  matrix 


Qp  ■  2N  x  2N  matrix  related  with  the  eigenvector 

diag  [e  ...,  e  Q2Nf^jj  (i  — 1 . 2N)  eigenvalue 


One  of  the  numerical  results  is  shown  in  Fig.  3-5.  We  take  N  =  b.  The  per¬ 
mittivity  of  esf,  <sj,  and  of  calm  ocean  are  calculated  by  the  Debye 
formula,  a  Marshall - Palmer  precipitation  dropsize  distribution  is  used  for 
liquid  and  solid  precipitation  in  both  regions  1  and  f.  Results  indicate  a 
significant  cooling  for  85.5  GHz  due  to  the  top  ice  layer,  even  at  moderate 
rain  rates.  Dotted  lines  show  the  case  without  the  top  ice  layer.  By  model¬ 
ing  the  underlying  media  as  a  layer  of  random  discrete  scatterers  and  choosim 
physically  reliable  parameters  for  the  size,  fractional  volume,  and  permit¬ 
tivity  of  the  scatterers  and  layer  thickness,  background  permittivity,  etc., 
the  emissivity  ep(6),  and  reflectivity  7^(0)  are  calculated  and  input  into  th< 
calculation  of  (3.14a)  and  (3.14b)  for  underlying  rough  sea  surface  with  foam 
FY  and  MY  sea  ice,  dry  and  wet  snow,  moistured  land  and  vegetation.  More 
results  and  detailed  discussions  will  be  given  in  the  future  paper. 

3.3.3  Statistical  Retrieval  of  Geophysical  Parameters 

By  the  Monte  Carlo  method,  we  choose,  in  a  statistical  sense,  most 
probable  combinations  of  geophysical  parameters  to  reproduce  the  available 
measurements  through  the  forward  calculations  by  our  theoretical  modeling. 

The  inverse  problem  is  whether  these  geophysical  parameters  can  be  inferred 
from  the  simulated  radiometric  measurements.  From  the  statistical  retrieva 
techniques  (Rogers,  1976;  Gaut  et.  al .  ,  1975;  Isaacs  and  Deblonde,  1987)  D- 
matrix  method,  the  estimated  parameters  of  p  is  generated  by 

p  -  D  «  4>  (3.1 

where  the  D-matrix  is  previously  computed  by 

D  -  C(p,*)  C'1  (£,*)  (3.1 


where  p,  d  are  the  parameter  and  data  vectors,  the  underbar  denotes  the 
vector,  <£  is  a  generalized  data  basis  function  and  dot  ined  as  <p  * 

col umn  (1,  (d . - d, )/d , ,  ....  (d  -d  ) /d  )  tor  m  channe !  r ad i ome  t i i c  da  t  a ,  vh 

m  m  m 

superbar  denotes  the  mean  value,  Cfp,^)  and  €<$,&)  ai<-  tin  'etiolation 
matrices  of  p  .  4>  and  ,  <fi ,  respectively. 


We  specify  the  most  interesting  parameters:  rain  rate,  emissivity  for 
both  vertical  and  horizontal  polarization  from  the  underlying  surface,  at¬ 
mospheric  opacity,  temperature  lapse  rate,  and  surface  temperature,  to  combine 
the  parameter  vector  p.  The  data  vector  d  is  obtained  from  the  simulated 
results  of  eight  channels  (four  frequencies,  19.35,  22.235,  37,  and  85.5  GHz) 
viewing  the  surface  at  a  nadiar  angle  of  59°.  By  using  100  observations  of 
ensemble  parameter  vector  pA  and  data  vector  dA  (set  A)  to  compute  D-matrix, 
another  ensemble  dg  of  set  B  is  then  used  to  estimate  the  parameter  p^,  and 
the  ground  truth  pg  will  be  compared  for  overall  error  analysis.  In  the 
presence  of  multiple  scattering  in  precipitation  and  surface  emission,  the 
relationship  between  d  and  p  is  much  more  sophisticated  than  a  simple  linear 
function.  However,  in  a  statistical  sense,  as  long  as  the  mean  value  and 
variance  of  the  ensemble  set  B  are  not  significantly  deviated  apart  from  a 
priori  set  A,  the  retrieval  techniques  of  D-matrix  is  working  fairly  well. 

Results  are  shown  in  Figs.  3-6  and  3-7  for  retrieval  of  rain  rate  and 
emissivity  of  underlying  surface.  It  is  shown  that  as  the  mean  value  and 
variance  of  set  B  is  close  to  those  in  a  priori  set  A,  retrievals  are  fairly 
good.  Large  deviation  from  a  priori  value  causes  large  error  in  retrieval 
since  D-matrix  relies  on  a  priori  information  as  given  in  (3.17). 
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The  geometry  of  the  theoretical  modeling. 
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Figure  3-7.  Retrieval  of  emissivities . 


3 . 4  Surface  Modeling 

The  Air  Force  Geophysics  Laboratory  (AFGL)  RADTRAN  code  is  used  to  cal¬ 
culate  transmission  and  brightness  temperature  spectra  for  atmospheric  slant 
paths  in  the  frequency  range  between  1  and  1000  GHz  (Falcone  et  al . ,  1979; 
Falcone  et  al.,  1982).  The  atmospheric  attenuation  submodels  of  the  clear 
atmosphere,  fog,  cloud,  and  rain  used  in  RADTRAN  have  been  thoroughly  docu¬ 
mented  (Falcone  et  al . ,  1982).  These  submodels  provide  the  necessary  data  to 
evaluate  gaseous  and  particulate  emission  and  other  potential  contributions  to 
slant  path  brightness  temperature  from  scattering  in  the  atmosphere  due  to  the 
presence  of  precipitation.  The  RADTRAN  code  is  frequently  applied  to  micro- 
wave  systems  effectiveness  simulation  studies  for  both  airborne  and  satellite- 
based  sensors  viewing  the  earth's  surface.  Such  applications  require  the 
capability  to  supplement  the  atmospheric  radiative  transfer  calculation  with 
realistic  models  of  the  microwave  properties  of  the  sensor  scene  f ield - of - vi ew 
allowing  for  variations  in  polarization  and  frequency  dependence  of  emissivity 


due  to  changes  in  surface  type.  Frequently,  this  information  is  provided  via 
available  empirical  data  sets.  However,  the  range  of  available  field  measure¬ 
ments  is  limited.  In  order  to  facilitate  the  evaluation  of  sensor  performance 
throughout  a  range  of  pertinent  geophysical  surface  scenarios,  deterministic 
models  of  microwave  surface  emission  are  desired  which  can  be  exercised  con¬ 
currently  with  the  models  for  atmospheric  transmission.  To  achieve  this  goal, 
a  set  of  microwave  surface  emission  models  has  been  implemented  to  accompany 
the  RADTRAN  atmospheric  properties  model  (Isaacs  et  al  .  ,  1988b).  In  addition 
to  obvious  applications  to  microwave  background  simulations,  application  of 
such  models  to  remote  sensing  simulation  studies  can  also  provide  a  theoreti¬ 
cal  basis  for  understanding  of  the  sensitivity  of  microwave  sensor  data  to 
variations  in  relevant  hydrological  parameters  such  as  soil  moisture,  snow, 
and  vegetative  moisture  content. 

3.4.1  Surface  Modeling  Approaches 


The  selection  of  a  simple  surface  modeling  approach  is  made  difficult  by 
the  complexity  of  geophysical  surfaces.  For  example,  homogeneous  dielectric 
slab  models  for  both  land  and  ocean  have  been  commonly  used  to  provide  the 
required  surface  emissivity  parameters  to  initiate  t rightness  temperature 
simulation  calculations.  It  Is  apparent  from  an  examination  of  recent  micro- 
wave  satellite  imagery  that  such  approaches  cannot  reproduce  the  fidelity  of 
the  complex  fields  of  observed  surface  properties.  This  is  due  to  the  neglect 
of  important  physical  mechanisms  such  as  scattering  by  such  approaches,  and 
their  failure  to  treat  the  spatial  inhomogeneities  in  dielectric  properties 
due  to  the  inherent  physical  structures  of  real  surfaces.  On  the  other  hand, 
it  is  necessary  to  consider  the  computational  level  of  effort  required  to 
model  all  of  the  physics,  even  if  it  were  well  understood. 

Related  to  the  choice  of  an  appropriate  model,  is  the  detail  of  surface 
type  characterization  desired,  i.e.  how  many  surface  types  to  treat.  In 
reality,  of  course,  there  is  a  continuum  of  geophysical  surface  types.  A 
sufficiently  general  model  can  attempt  to  simulate  some  of  the  behavior 
exhibited  by  subsets  of  surface  types  within  this  continuum  by  choosing  ap¬ 
propriate  parameter izations  of  relevant  surface  properties  and  varying  them 
within  representative  ranges.  The  choice  of  surface  types  employed  within  the 
RADTRAN  surface  modeling  package  was  based  both  on  the  desire,  to  treat  a  com- 


prehensive  set  of  surfaces  and,  to  some  extent,  on  the  requirements  of  poten¬ 
tial  model  users  with  specific  surface  related  simulation  applications.  The 
surface  types  selected  are:  (a)  calm  and  rough  ocean,  (b)  first  year  (FY)  and 
multiyear  (MY)  sea  ice,  (c)  wet  and  dry  snow  over  land,  (d)  wet  and  dry  soil, 
(e)  vegetation,  and  (f)  land.  The  land  surface  type  provides  a  background  for 
snow,  soil,  and  vegetation  models  in  addition  to  its  potential  role  as  a 
distinct  surface  type  itself. 

Both  the  calm  ocean  and  land  are  modeled  as  simple  dielectric  slabs.  The 
other  surface  types,  however,  clearly  require  a  more  sophisticated  modeling 
treatment.  As  the  following  discussion  will  indicate,  it  is  not  appropriate 
to  treat  all  of  the  surface  types  delineated  above  by  a  single  formalism. 
Therefore,  two  distinct  approaches  have  been  applied  in  the  development  of 
these  surface  emission  models:  that  based  on  wave  theory  for  random  discrete 
scatterers  and  that  based  on  radiative  transfer  theory  for  continuous  random 
media.  In  the  sections  to  follow,  the  former  approach  is  applied  to  modeling 
the  ocean  surface,  sea  ice,  and  snow,  while  the  latter  is  applied  to  both  soil 
and  vegetation.  These  approaches  are  summarized  in  Table  3-6. 

Table  3-6.  Surface  Model  Types  and  Modeling  Approaches 


Model 

Surface  Type 

Modeling  Approach 

1 

Calm  ocean 

Dielectric  slab 

2 

Rough  ocean 

Random  discrete 

scatterers 

3 

FY  sea  ice 

Random  discrete 

scatterers 

4 

MY  sea  ice 

Random  discrete 

scatterers 

5 

Dry  snow 

Random  discrete 

scatterers 

6 

Wet  snow 

Random  discrete 

scatterers 

7 

Land 

Dielectic  slab 

8 

Wet  soil 

Continuous  random  medium 

9 

Vegetation 

Continuous  random  medium 

In  the  wave  theory  approach  for  random  discrete  scatterers,  one  or  more 
layers  is  defined  consisting  of  a  dielectric  medium  with  either  taniform 
properties  or  containing  a  random  distribution  of  discrete  dielectric  spheres 
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with  distinct  dielectric  properties.  These  latter  inclusions  give  the  medium 
scattering  properties  which  by  appropriate  choice  of  the  background  and  in¬ 
clusion  permittivities  can  be  tuned  to  exhibit  the  observed  behavior  of  sea 
ice  and  dry  snow,  for  example.  The  radiative  transfer  approach  for  continuous 
random  media  approaches  the  problem  of  modeling  the  surface  from  a  different 
perspective.  Some  surfaces  are  spatially  inhomogeneous  in  their  dielectric 
properties,  yet  the  inhomogenei ties  are  not  due  to  discrete  spherical 
scatterers.  The  approach  provides  an  alternative  treatment  in  which  the 
permittivity  is  varied  continuously  throughout  the  medium.  Furthermore,  these 
spatial  variations  are  parameterized  in  such  a  manner  that  the  relative 
effects  of  variations  in  the  vertical  and  horizontal  physical  structure  of  the 
medium  can  be  modeled.  The  application  of  this  approach  to  soils  and 
vegetation  is  obvious . 


3.4.2  Application  of  Wave  Theory  Models  for  Discrete  Scatterers 


A  variety  of  surface  types  can  be  modeled  as  either:  (a)  a  single  layer 
with  uniform  dielectric  properties,  (b)  a  single  layer  including  a  random 
distribution  of  discrete  scatterers,  or  (c)  two  bounded  layers,  each  with  one 
of  the  above  characteristics.  The  most  straightforward  example  of  the  first 
type  is  the  calm  ocean  surface  which  has  been  modeled  in  this  study  as  a 
homogeneous  half  space  medium  with  microwave  dielectric  properties  calculated 
by  the  Debye  formula.  Figure  3-8  illustrates  the  resultant  variation  of 
emissivity  with  frequency  for  each  polarization  and  a  look  angle,  of  54°.  This 
angle  was  chosen  to  correspond  to  that  of  a  conically  scanning  radiometer  such 
as  the  NIMBUS  7  Scanning  Multichannel  Microwave  Radiometer  (SMMR)  or  DMSP 
Special  Sensor  Microwave  Imager  (SSM/I).  The  emissivity  of  a  rough  sea 
surface  is  greater  than  that  of  a  smooth  sea  because  the  wind-driven  rough  sea 
surface  is  covered  with  a  layer  of  foam  and  white  caps.  Calm  seas  are  defined 
as  those  for  which  the  ocean  surface  wind  speed  is  less  than  about  7  ms’^, 
while  rough  seas  are  present  when  surface  wind  speeds  exceed  this  critical 
value.  The  effect  of  white  water  (foam  and  white  caps)  have  been  modeled  as  a 
layer  of  porous  dielectric  (Droppleraan,  1970)  and  a  series  of  thin  water  films 
(Rosenkranz,  1972).  In  this  study,  the  wind-driven  rough  sea  is  modeled  by  a 
layer  of  random  discrete  scatterers  with  a  periodic  rough  surface.  The 
dependence  of  emissivity  on  look  angle  is  shown  in  Figure  3-9.  Numerical 
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results  are  consistent  with  observations  (Rosenkranz,  1972).  Increased  emis- 
sivity  due  to  the  foam  layer  is  about  0.09  at  a  look  angle  of  11°,  and  the 
resultant  brightness  temperature  increase  is  about  11°,  which  corresponds  to  a 
wind  speed  of  about  10  m/sec.  Increase  of  the  fractional  volume  of  foam  would 
also  increase  the  brightness  temperature. 

Snowpack  over  land  is  modeled  as  a  bounded  layer  of  random  discrete 
scatterers.  Dry  snow  is  composed  of  ice  particles  and  air,  while  wet  snow 
consists  of  small  water  drops,  ice  particles,  and  air.  The  wave  theory  ap¬ 
proach  for  a  bounded  layer  of  random  discrete  scatterers  is  used  to  calculate 
the  emissivity  (Jin  and  Kong,  1984).  This  approach  includes  first  and  second 
order  scattering  effects  and  is  especially  applicable  to  modeling  densely 
distributed  random  discrete  scatterers.  Figure  3-10  shows  the  frequency 
dependence  of  snowpack  emissivity  at  a  look  angle  of  54°.  A  snow  depth  of  60 
cm  is  assumed.  As  frequency  increases,  the  effect  of  scattering  is 
enhanced.  This  results  in  a  decrease  of  emissivity.  Figure  3-11  also  shows 
the  variation  in  snowpack  emissivity  with  wetness.  As  wetness  increases,  the 
emissivity  of  the  snowpack  gradually  increases  due  to  increased  thermal  emis¬ 
sion.  Snowpack  emissivity  as  a  function  of  angle  at  18.0  OH z  is  illustrated 
in  Figure  3-11.  For  comparison,  the  angular  variation  of  emissivity  for  a 
snow-free  land  surface  is  included.  There  is  a  significant  decrease  in  the 
vertical  emissivity  due  to  the  presence  of  snow.  This  difference  is  modulated 
by  variations  in  both  snow  depth  and  wetness. 

The  approach  can  also  be  applied  to  calculate  the  emissivities  of  first 

year  (FY)  and  multiyear  (MY)  sea  ice.  It  is  known  that  there  is  a  strong 

contrast  between  the  emissivity  of  sea  ice  and  that  of  the  ocean.  The  erai  s - 

sivity  of  sea  ice  depends  on  age,  thickness,  salinity,  density,  composition 

and  surface  characteristics  (Comiso,  1983).  To  some  extent,  thickness,  t,  and 

density,  p,  are  correlated  with  ice  age.  First  year  sea  ice  has  a  layer  of 

3 

thin  congelation  ice  (  20  cm,  p  ~  0.85  gm/cm  )  over  a  thick  congelation  ice 

3 

layer  (p  =  0.92  gm/cm  ) .  Multiyear  sea  ice  has  a  recrystallized  ice  layer 

(t  -  50  cm)  over  a  congelation  ice  layer  with  thickness  up  to  500  cm.  The 

•> 

recrystallized  ice  layer  of  MY  sea  ice  (p  ~  0.7  gm/cm  )  has  a  large  fractional 
volume  of  air  pockets  which  cause  scattering  and  darken  thermal  emission  from 
the  underlying  medium.  In  our  calculations,  MY  sea  i  <_  f  was  modeled  as  a 
bounded  layer  of  random  discrete  scatterers  (air  pockets  embedded  in  homop.ene- 


ous  ice  medium)  over  a  homogeneous  ice  medium,  while  FY  sea  ice  is  modeled  as 
two  layers  of  homogeneous  ice  media  with  appropriate  complex  dielectric 
constants.  Calculated  emissivities  as  functions  of  frequency  and  look  angle 
for  FY  and  MY  sea  ice  are  shown  in  Figures  3-12  and  3-13,  respectively. 
Results  are  fairly  consistent  with  the  observations  and  experimental  data 
(Comiso,  1983;  Simonett,  1983). 

3.4.3  Application  of  Radiative  Transfer  Theory  Models  for  Continuous 


Random  Media 


Experimental  measurements  of  the  thermal  emission  from  wet  soil  and 
vegetative  canopies  have  been  studied  by  a  number  of  investigators  (Ulaby  et 
al .  ,  1983;  Ulaby  and  Jedlicka,  1984,  Schmugge  et  al.,  1977,  Schmugge,  1983)  by 
means  of  truck- mounted ,  airborne,  and  satellite  radiometers  at  low  frequen¬ 
cies,  such  as  1.4  and  5  GHz.  Measurements  at  higher  frequencies  have  not  been 
well  studied  due  to  some  difficulties  (Ulaby  et  al . ,  1983).  Theoretical  anal¬ 
yses  are  usually  restricted  to  finding  an  effective  permittivity  for  wet  soil 
and  vegetation  at  these  frequencies.  Modeling  using  the  radiative  transfer 
equation  including  multiple  scattering,  has,  however,  also  been  investigated 
(Wang  and  Schmugge,  1980;  Wang  et  al . ,  1984). 

In  the  RADTRAN  surface  emissivity  model  implementation  wet  soil  and 
vegetation  are  modeled  as  layers  of  continuous  random  medium  bounded  by  an 
underlying  homogeneous  medium.  Random  dielectric  fluctuations  in  these  media 
are  characterized  by  their  respective  spatial  correlation  functions.  The  mean 
permittivity  and  variance  are  calculated  by  using  the  appropriate  multiphase 
mixture  formula  (Ulaby  and  Jedlicka,  1984;  Wang  and  Schmugge,  1980).  A 
numerical  Gaussian  quadrature  method  similar  to  that  described  in  (Jin  and 
Isaacs,  1987)  is  then  used  to  solve  the  radiative  transfer  equation  to 
calculate  emissivity. 

The  theoretical  model  for  wet  soil  or  vegetative  canopy  consists  of  a 
layer  of  random  medium  with  permittivity  <(r),  with  mean  value,  t  and 
fluctuating  part,  «jp(r): 


<(r)  -  <  +  «f(r) , 

ra  t 

«f(r)>  -  0. 


(3.18) 
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The  mean  value  of  the  permittivity,  <  may  be  related  to  the  relevant  mixture 
components  and  the  physicaL  structure  of  the  layer.  In  our  soil  model,  the 
empirical  mixture  formula  given  by  Wang  and  Schmugge  (Wang  and  Schmugge ,  1980) 
is  used: 


For  Wc  <  Wt 


£m  -  Vx  +  fp  *  Wc)  £b  +  ^  -  P>  f rock 


(3. 19a) 
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For  Wc  >  Wt 
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where  Wc  is  the  volumetric  moisture,  and  Wc  is  the  transient  moisture.  The 
transient  moisture  denotes  the  turning  point  for  transition  of  the  dielectric 
permittivity  from  slow  to  steep  behavior,  £x  stands  for  the  dielectric 
constant  of  the  initially  absorbed  water,  7  is  a  parameter  (about  0.3  to  0.5), 
£ice*  'water'  £rock'  anc*  £b  arc  dielectric  constants  for  ice,  water,  rock, 

and  background  medium,  respectively  Equation  3.19  is  an  empirical  formula 
which  is  based  on  the  behavior  of  the  biphase  dielectric  properties  of  water 
in  soils,  i.e.,  below  the  transient  moisture,  the  dielectric  property  of  the 
soil  water  is  thought  to  be  like  that  of  ice,  while  above  W  ,  the  water 
dielectric  constant  is  used.  In  our  calculations  the  following  constants  are 
used:  Wt  -  0.25,  7  -  0.95,  P  -  0.5,  erock  -  (5,  0),  «ice  -  (3,  0.01),  and 
£water  Ps  calculate<^  by  the  Debye  formula  which  depends  on  the  temperature  and 
frequency  (Sadiku,  1985)  . 

For  a  vegetative  canopy,  the  four-phase  refractive  model  of  the  mixture 
following  (Schmugge  et  al . ,  1977)  is  used  as 
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where  V  ,  Vgw,  and  V^w  are  the  volumetric  fraction  of  air,  free  water,  and 
bound  water  in  the  vegetation,  respectively;  nt^  is  the  volumetric  wetness  of 
the  vegetation,  is  taken  as  that,  of  ice,  and  the  bulk  dielectric  constant, 

£bulk  ta^en  as  (5.5,  0.1).  Therefore,  the  mean  value ,  t  in  the  layer  of 
vegetation  is 
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where  Vv  is  the  volumetric  layer  fraction,  taken  to  be  0.0033  and  the  quantity 
(1  -  nty  -  Vz)  =  0.1  (Ulaby  and  Jedlicka,  1984).  The  variance  6  is  calculated 
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The  above  formulas  (3.19)  and  (3.20)  and  parameters  are  consistent  with  avail¬ 
able  measurements  of  the  permittivity  at  1.4  and  5  GHz  (Ulaby  and  Jedlicka, 
1984;  Wang  and  Schmugge ,  198C) .  These  expressions  were  extended  to  higher 
frequencies  (19.35,  22.235,  and  37  GHz)  for  calculation  of  the  mean  value  «  , 
based  on  the  observation  that  the  dielectric  constants  of  ice  and  rock  are  not 
very  sensitive  to  changes  in  the  frequency  in  this  region,  while  the  dielec¬ 
tric  constant  of  water  can  be  calculated  by  the  Debye  formula.  It  has  been 
assumed  that  the  possible  anisotropic  and  inhomogeneous  properties  of  vegeta¬ 
tion  are  not  significant  enough  to  effect  this  approximation. 

The  correlation  function  of  dielectic  fluctuation  in  random  media  is 
usually  assumed  to  be  Gaussian  in  the  transverse  direction  and  exponential  in 
the  vertical  (Tsang  and  Kong,  1976). 
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The  transverse  and  vertical  correlation  lengths  i  ,  i  ,  and  variance,  6,  which 
characterize  the  magnitude  of  scattering  effects,  are  dependent  on  the  size, 
shape,  and  spatial  distribution  of  the  scatterers.  An  analysis  applying  these 
concepts  to  simple  dry  snowpack  can  be  found  in  (Vallese  and  Kong,  1981).  For 
the  more  complicated  case  of  a  vegetative  canopy,  it  is  possible  to  define  the 
optimal  correlation  lengths  which  give  the  proper  scattering  intensity 
magnitudes . 

By  using  the  radiative  transfer  equation  in  region  1,  and  applying  the 
Gaussian  quadrature  method  (Jin  and  Isaacs,  1987;  Tsang  and  Kong,  1976), 
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where  0  and  a  denote  the  vertical  and  homogeneous  polarization  v  or  h,  and 
and  kg  are  the  absorption  and  extinction  coefficients,  respectively,  in 
region  1.  Sixteen  quadrature  points  are  used  in  the  calculation.  The 
scattering  phase  functions  are  as  follows  (Tsang  and  Kong,  1976). 
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The  boundary  conditions  are: 
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Therefore,  the  brightness  temperatures  at  the  altitude  Z  -  D  are  given  by: 
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where  8 £0  are  determined  by  Snell's  law 
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(3.28) 


Eqn.  (3.27)  provides  the  required  brightness  temperature  solutions  above  an 
atmosphere  with  gas  absorption  coefficient  equal  to  kg(z) . 

The  roughness  of  the  underlying  medium,  region  2,  is  treated  following 
(Choudhury  et  al. ,  1979)  where  only  the  coherent  part  of  the  scattered  field 
is  taken  into  account 


(o)., 1\  -hcos20 
T/9i2(*7  “  T/912 


(3.29) 


where  is  c*ie  re^ lectivity  for  a  smooth  surface,  and  h  is  the  effective 

roughness  which  is  generally  smaller  than  the  physical  roughness.  Tsang  and 
Newton  (Tsang  and  Newton,  1982)  obtained  a  more  sophisticated  formula  which 
includes  the  incoherent  part  and  calculated  relevant  results  for  1.4  GHz  at  a 
nadir  angle  for  a  half- space  model.  In  the  two- layer  model  with  a  rough 
boundary  at  z  -  -d,  the  incidence  of  thermal  emission  should  be  accounted  for 
from  all  angles;  the  approximated  calculation  near  the  nadir  angle  cannot  be 
used.  Therefore,  our  model  uses  only  the  effective  h  to  take  account  of  the 
roughness  contribution  as  given  in  (3.29).  In  fact,  for  fully-grown  vege¬ 
tation,  the  vegetative  canopy  will  likely  shield  the  soil  contribution 
completely,  and  the  contribution  of  (3.29)  from  region  2  is  probably  not 
significant  at  all. 

Frequency  dependent  vertical  and  horizontal  emissivities  for  dry  and  wet 
soils  are  illustrated  in  Figure  3-14.  Corresponding  look  angle  dependent 
emissivity  results  evaluated  at  a  frequency  of  18  GHz  are  shown  in 
Figure  3-13.  The  model  used  consists  of  tvo  layers:  an  upper  layer  with  a 


j"  V*  V"  "V  V  V  **>  V 


; -  'vvv:  ^ 


depth  of  50  cm  and  variable  soil  moisture  and  a  lower  layer  with  a  fixed 
wetness  of  0.2.  The  dry  soil  has  a  moisture  volume  of  0.1  while  the  wet  soil 
has  a  moisture  volume  of  0.5.  Calculations  are  shown  for  correlation  lengths 

of  (2  2  )  equal  to  (0.1,  1.0),  respectively.  The  observed  decrease  of 

P  z 

brightness  temperature  with  increased  soil  moisture  for  both  polarizations  is 
illustrated.  Notably,  the  magnitude  of  emissivity  is  much  more  dependent  on 
the  specified  correlation  lengths  at  low  values  of  soil  moisture  where 
scattering  is  important. 

Applying  the  model  to  vegetation,  specification  of  the  appropriate  cor¬ 
relation  lengths  becomes  significant.  The  correlation  lengths  are  defined 
such  that  as  the  distance  between  two  points  becomes  larger  than  the  correla¬ 
tion  length,  their  correlation  will  be  less  than  e'1.  Therefore,  the  cor¬ 
relation  lengths  actually  correspond  to  a  measure  of  the  sc atterer  size, 
shape,  and  spatial  distribution.  Increasing  2  or  2  means  that  the  scat- 
terers  are  distributed  more  horizontally  or  vertically,  respectively.  There¬ 
fore,  the  appropriate  2  or  2  for  different  vegetation  canopies  are  important 
at  high  frequencies.  If  scatterers  are  spherical,  the  correlation  length  is 
approximately  equal  to  the  radius.  For  fully  grown  corn,  for  example,  2 z  is 
larger  than  2p  due  to  the  predominantly  vertical  physical  structure.  If  the 
random  medium  has  a  laminar  structure,  i.e.  one  where  vegetation  is  oriented 
predominantly  in  the  horizontal  direction  then  2  -*<*>.  For  this  geometry,  the 

behavior  of  emissivity  simulates  that  of  a  stratified  medium. 

Figures  3-16  and  3-17  illustrate  emissivity  results  for  four  different 
vegetation  models  based  on  two  vegetation  moisture  volumes  (01,  0.5)  and  two 
thicknesses  (20  and  200  cm,  respectively) .  Calculations  are  shown  for  cor¬ 
relation  lengths  of  (2  2  )  equal  to  (0.1,  1  0),  respectively.  In  general, 

the  behavior  of  the  emissivity  results  can  be  understood  in  terms  of  the  in¬ 
fluence  of  the  scattering  mechanism.  Scattering  is  enhanced  for  higher 
vegetative  moisture  volumes  and  thicknesses  and  tends  to  depolarize  the  sur¬ 
face  emission.  At  higher  frequencies,  this  depolarization  due  to  vegetation 
becomes  more  significant.  Thus  referring  to  Figures  3-16  and  3-17,  vegetation 
model  3  (thickness,  20u  cm;  moisture  volume,  0.5}  exhibits  the  lowest  surface 
emission  polarization  across  this  frequency  range,  while  model  2  (thickness, 

20  cm;  moisture  volume,  0.1)  exhibits  the  highest.  These  vegetation  models 
are  characterized  by  high  and  low  scattering,  respectively,  with  models  1  and 
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4  exhibiting  intermediate  scattering  properties  As  can  be  seen  in 
Figure  3-17,  complete  depolarization  is  reached  at  progressively  lower  fre¬ 
quencies  as  the  scattering  character  of  the  vegetation  increases  as  defined 
above.  In  fact,  model  3  is  unpolarized  at  all  frequencies  within  this  domain, 
while  model  2  is  totally  polarized. 

Using  the  vegetation  surface  emissivity  models,  brightness  temperatures 
were  evaluated  for  comparison  to  available  field  data.  All  ground  truth 
measurements  are  given  by  Wang  et  al.  (Wang  et  al.,  1984).  Figures  3-18  and 
3-19  compare  model  results  with  measurements  at  1.4  GHz  for  fully-grown  corn 
and  cut  corn;  Figure  3-21  shows  the  results  at  5  GHz  for  dry  corn.  In 
general,  the  model  results  at  these  frequencies  are  quite  good.  This  provides 
confidence  to  extend  the  model  to  higher  frequencies.  Calculations  also 
indicate  that  scattering  effects  are  negligible,  especially  at  1.4  GHz,  since 
the  brightness  temperatures  Tg  are  not  sensitive  to  any  change  of  £  ,  £z,  or 
letting  5-0. 

At  higher  frequencies,  furthermore,  the  mean  values,  and  de¬ 
crease.  Therefore  the  emissivities  of  the  vegetative  canopy  become  larger  and 
resultant  brightness  temperatures  Tfi  will  be  higher  than  those  at  low  frequen¬ 
cies.  As  the  observation  angle  increases,  the  random  layer  of  scatterers 
exhibits  more  backscattering,  which  in  turn  darkens  thermal  emission.  There¬ 
fore,  at  high  frequencies,  brightness  temperature  decreases  monotonically  as 
the  observation  angle  increases.  This  behavior  for  higher  frequencies  is 
different  from  that  at  low  frequencies  (cf.  Figures  3-18  through  3-19).  At 
larger  angles,  warming  effects  due  to  atmospheric  gas  absorption  in  the 
atmospheric  layer  above  the  surface  also  become  significant.  The  effect  of 
multiple  scattering  also  increases  with  increasing  frequency.  This  effort 
significantly  darkens  the  thermal  emission.  The  appropriate  choice  of  cor¬ 
relation  length  and  their  relative  magnitudes  also  affects  resultant  bright¬ 
ness  temperatures.  As  i  »  £  ,  the  laminar  structure  will  cause  more  back- 

P  ^ 

scattering  and  significantly  cool  the  thermal  emission.  For  a  fixed  volu¬ 
metric  moisture  fraction,  increases  in  the  magnitudes  of  £pl  £z  will  also 
enhance  backscattering  and  cool  (i.e.,  decrease)  the  thermal  emission.  The 
larger  the  ratio  i/£  the  warmer  the  brightness  temperature.  Conversely, 

z  p 

decreasing  the  correlation  length  gives  warmer  brightness  temperatures,  since 
there  is  less  scattering.  The  modeled  brightness  temperature  is  also  a  func- 


tion  of  the  thickness  and  wetness  of  the  vegetation.  As  the  vegetation  wet¬ 
ness  increases,  the  brightness  temperature  will  increase  due  to  more  thermal 
emission.  For  additional  increases  in  wetness,  saturation  is  reached.  This 
characteristic  behavior  is  enhanced  for  thicker  vegetation.  Once  saturation 
is  reached,  brightness  temperatures  decrease  since  more  scattering  and 
attenuation  cool  the  emission.  This  behavior  is  different  from  that  at  low 
frequencies,  where  brightness  temperature  monotonically  increases. 
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3.4.4  Comparison  with  SMMR  Data 

3 . 4 . 4 . 1  Data  Sample 

Nimbus  7  Scanning  Multichannel  Microwave  Radiometer  (SMMR)  data  was  used 
to  provide  a  comparison  to  RADTRAN  simulated  sensor  data  for  a  selected  subset 
of  relevant  surface  models.  The  choice  of  surface  types  simulated  was  deter¬ 
mined  by  the  availability  of  suitable  SMMR  data  sets.  The  SMMR  observes  ten 
channels  (five  dual-polarized  frequencies:  6.6,  10.7,  18,  21,  and  37  GHz)  and 
scans  conically  with  an  angle  of  incidence  on  the  surface  of  the  earth  of  50 
degrees.  The  SMMR  is  described  in  Gloersen  and  Barath  (1977)  and  NASA 
(1978).  For  selected  Nimbus- 7  passes,  therefore,  SMMR  data  were  acquired  from 
NSSDC  (National  Space  Science  Data  Center)  in  the  TCT  (Calibrated  Temperature 
Tape)  format.  The  TCT  format  presents  the  data  as  brightness  temperatures  for 
actual  SMMR  footprints,  with  the  FOV  locations  determined  from  the  pierce 
points  corresponding  to  the  four  antenna  angles .  The  locational  accuracy  is 
good,  usually  within  half  a  footprint  (A.  Chang,  personal  communication). 

SMMR  data  at  frequencies  of  6.6,  10.6,  21,  and  37  GHz  were  sampled  for 
selected  test  areas  from  brightness  temperature  maps  created  from  the  SMMR 
data  tapes.  Both  vertical  and  horizontal  polarizations  were  used.  Test  areas 
were  chosen  to  provide  data  on  specific  earth  surface  types  including  the  calm 
ocean  (Case  1),  moderate  vegetation  (Case  2),  and  heavy  vegetation  (Case  3). 
All  brightness  temperatures  were  spatially  smoothed  to  the  resolution  of  the 
6.6  GHz  channels.  Spatial  resolution  at  this  frequency  is  approximately 
160  km.  It  was  not  necessary  to  insure  that  the  observed  field-of  view  was 
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spatially  uniform  since  it  is  our  purpose  to  provide  simulation  models  which 
are  consistent  with  data  averaged  over  typical  sensor  footprints. 

3.4.4. 2  Results  of  Data  ComDarison 


Simulated  SMMR  brightness  temperature  data  were  calculated  using  the 
RADTRAN  atmospheric  transmission/brightness  temperature  model.  Appropriate 
atmospheric  profile  input  data  for  temperature  and  water  vapor  were  obtained 
from  a  summer  midlatitude  model.  This  model  also  provided  surface  temperature 
data.  SMMR  data  collected  over  the  north  central  United  States  provides  a 
surface  background  of  moderate  vegetation.  Three  subregions  were  selected  and 
frequency  dependent  SMMR  brightness  temperature  data  were  plotted.  These  data 
for  each  area  are  illustrated  in  Figure  3-21.  Both  vertical  and  horizontal 
polarizations  are  shown  for  each  area  and  frequency.  The  data  for  area  1 
shows  a  vertical  brightness  temperature  near  253  K  which  is  near  constant  with 
frequency.  The  horizontal  brightness  temperature  for  this  case  is  about  25  K 
colder  with  the  polarization  difference  decreasing  with  frequency.  The  data 
from  areas  2  and  3  are  qualitatively  similar  except  that  brightness  tempera¬ 
tures  are  about  18  K  higher.  The  corresponding  model  simulations  are  also 
presented  in  Figure  3-22  for  comparison.  The  vegetation  model  was  used  to 
evaluate  surface  emissivity.  Four  vegetation  types  are  shown  corresponding  to 
variations  in  the  thickness  of  the  vegetation  and  its  moisture  content.  T^e 
least  amount  of  scattering  is  given  by  the  vegetation  2  model.  This  can  be 
diagnosed  from  the  large  polarization  differences  at  low  frequencies  which 
decreases  with  frequency.  This  result  for  a  vegetation  depth  of  20  cm  and 
vegetative  moisture  content  of  0.1  is  very  similar  to  that  for  bare  moist 
soil,  e.g.  the  vegetation  does  little  to  mask  the  soil  emission.  The  most 
scattering  case  is  for  vegetation  3  with  a  thickness  of  200  cm  and  a  moisture 
value  of  0.5.  In  this  model  the  polarized  surface  is  completely  masked  by 
scattering  in  the  vegetative  layer  and  the  polarization  difference  is  very 
small.  An  examination  of  Figure  3-22  indicates  that  the  domain  of  the  model 
simulations  brackets  the  actual  SMMR  data.  That  is  on  average  the  properties 
of  the  vegetation  in  the  SMMR  fields-of-view  plotted  in  Figure  3-21  are 
between  those  of  light  and  heavy  vegetation.  No  attempt  was  made  to  tune  the 
vegetation  emission  models  to  fit  the  data,  however,  it  is  clear  that  some 
intermediate  degree  of  scattering  would  suffice.  This  could  be  achieved  by 
adjusting  the  thickness  or  moisture  content. 
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The  second  set  of  comparisons  uses  SMMR  data  from  the  southeastern  United 
States.  Again  three  subareas  were  chosen.  Areas  1  and  2  were  selected  over 
the  calm  ocean.  Calmness  was  judged  from  the  magnitude  of  37  GHz  brightness 
temperature  and  precipitating  areas  were  avoided  by  examining  the  37  GHz 
polarization  difference.  The  third  area  was  chosen  over  a  heavily  vegetated 
region  of  central  Florida.  SMMR  brightness  temperatures  for  these  areas  are 
plotted  as  a  function  of  frequency  in  Figure  3-22.  It  can  be  seen  that  the 
signatures  of  the  first  two  areas  and  that  of  area  three  are  quite  distinct. 
Furthermore  areas  1  and  2  are  very  similar.  The  simulated  brightness  tempera¬ 
tures  for  the  oceanic  cases  are  also  shown  in  Figure  3-22.  Here,  both  calm 
and  rough  sea  models  are  illustrated.  Comparison  of  these  results  with  the 
data  for  areas  1  and  2  shows  excellent  agreement  for  the  calm  ocean  simula¬ 
tion.  Indeed,  surface  wind  speeds  for  this  time  and  area  are  reported  to  be 
low  (i.e.  less  than  3  m/s).  Simulation  results  for  the  area  3  are  shown  in 
the  upper  part  of  Figure  3-22.  Again  four  possible  models  of  the  surface 
vegetation  are  illustrated  which  bracket  the  data.  The  small  polarization 
difference  of  the  data  is  consistent  with  the  heavy  vegetation  model  with  high 
scattering,  however,  the  magnitudes  of  the  brightness  temperatures  for  data 
and  model  are  offset.  This  is  most  likely  due  to  a  difference  between  the 
actual  and  modeled  surface  temperatures. 


Summary  and  Conclusions 


A  surface  property  modeling  subpackage  has  been  described  which  supple¬ 
ments  the  atmospheric  transmission  and  radiance  simulation  capability  provided 
by  the  AFGL  RADTRAN  computer  algorithm.  A  variety  of  geophysical  surface 
types  are  treated  by  the  subpackage  which  calculated  frequency  and  polariza¬ 
tion  dependent  surface  emissivities  for  use  by  RADTRAN  in  support  of  earth 
viewing  microwave  remote  sensing  system  simulations.  The  theoretical  basis 
for  the  emissivity  calculation  depends  on  surface  type.  Two  distinct 
approaches  are  implemented:  that  based  on  wave  theory  for  discrete  scatterers 
and  that  based  on  radiative  transfer  theory  for  continuous  random  media.  The 
former  approach  is  used  to  moJel  the  ocear.  curf a'"' ,  v-ai  Lous  forms  of  sea  ice , 
and  snow  over  land,  while  the  latter  is  used  to  treat  soils  and  vegetation. 

Calculations  of  polarized  surface  emissivity  using  these  models 
illustrate  the  diverse  frequency  dependent  signatures  of  the  various  surface 


types.  The  angular  dependence  of  surface  emissivity  is  also  illustrated. 

Using  these  surface  emissivities  as  input  to  the  RADTRAN  atmospheric  trans¬ 
mission  submodels,  simulated  brightness  temperature  spectra  have  been  computed 
for  a  variety  of  surface  types.  These  simulations  have  been  compared  with 
Nimbus-7  SMMR  data  collected  over  selected  areas.  The  comparisons  illustrate 
that  the  simulation  algorithm  brightness  temperature  values  based  on  the 
calculated  surface  emissivities  provide  a  reasonable  qualitative  description 
of  the  spectral  properties  of  the  SMMR  brightness  temperatures.  Further 
comparisons  are  necessary  to  validate  the  models  especially  for  moist  soil, 
snow,  and  sea  ice. 

In  summary,  the  emissivity  modeling  subpackage  significantly  enhances  the 
capability  of  the  RADTRAN  code  to  provide  realistic  brightness  temperature 
simulations  for  sensors  viewing  geophysical  surfaces. 
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Figure  3-8. 


Frequency  dependence  of  ocean  surface  emissivity  for  calm 
(u  <  7ms"^)  and  rough  (u  >  7ms’^)  ocean  at  a  look  angle  of  54' 


Figure 


-13.  Look  angle  dependence  of  first  year  and  multiyear  sea  ice 
emissivity  at  a  frequency  of  18.0  GHz  (Ts  -  273. K). 
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Figure  3-14. 


Frequency  dependence  of  dry  and  wet  soil  emissivity  at  a 
angle  of  54 0 . 


4.  SSM/I  FOOTPRINT  CHARACTERIZATION  USING  OTHER  DMSP  SENSORS 
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4 . 1  Background 

Observed  brightness  temperatures  at  window  frequencies  used  in  the 
retrieval  of  temperature  and  water  vapor  near  the  earth's  surface  are  quite 
sensitive  to  the  eraissivities  of  the  underlying  surface.  For  example,  the 
accuracy  of  water  vapor  retrievals  near  the  surface  is  much  better  in  general 
over  the  oceans  than  over  land.  This  is  due  to  the  greater  contrast  to  water 
vapor  emission  provided  by  the  less  emissive  ocean  surface  background.  For 
this  reason,  the  inclusion  of  data  from  the  SSM/I  microwave  imager  at  spatial 
resolutions  of  from  12.5  to  50  km  can  potentially  improve  the  retrievals  of 
temperature  and  moisture  by  helping  to  characterize  the  emissive  properties  of 
the  fields-of-view  for  the  SSM/T-1  microwave  temperature  and  SSM/T-2  microwave 
moisture  sounders.  These  sensors  have  spatial  resolutions  of  about  180  and  50 
km,  respectively. 

Likewise,  retrievals  are  to  be  performed  for  a  variety  of  meteorological 
parameters  using  data  from  the  SSM/I  sensor  itself.  Again,  over  spatial 
scales  of  up  to  50  km  there  is  considerable  variation  in  the  properties  of 
both  the  atmosphere  and  the  underlying  surface.  As  an  imager,  the  SSM/I  seeks 
to  measure  both  atmospheric  and  surface  parameters.  The  size  of  the  field-of- 
view  would  not  be  a  problem  if  it  was  smaller  than  the  characteristic  scales 
of  the  phenomenon  to  be  measured.  Unfortunately  this  is  not  the  case.  De¬ 
sired  atmospheric  parameters  such  as  cloud  liquid  water,  precipitation,  and 
water  vapor  commonly  vary  on  scales  of  as  little  as  a  few  kilometers. 
Individual  convective  elements,  for  example,  are  typically  a  kilometer  or  so 
in  horizontal  extent.  Likewise  surface  properties  such  as  soil  moisture, 
vegetation,  snow  cover,  and  sea  ice  vary  on  many  scales.  It  must  be  under¬ 
stood  that  retrievals  of  these  quantities  based  on  simulations  which  assume 
uniform  conditions  within  the  SSM/I  sensor  field-of -view  will  be  biased  when 
observing  non-uniform  conditions.  It  is  likely  that  nonuniforra  conditions 
will  be  common. 

Other  sources  of  data  exist  to  help  in  characterizing  the  nonuniformity 
within  the  SSM/I  field-of-view.  Colocated  with  the  SSM/I  aboard  the  DMSP 
spacecraft,  the  Operational.  Linescan  System  (OLS)  provides  both  visible  and 
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infrared  imagery  at  high  spatial  resolution.  Currently  there  are  no  plans  to 
use  these  data  along  with  that  from  the  microwave  imager,  although  there  are 
advantages  from  the  perspectives  of  both  meteorological  analysis  and  SSM/1 
retrievals  (Isaacs  and  Barnes,  1985).  With  much  higher  spatial  resolution, 
the  visible  and  infrared  data  from  the  OLS  imagery  can  be  used  to  provide 
guidelines  on  the  uniformity  of  the  SSM/I  footprints.  In  those  areas  where 
the  contributions  from  the  atmosphere  to  microwave  brightness  temperature  are 
small  (i.e.,  nonprecipitating  situations),  daytime  visible  data  should  be  able 
to  provide  guidelines  on  the  uniformity  of  the  surface  observed  within  an 
SSM/I  f ield-of -view .  This  capability  will  be  degraded  when  clouds  obscure 
portions  of  the  field-of-view. 

In  cases  of  precipitation,  the  availability  of  visible  and  infrared  data 
should  considerably  enhance  the  information  content  of  the  microwave  retrieval 
by  providing  guidelines  to  ascertain  the  spatial  uniformity  of  the  rainfall. 

It  is  well  known  that  a  variety  of  methods  exist  to  monitor  rainfall  which  are 
based  on  the  use  of  visible  and  infrared  data  alone  (see  Barrett  and  Martin, 
1981).  These  methods  differ  from  those  used  with  microwave  data  since  they 
are  based  on  inferences  derived  from  empirical  relationships  between  the 
properties  of  the  clouds  accompanying  the  precipitation  and  the  magnitude  of 
the  rainfall.  The  microwave  methods  result  from  actual  interaction  between 
the  microwave  radiation  and  precipitation  and,  but  for  the  larger  field-of- 
view,  should  potentially  provide  more  quantitative  rainfall  rate  retrievals. 
Nevertheless,  the  visible  and  infrared  data  has  a  spatial  resolution  advantage 
over  the  microwave  data.  Thus  a  combination  of  the  two  data  sources  (i.e., 
microwave  and  vis/IR)  should  be  exploited. 

Based  on  the  considerations  outlined  above,  we  have  undertaken  to  in¬ 
vestigate  the  application  of  visible  and  infrared  imagery  from  the  OLS  to  aid 
in  the  characterization  of  the  uniformity  of  the  field-of-view  of  the  micro- 
wave  imager.  This  investigation  has  focused  on  two  questions:  (1)  can  high 
resolution  visible  and  infrared  imagery  provide  information  on  the  uniformity 
of  the  microwave  field-of-view  either  when  observing  the  surface  or  atmo- 
sphe  ric  phenomenon  such  as  precipitation?,  and  (2)  can  radiative  transfer 
models  of  the  surface  emissivity  and  precipi tat  ion  be  used  along  with  informa¬ 


tion  obtained  from  the  OLS  imagery  on  the  nonuni iormity  of  the  observed  scene 
to  simulate  the  averaging  process  which  occurs  within  the  microwave  field-of- 
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3.4  of  this  reDort 
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In  order  to  evaluate  the  feasibility  of  characterizing  the  uniformity  of 
the  microwave  imager  f ield-of-view  using  visible  and  infrared  imagery  and 
assess  the  resulting  benefits  for  non-uniform  scene  modeling,  the  use  of  real 
data  was  desired.  Since  SSM/I  data  was  not  yet  available,  this  study  relied 
on  archived  Scanning  Multichannel  Microwave  Radiometer  (SMMR)  data  from  the 
Nimbus-7  spacecraft.  The  SMMR  data  was  used  in  conjunction  with  concurrent 
OLS  data  obtained  in  photographic  positive  imagery  format.  The  data  subset 
was  previously  used  to  investigate  the  capabilities  of  combining  these  sensors 
for  the  purpose  of  analyzing  precipitation  situations  for  tactical  forecast 
situations. 

The  SMMR  data  used  in  the  analysis  as  an  analog  for  5SM/1  imagery  were 
obtained  from  the  NASA  Goddard  Spaceflight  Center  in  the  CELL-ALL  format.  The 
CELL-ALL  listings  provide  horizontal  and  vertical  polarized  brightness 
temperatures  for  each  of  the  five  SMMR  frequencies.  The  spatial  resolution 
varies  with  frequency  from  150  km  for  the  6.6  GHz  channel  to  60  km  for  the 
highest  frequency  SMMR  channel  at  37  GHz.  The  37  CHz  data  from  the  SSM/I,  for 
example,  should  be  available  at  a  spatial  resolution  of  about  25  km.  It 
should  be  noted  that  the  data  used  is  considered  to  be  preliminary,  unverified 
data,  i.e.,  final  antenna  corrections  have  not  been  applied,  however,  it 
should  be  more  than  adequate  for  these  purposes . 

Since  digital  data  is  not  readily  available  from  the  operational  DMSP 
OLS,  digital  analysis  has  been  performed  using  simulated  OLS  data.  The 
simulation  process  is  based  on  transformation  of  Landsat  MSS  data  obtained  in 
digital  format  using  a  methodology  described  in  Isaacs  et  al.,  (1986)  and 
briefly  reviewed  below.  The  Landsat  data  has  higher  spatial  and  spectral 
resolution  than  the  DMSP  data.  In  order  to  simulate  DMSP  data,  the  Landsat 
data  is  spatially  averaged  to  the  desired  resolution  of  DMSP  LF  or  LS  data. 
Since  the  four  MSS  spectral  bands  cover  the  wavelength  domain  between  0.4  and 
1.1  pm  spanned  by  the  DMSP  OLS  visible  data,  the  MSS  data  can  also  be 
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spectrally  averaged  to  provide  the  DMSP  bandpass  to  high  accuracy.  An  example 
of  this  transformation  is  illustrated  in  Figure  4-1.  Figure  4-la  is  the 
original  scene  depicted  from  Landsat  MSS  bank  4  (0.4-0. 5  /ira) .  Figure  4-lb  is 
the  same  scene  transformed  to  simulate  DMSP  OLS  LS  data  using  all  four  MSS 
bands  and  spatially  smoothed  to  provide  the  correct  spatial  resolution.  The 
width  of  the  original  Landsat  data  was  chosen  to  be  600  pixels.  This  cor¬ 
responds  to  approximately  50  km  or  the  width  of  the  field-of -view  for  the  37 
GHz  channels  of  the  SSM/I .  As  can  be  seen  from  the  resolution  element  size 
for  the  simulated  OLS  LF  data  in  Figure  4-lb,  the  microwave  field-of -view 
corresponds  to  about  one  hundred  visible  f ields-of -view . 


To  investigate  the  non-uniformity  of  the  microwave  footprint  using 
visible  and  infrared  data,  we  must  assume  that  spectral  signatures  which  are 


manifested  in  the  visible  and  infrared  data  also  result  in  microwave  non¬ 


uniformity.  This  assumption  should  be  quite  reasonable.  Essentially,  for 
visible  data  we  assume  that  those  surface  types  which  differ  in  visible 
reflectance  also  differ  in  microwave  emissivity.  The  gross  categories  of 
soil,  vegetation,  water,  snow,  ice,  etc.  obey  this  simple  rule.  For  the 
infrared  data,  the  assuption  is  that  types  with  differing  infrared  emissivity 
and/or  emission  temperature  also  differ  in  microwave  emissivity.  For  surface 
features,  this  may  not  always  be  true.  For  example,  bare  soil  and  snow  cover 
may  have  the  same  infrared  signature  (i.e.  same  surface  temperature  and 
emissivity) ,  but  have  different  microwave  signatures  due  to  microwave  signa¬ 
tures  of  water  bodies  and  land  surfaces  will  differ  due  to  varying  thermal 
heat  capacities,  while  their  respective  microwave  emissivities  are  quite 
different  due  to  distinctive  permitivities .  For  atmospheric  rather  than 
surface  sources  of  non-uniformity,  such  as  precipitating  cells,  the  correla¬ 
tion  between  visible/infrared  and  microwave  signatures  will  insure  that  the 
adopted  assumption  is  valid.  It  is  understood  that  clear,  non-precipitating 
areas  will  rarely  be  confused  with  cloudy,  precipitating  areas  in  visible  or 


infrared  data. 


To  characterize  the  non-uniformity  of  the  microwave  footprint,  therefore, 
one  must  quantify  the  non-uniformity  of  essentially  a  problem  in  image  proces¬ 
sing.  Classically,  techniques  to  infer  non-uniformity  have  been  referred  to 
as  texture  analysis  methods.  A  number  of  approaches  can  be  used  for  texture 
analysis  including:  (a)  examination  of  the  spatial  power  spectrum  of  an  image 
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through  Fourier  decomposition,  (b)  edge  enhancement,  and  (c)  spatial  co¬ 
herence.  We  have  chosen  to  investigate  the  spatial  coherence  approach  since 
it  has  application  in  the  determination  of  cloud  properties  from  high  resolu¬ 
tion  scanner  data  (Coakley  and  Bretherton,  1982;  Coakley,  1983).  This  data 
analysis  method  might  thus  be  applicable  to  the  determination  of  both  cloud 
and  surface  properties  from  OLS  data. 

Simply  stated  the  approach  quantifies  what  the  eye  observes  when  viewing 
an  image  or  a  portion  thereof.  The  statistic  evaluated  is  the  local  standard 
deviation  (LSD)  of  radiance  (or  gray  shade)  values  for  the  kth  n  x  n  subset  of 
the  field-of-view.  The  LSD  is  calculated  as: 


1  n  n  -  2  1/2 

LSDk  “  '  ~2  *  =  (Ii,j  -  V  1  < 

n  i-l  j-1  ,J 

where  I-  j  is  the  radiance  of  an  individual  pixel  within  the  n  x  n  subset 
array,  and  1^  is  the  men  radiance  of  the  kth  subset  given  simply  by: 


(4.1) 


Xk  “  -  2  2  S  Xi  1 

n  i-1  j-1 


(4.2) 


To  characterize  a  field-of-view  consisting  of  a  total  of  K  individual  n  x  n 
subsets  (i.e.  a  total  of  K  x  n  x  n  pixels),  the  local  standard  deviation 
(LSD^)  and  mean  (1^)  are  calculates  for  each  subset.  A  plot  of  LSD^ 
vs.  1^  will  then  give  the  following  information: 

(a)  If  the  field-of-view  is  uniform,  all  LSD^  values  will  be 

identically  zero,  and  all  I  values  will  be  equal  to  the  constant 
scene  radiance . 


(b)  If  there  are  a  few  distintive  regions  of  constant  radiance,  LSD^ 

values  will  be  generally  small,  and  1^  values  will  assume  those  of 
the  discernable  subregions. 


(c)  For  a  general  case,  LSD  values  and  1^  values  will  be  continously 
variable . 


For  realistic  geophysical  f iel ds - of -view ,  the  idealized  characterizations 
of  the  spatial  coherence  approach  described  above  will  be  modified.  In 
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general,  we  would  expect  to  obtain  clusters  of  points  with  the  (LSD.l) 
domain.  Simply  cluster  analysis  methods  could  then  be  used  to  identify 
significant  clusters  and  define  criteria  for  inclusion  of  individual  pixels. 
Based  on  the  cluser  analysis,  the  microwave  f ield-of -view  can  be  divided  into 
a  small  number  of  each  FOV  assignable  to  each  category  could  be  determined. 
Figure  4-2  illustrates  an  analogous  application  to  the  determination  of  frac¬ 
tional  cloud  cover  using  infrared  scanner  data  from  Coakley  and  Bretherton 
(1982).  The  local  standard  deviation  vs  mean  radiance  (expressed  as  bright¬ 
ness  temperature)  plot  forms  an  arch.  The  cluster  at  the  foot  of  the  arch 
with  highest  temperature  and  lowest  LSD  can  be  identified  with  uniform  sub- 
regions  radiating  at  the  surface  temperature,  while  the  other  foot  at  lower 
temperatures  can  be  identified  with  the  top  temperature.  The  other  pixels 
with  intermediate  temperatures  and  nonzero  local  standard  deviations  in  the 
arch  are  partially  cloudy  fields -of -view . 


UNIFORM  REGION 


HETEROGENEOUS  REGION 


Fig.  4-2  Spatial  Coherence  Approach  for  Characterizing  Field-of -View 

Uniformity:  (a)  uniform  FOV,  (b)  FOV  with  a  few  distinct  regions 
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4 . 3  Modeling  Nonuniform  Fields-of -View 

The  imagery  analysis  techniques  described  in  the  previous  section  pro¬ 
vides  a  basis  to  obtain  information  on  the  non-uniformity  of  the  microwave 
field-of -view  based  on  the  concurrent  use  of  visible  and  infrared  imagery  data 
with  higher  spectral  resolution.  For  example,  based  on  the  texture  analysis 
arguments  presented,  visible  data  can  be  used,  first  to  decide  whether  the 
microwave  field-of -view  is  non-uniform,  and,  if  so,  to  segment  the  FOV  into  a 
reasonable  number  of  individual  subregions  which  have  similar  characteris¬ 
tics.  These  subregions  can  then  be  assigned  the  emissive  or  radiative  charac¬ 
teristics  of  known  surface  or  atmospheric  categories  (i.e.  rain,  soil,  water, 
snow,  etc.),  given  that  sufficient  a  priori  data  is  available. 

In  this  section  we  discuss  how  fields-of -view  which  have  been  determined 
to  be  non-uniform  can  be  modeled.  The  information  assumed  available  through 
the  concurrent  availability  of  high  resolution  visible  and  infrared  data 
includes  the  fraction  of  the  microwave  field-of -view  which  belongs  to  a  given 
type  category,  F^.  It  is  assumed  that  we  can  model  the  emissive  or  radiative 
properties  of  these  categories  using  one  of  the  models  described  in  Section  3 
of  this  report. 

If  it  is  assumed  that  the  ith  microwave  field  of  view  at  frequency  j 

consists  of  N-  x  M.  pixels,  each  with  surface  emissivity  and  temperature  given 
J  J  £ 

by  e-  and  T  ,  respectively,  then  the  area  averaged  brightness  temperature,  Tz, 
J  5  J 

will  be  given  by: 


v- 
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J  J 
Z  Z 
n-1  m-1 


e.nm  T  /(N  .  M  ) 
J  s  j  j 


(4.3) 


For  a  completely  uniform  field-of -view,  the  surface  emissivities  for  each 
pixel  will  be  identical,  i.e.: 


mn 

e  .  -  e  . 

J  J 


(4.4) 


Thus,  for  a  uniform  surface  temperature,  T  ,  the  ith  field-of -view  brightness 
temperature  will  be: 
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If  the  f ield-of-view  Is  characterized  as  nonuniform,  the  pixel  by  pixel 
distribution  given  by  equation  (4.1)  is  too  general  for  practical  use. 

Instead,  a  few  surface  type  categories  such  as  those  listed  in  Table  3-6  of 
Section  3.4  can  be  used.  Assuming  that  L  categories  are  utilized  and  that  the 
fraction  of  the  total  pixels  in  each  category  is  F^,  equation  (4.3)  can  be 
written  as: 


F.  e.  T 
its 


(4.6) 


where : 


2  F  -  1 


(4.7) 


Equation  4.6  assumes  that  the  surface  temperature  is  un.ifoim  throughout  the 
ith  f ield-of-view.  If  this  is  not  the  case  and  surface  temperatures  are 
deducible  on  a  pixel  by  pixel  basis,  then  (4.6)  can  be  generalized  to  give: 


M  N 

L  it 

T  1  -  2  e  2  E  T  “n 

^  i- 1  1  m-1  n-1  S 

where  and  Mjp  are  the  number  of  pixels  in  each  type  category  and: 


M.N.  -  M.N. 
*  *  j  J 


(4.8) 


(4.9) 


The  desired  characteristics  of  the  microwave  f ield-of-view  uniformity  can 
be  obtained  to  some  degree  by  exploiting  the  information  content  of  high 
spatial  resolution  imager  data.  A  f ield-of -view  analysis  technique  based  on 
these  concepts  is  currently  under  investigation  (see  Isaacs  et  al.,  1988c). 


.  -r.  ■  "»*  W"  V*  •*  . 


MW! 


5. 


CONCLUSIONS  AND  RECOMMENDATIONS 


5 . 1  Summary  and  Conclusions 

This  report  has  documented  the  following  research  accomplishments  and 

milestones : 

•  Developed  RADTRAN  based  models  for  the  simulation  of  microwave  and 
millimeter  wave  satellite  sensor  data  applicable  to  the  SSM/I ,  SSM/T-1, 
and  SSM/T-2  (non  scattering  cases,  i.e.  no  precipitation). 

•  Extended  simulation  capabilities  to  treat  multiple  scattering  by  clouds 
and  precipitation  including  effects  of  polarization. 

•  Developed  parameterization  of  precipitation  extinction,  scattering,  and 
angular  scattering  properties  (has  been  adopted  for  use  in  FASCODE!). 

•  Developed  surface  polarization/emissivity  subpackage  capable  of  treating 
ocean,  sea  ice,  snow,  soil  wetness,  and  various  types  of  vegetation. 

•  Evaluated  accuracy  of  183  GHz  water  vapor  retrievals  and  SSM/I 
precipitation  retrievals  using  statistical  approach. 

•  Identified  and  evaluated  the  effect  of  cloud  on  millimeter  wave  moisture 
retrieval . 

•  Developed  and  tested  unified,  non-linear,  physically  based  retrieval 
method  to  obtain  temperature  and  water  vapor  profiles  and  surface 
temperature  and  emissivity  from  the  SSM/T-1,  SSM/T-2,  and  SSM/I 
instrument  suite. 

'  Demonstrated  improvement  of  above  retrievals  over  purely  statistical 
results  for  simulated  data  sets. 

•  Developed  image  processing  based  approaches  to  characterize  microwave 

f ields-of -view  using  high  spatial  resolution  visible/infrared  imager  data 
(add-on  task) . 
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In  this  report  we  have  focused  on  the  retrieval  of  temperature  and  mois¬ 
ture  using  the  DMSP  microwave  temperature  sounder  (SSM/T-1)  and  millimeter 
wave  moisture  sounder  (SSM/T-2) .  It  has  been  demonstrated  that  in  cloud  free 
conditions,  the  retrieval  of  moisture  using  a  purely  statistical  approach  is 
improved  considerably  when  data  from  the  temperature  sounder  is  combined  with 
that  of  the  moisture  sounder.  The  moisture  sounder  data  appears  to  do  little 
to  improve  the  statistical  retrieval  of  temperature.  Furthermore,  and  sig¬ 
nificantly,  it  has  been  demonstrated  that  additional  information  can  be  ob¬ 
tained  from  the  brightness  temperature  data  of  both  instruments  by  adjusting 
the  statistical  retrieval  first  guess  profiles  of  temperature  and  moisture 
using  physically  based  considerations.  The  effect  of  this  optional  physical 
retrieval  step  is  to  increase  the  accuracy  of  inferred  surface  temperatures 
and  moisture  near  the  surface. 

The  physically  based  correction  is  in  effect,  a  single  iteration  of  a 
general  physical  inversion  scheme  such  as  that  proposed  by  Chahine .  The 
method  used  is  to  simulate  sensor  brightness  temperatures  based  on  the  first 
guess  profiles  of  temperature  and  moisture  and  guesses  for  the  surface 
temperature  and  emissivity.  The  simulation  is  performed  using  a  numerical 
scheme  to  integrate  the  equation  of  radiative  transfer.  The  moisture  and 
temperature  profile  first  guesses  are  simply  the  statistical  retrievals 
obtained  from  the  operational  statistically  based  algorithms  for  the  SSM/T-1 
and  that  proposed  for  the  SSM/T-2.  In  our  tests  we  used  an  analogous  "D" 
matrix  approach  based  on  the  work  of  Gaut  et  al.  (1975)  to  provide  this  first 
step.  Once  the  sensor  data  is  simulated  based  on  the  first  guesses,  it  is 
compared  with  the  actual  data.  If  the  differences  are  less  than  an  acceptable 
tolerance  (determined  by  sensor  and  numerical  noises),  the  statistical  first 
guess  is  accepted  as  the  final  retrieval.  If  not,  the  difference  between  the 
simulated  and  actual  brightness  temperatures  (i.e.,  the  residual)  is  used  to 
correct  the  first  guess  profiles.  The  scheme  used  to  perform  this  adjustment 
is  that  discussed  by  Smith  et  al.  (1985).  The  effect  of  this  correction  is  to 
decrease  the  residual,  and  thus  to  make  the  retrieved  profiles  more  consistent 
with  the  observed  brightness  temperatures.  In  order  to  increase  the  computa¬ 
tional  efficiency  of  the  overall  scheme,  we  have  performed  our  calculations 
based  on  a  single  corrective  step. 


In  order  to  optimally  obtain  these  quantities,  we  believe  that  it  is 
necessary  to  consider  a  unified  retrieval  scheme  combining  the  moisture 
sounder  data  with  that  from  the  SSM/T-1  microwave  temperature  sounder,  the 
SSM/I  microwave  imager,  and  the  OLS .  Although  there  are  distinct  advantages 
to  this  type  of  approach,  it  has  not  to  our  knowledge  been  proposed  for 
operational  use.  On  the  contrary,  it  is  planned  that  data  from  the  three 
microwave  mission  sensors  cited  (i.e.,  SSM/T-1,  SSM/T-2,  SSM/1). 

This  report  prescribes  a  unified  retrieval  approach  tailored  to  the  DMSP 
meteorological  sensor  suite.  All  of  the  available  data  is  utilized  in  a 
multispectral  sense.  The  existing  operational  statistical  retrieval  capa¬ 
bility  is  exploited  to  provide  parameter  first  guesses.  First  guess  derived 
brightness  temperature  simulations  are  then  compared  with  the  original  data 
and  adjustments,  based  on  sensor  channel  residuals,  are  made  to  the  first 
guesses  as  required.  A  physically  based  simultaneous  retrieval  provides  the 
appropriate  parameter  adjustments.  Results  of  a  retrieval  simulation  cal¬ 
culation  illustrate  the  improvements  possible  to  the  statistical  first  guesses 
by  utilizing  all  sensor  data  multispectrally  and  to  the  first  guess  parameters 
by  employing  the  physical  adjustment  step.  Recognizing  the  potential  effect 
of  clouds  on  the  millimeter  wave  data,  a  method  is  outlined  to  introduce 
required  first  guess  cloud  information  by  image  processing  visible/infrared 
imager  data.  This  procedure  provides  the  opportunity  to  characterize  the 
uniformity  (both  atmospheric  and  surface)  of  the  relatively  large  microwave 
FOV.  This  provides  potential  insights  into  both  cloud  and  surface  type 
classification  with  implications  for  determining  their  emissivities . 

While  the  retrieval  approach  outlined  has  been  formulated  as  a  stand 
alone  system,  its  most  important  application  may  be  to  provide  data  for  use  in 
a  numerical  weather  prediction  model.  In  this  context  it  is  noted  that 
modifications  to  the  procedure  outlined  are  desirable.  Operationally,  for 
example,  it  may  be  advantageous  to  obtain  some  first  guess  elements  such  as 
the  temperature  and  moisture  profiles  from  model  predicted  fields.  Further¬ 
more,  in  the  context  of  a  prediction  model  application,  forecast  error  co- 
variances  may  be  used  as  constraints  on  the  adjustment  process  (J .  Eyre, 
personal  communication) .  These  changes  are  not  inconsistent  with  the 
retrieval  described  and  should  in  fact  improve  its  effectiveness.  Further¬ 
more,  the  method  provides  estimates  of  each  retrieval's  accuracy. 
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5 . 2  Discussion  of  Future  Work 

The  work  accomplished  during  the  tenure  of  this  contract  has  provided  a 
number  of  useful  tools  and  opened  a  few  new  avenues  of  investigation.  The 
testing  of  the  unified  retrieval  algorithm  has  not  been  comprehensive  in  the 
sense  that  only  a  few  climatological  scenarios  were  investigated.  This  should 
be  improved  by  undertaking  more  representative  tests.  An  important  considera¬ 
tion  along  these  lines  is  the  development  and  testing  of  the  image  processing 
based  cloud  analysis  proposed  to  be  applied  to  OLS  data  to  provide  cloud  in¬ 
formation  in  the  unified  retrieval.  All  of  the  testing  reported  here  was  done 
essentially  for  clear  fields-of-view.  Therefore  this  remains  an  important 
area  for  future  research.  This  work  should  be  coordinated  with  ongoing 
studies  related  to  the  operational  nephanalysis  in  order  to  coordinate  the 
type  of  cloud  information  desired  and  the  techniques  applied  to  obtain  them. 
Ultimately,  the  nephanalysis  should  be  a  by-product  of  the  unified  retrieval 
scheme.  The  use  of  supporting  data  sets  from  the  non-OLS  sensors  should 
improve  the  product  currently  available. 

Many  of  the  models  developed  during  the  course  of  this  research  such  as 
the  surface  models,  the  precipitation  optical  properties  models  and  the 
retrieval  code  could  be  usefully  combined  into  an  overall  microwave  simula¬ 
tion/retrieval  code.  We  have  not  accomplished  this  during  the  course  of  this 
work,  relying  instead  on  segmenting  the  calculation  and  providing  necessary 
input  data  to  subsequent  code  elements.  It  would  be  useful  to  other  users  in 
the  DoD  modeling  community  to  perform  this  code  consolidation  at  some  time  in 
the  future . 

Finally,  there  are  a  number  of  subtleties  related  to  the  formulation  of 
the  retrieval  code  which  could  be  improved.  The  unified  approach  would  be 
enhanced  if  these  were  pursued  in  the  future  in  a  systemmatic  way. 
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